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PREFACE 

In this working group we have investigated the propects for Higgs boson searches at the Tevatron and 
LHC and, in particular, the potential of these colliders to determine the Higgs properties once these 
particles have been found. The analyses were done in the framework of the Standard Model (SM) and 
its supersymmetric extensions as the minimal (MSSM) and next-to-minimal (NMSSM) supersymmetric 
extensions. The work for the discovery potential of the LHC mainly concentrated on the difficult regions 
of previous analyses as those which are plagued by invisible Higgs decays and Higgs decays into su- 
persymmetric particles. Moreover, the additional signatures provided by the weak vector-boson fusion 
process (WBF) have been addressed and found to confirm the results of previous analyses. A major 
experimental effort has been put onto charged Higgs boson analyses. The final outcome was a significant 
improvement of the discovery potential at the Tevatron and LHC than previous analyses suggested. 

For an accurate determination of Higgs boson couplings, the theoretical predictions for the signal 
and background processes have to be improved. A lot of progress has been made during and after this 
workshop for the gluon-fusion gg — > H + (0, 1,2 jets) and the associated ttH production process. 
A thorough study of the present theoretical uncertainties of signal and background processes has been 
initialized, culminating in a list of open theoretical problems. A problem of major experimental interest 
is the proper treatment of processes involving bottom quark densities, which is crucial for some important 
signal and background processes. Further theoretical improvements have been achieved for the MSSM 
Higgs boson masses and Higgs bosons in the NMSSM. 

This report summarizes our work. The first part deals with theoretical developments for the sig- 
nal and background processes. The second part gives an overview of the present status of Higgs boson 
searches at the Tevatron. The third part analyzes invisible Higgs boson decays at the LHC and the forth 
part the Higgs boson search in the WBF channel. Part 5 summarizes the progress that has been achieved 
for A/H — > t + t~ decays in the MSSM. In part 6 the status of the Higgs boson search in ttH production 
is presented. Finally, part 7 describes the charged Higgs boson analyses in detail. 

Acknowledgements. 
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A. Theoretical Developments 

S. Balatenychev, G. Belanger, F. Boudjema, A. Cottrant, M. Carena, S. Catani, V. Del Duca, D. de Florian, 
M. Frank, R. Godbole, M. Grazzini, S. Heinemeyer, W. Hollik, V. Ilyn, W. Kilgore, R. Lafaye, S. Moretti, 
C. Oleari, A. Pukhov, D. Rainwater, DP Roy, C. Schmidt, A. Semenov, M. Spira, C. Wagner, G. Weiglein 
and D. Zeppenfeld 



Abstract 

New theoretical progress in Higgs boson production and background processes 
at hadron colliders and the relations between the MSSM Higgs boson masses 
is discussed. In this context new proposals for benchmark points in the MSSM 
are presented. Additional emphasis is put on theoretical issues of invisible 
SUSY Higgs decays and multiple Higgs boson production within the NMSSM. 



1. Higgs boson production at hadron colliders: signal and background processes^ 
1.1 Introduction 

The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersymmetric extensions. 
Thus, the search for Higgs bosons is one of the most important endeavors at future high-energy experi- 
ments. In the SM one Higgs doublet has to be introduced in order to break the electroweak symmetry, 
leading to the existence of one elementary Higgs boson, H [JXJ| _ The scalar sector of the SM is uniquely 
fixed by the vacuum expectation value v of the Higgs doublet and the mass m# of the physical Higgs 
boson [gj. The negative direct search for the Higgsstrahlung process e + e~ — > ZH at the LEP2 collider 
poses a lower bound of 114.1 GeV on the SM Higgs mass while triviality arguments force the 

Higgs mass to be smaller than ~ 1 TeV [|]]. 

Since the minimal supersymmetric extension of the Standard Model (MSSM) requires the in- 
troduction of two Higgs doublets in order to preserve supersymmetry, there are five elementary Higgs 
particles, two CP-even (h, H), one CP-odd (A) and two charged ones (H^). At lowest order all cou- 
plings and masses of the MSSM Higgs sector are fixed by two independent input parameters, which 
are generally chosen as tan/3 = vi/v\, the ratio of the two vacuum expectation values v% t 2, and the 
pseudoscalar Higgs-boson mass mA- At LO the light scalar Higgs mass has to be smaller than 
the Z-boson mass mz- Including the one-loop and dominant two-loop corrections the upper bound 
is increased to < 135 GeV The negative direct searches for the Higgsstrahlung processes 

e + e~ — ► Zh, ZH and the associated production e + e~ — ► Ah, AH yield lower bounds of m^H > 91.0 
GeV and itla > 91.9 GeV. The range 0.5 < tan (3 < 2.4 in the MSSM is excluded by the Higgs searches 
at the LEP2 experiments 

The intermediate mass range, mjf < 196 GeV at 95% CL, is also favored by a SM analysis of 
electroweak precision data [|3],|]]. In this contribution we will therefore concentrate on searches and 
measurements for run < 200 GeV. The Tevatron has a good chance to find evidence for such a Higgs 
boson, provided that sufficient integrated luminosity can be accumulated [10]. The Higgs boson, if it 
exists, can certainly be seen at the LHC, and the LHC can provide measurements of the Higgs boson 
mass at the 10~ 3 level [11, 12], and measurements of Higgs boson couplings at the 5 to 10% level [13]. 
Both tasks, discovery and measurement of Higgs properties, require accurate theoretical predictions of 
cross sections at the LHC, but these requirements become particularly demanding for accurate coupling 
measurements. 

In this contribution we review the present status of QCD calculations of signal and background 
cross sections encountered in Higgs physics at hadron colliders. Desired accuracy levels can be esti- 
mated by comparing to the statistical errors in the determination of signal cross sections at the LHC. 
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For processes like H — > 77, where a very narrow mass peak will be observed, backgrounds can be 
accurately determined directly from data. For other decay channels, like H — > bb or H — > rr, mass res- 
olutions of order 10% require modest interpolation from sidebands, for which reliable QCD calculations 
are needed. Most demanding are channels like H — > W + W~ — ► l + l~p" T , for which broad transverse 
mass peaks reduce Higgs observation to, essentially, a counting experiment. Consequently, requirements 
on theory predictions vary significantly between channels. In the following we discuss production and 
decay channels in turn and focus on theory requirements for the prediction of signal and background 
cross sections. Because our main interest is in coupling measurements, we will not consider diffractive 
channels in the following, which are model-dependent and have large rate uncertainties [ |l4| ] ; potentially, 
they might contribute to Higgs discovery if, indeed, cross sections are sufficiently large. 



1.2 Gluon fusion 



The gluon fusion mechanism gg —> qb provides the dominant production mechanism of Higgs bosons 
at the LHC in the entire relevant mass range up to about 1 TeV in the SM and for small and moderate 



values of tan (3 in the MSSM [15]. At the Tevatron this process plays the relevant role for Higgs masses 



between about 130 GeV and about 190 GeV [|10p. The gluon fusion process is mediated by heavy quark 
triangle loops and, in the case of supersymmetric theories, by squark loops in addition, if the squark 
masses are smaller than about 400 GeV [16], see Fig. [[]. 



g onr?nnn 



9 ^QJLQJLQJ 




Fig. 1 : Typical diagram contributing to gg — > <j> at lowest order. 

In the past the full two-loop QCD corrections have been determined. They increase the production 
cross sections by 10-80% Q17[|, thus leading to a significant change of the theoretical predictions. Very 



recently, Harlander and Kilgore have finished the full NNLO calculation, in the heavy top quark limit []18 



19]. This limit has been demonstrated to approximate the full massive K factor at NLO within 10% 



for the SM Higgs boson in the entire mass range up to 1 TeV [2C]. Thus, a similar situation can be 
expected at NNLO. The reason for the quality of this approximation is that the QCD collections to the 
gluon fusion mechanism are dominated by soft gluon effects, which do not resolve the one-loop Higgs 
coupling to gluons. Fig. || shows the resulting iT-factors at the LHC and the scale variation of the K- 
factor. The calculation stabilizes at NNLO, with remaining scale variations at the 10 to 15% level. These 
uncertainties are comparable to the experimental errors which can be achieved with 200 fb _1 of data at 
the LHC, see solid lines in Fig. ||. The full NNLO results confirm earlier estimates which were obtained 



in the frame work of soft gluon resummation Q20| ] and soft approximations Q21| , |22|] of the full three-loop 



result. The full soft gluon resummation has been performed in Ref. []23[]. The resummation effects 
enhance the NNLO result by about 10% thus signaling a perturbative stabilization of the theoretical 
prediction for the gluon-fusion cross section. 

In supersymmetric theories the gluon fusion cross sections for the heavy Higgs, H, and, for small 
rriA, also for the light Higgs, h, may be dominated by bottom quark loops for large values of tan > 10 
so that the heavy top quark limit is not applicable. This can be clearly seen in the NLO results, which 
show a decrease of the K factor down to about 1.1 for large tan j3 [|l7|]. This decrease originates from 
an interplay between the large positive soft gluon effects and large negative double logarithms of the 
ratio between the Higgs and bottom masses. In addition, the shape of the pr distribution of the Higgs 
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Fig. 2: Scale dependence of the K -factor at the LHC. Lower curves for each pair are for fin — luiu, = mn/2, upper 
curves are for (1r = run /2, (If = 2mn- The K-factor is computed with respect to the LO cross section at fin — fiF = mn- 
From fle/ /|7J/. 

boson may be altered; if the bottom loop is dominant, the pt spectrum becomes softer than in the case 
of top-loop dominance. These effects lead to some model dependence of predicted cross sections. 

Let us now turn to a discussion of backgrounds for individual decay modes. 

(i) <p — > 77: At the LHC the SM Higgs boson can be found in the mass range up to about 150 GeV by 
means of the rare photonic decay mode H — » 77 [11, 12 1. The dominant Higgs decays H — > bb, t + t~ 
are overwhelmed by large QCD backgrounds in inclusive searches. The QCD 77 background is known 
at NLO, including all relevant fragmentation effects. The present status is contained in the program 
DIPHOX [|25|]. The loop mediated process gg — > 77 contributes about 50% to the 77 background and 
has been calculated at NLO very recently fl26|]. However, a numerical analysis of the two-loop result is 
still missing. 

Once the experiment is performed, the diphoton background can be determined precisely from the 
data, by a measurement of da/dm^ on both sides of the resonance peak. The NLO calculations are use- 
ful, nevertheless, for an accurate prediction of expected accuracies and for a quantitative understanding 
of detector performance. 

( ii) H — > W + W~ : This mode is very important for Higgs masses above VF-pair but below Z-pair thresh- 
old, where B(H — » WW) is close to 100%. In order to suppress the ti — > bbW + W~ background for 
W + W~ final states, a jet veto is crucial. However, gluon fusion receives sizeable contributions from 
real gluon bremsstrahlung at NLO, which will also be affected by the jet veto. These effects have re- 
cently been analyzed in Ref. [27], in the soft approximation to the full NNLO calculation. A veto of 



additional jets with pTj > 15 GeV, as e.g. envisioned by ATLAS Q12[], reduces the NNLO iv~-factor to 
about K = 0.80, i.e. one loses more than 60% of signal events. In addition the scale dependence of 
the cross section starts to grow with such stringent veto criteria. These effects need to be modeled with 
a NLO Monte Carlo program for H + jet production in order to reach a reliable quantitative result for 
the signal rate. Since stop and sbottom loops are sizeable in supersymmetric theories for squark masses 

2 It should be noted that for this strong cut in prj the NNLO result may be plagued by large logarithms of this cut, which 
have to be resummed, see IE3|]. 
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Fig. 3: Expected relative error on the determination of B a for various Higgs search channels at the LHC with 200 fb~ of 
data Solid lines are for inclusive Higgs production channels which are dominated by gluon fusion. Expectations for 

weak boson fusion are given by the dashed lines. Dotted lines are for tiH production with H —* bb /^/ (black) and H — > 
W + W~ /^/ (red) and assume 300 ft)- 1 of data. 

below about 400 GeV, their inclusion is important in these investigations. 

From the perspective of background calculations, H — ► WW is the most challenging channel. 
Backgrounds are of the order of the signal rate or larger, which requires a 5% determination or better for 
the dominant background cross sections in order to match the statistical power of LHC experiments. In 
fact, the large errors at mjj < 150 GeV depicted in Fig. || (gg —> H —> WW curve) are dominated by an 
assumed 5% background uncertainty. Clearly, such small errors cannot be achieved by NLO calculations 
alone, but require input from LHC data. Because of two missing neutrinos in the W + W~ — > l + l~fi T 
final state, the Higgs mass cannot be reconstructed directly. Rather, only wide {l + l~\$ T ) transverse 
mass distributions can be measured, which do not permit straightforward sideband measurements of the 
backgrounds. Instead one needs to measure the normalization of the backgrounds in signal poor regions 
and then extrapolate these, with the help of differential cross sections predicted in perturbative QCD, to 
the signal region. The theory problem is the uncertainty in the shape of the distributions used for the 
extrapolation, which will depend on an appropriate choice of the "signal poor region". No analysis of 
the concomitant uncertainties, at LO or NLO QCD, is available to date. 

After the jet veto discussed above, the dominant background processes are pp — ► W + W~ and 
(off-shell) tt production [11, 12]. W + W~ production is known at NLO [28] and available in terms 
of partem level Monte Carlo programs. In addition, a full NLO calculation including spin correlations 
of the leptonic W, Z decays, in the narrow width approximation, is available [29]. For Higgs boson 
masses below the W + W~(ZZ) threshold, decays into WW*(ZZ*) are important [15, 3^]. Since hadron 
colliders will be sensitive to these off-shell tails, too, the backgrounds from VV* production become 
relevant. There is no NLO calculation of VV* background processes available so far, so that it is not 
clear if NLO effects will be significant in the tails of distributions needed for the Higgs search in these 
cases. Moreover, for WW* production the inclusion of spin correlations among the final state leptons is 



mandatory \pl\j. 

Top quark backgrounds arise from top-pair and tWb production. Recently, a new theoretical 
analysis of pp — > fM£M has become available including full lepton correlations and off-shell effects of 



the final state top quarks arising from the non-zero top decay width [32]. This calculation automatically 
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includes pp — ► tbW and those contributions to pp — ► bbW + W~ , which are gauge-related to i&VF 
couplings and describes the relevant tails for the Higgs search at LO. It is now necessary to investigate 
the theoretical uncertainties of this background. A NLO calculation of off-shell top-pair production may 
well be needed to reach the required 5% accuracy for extrapolation to the Higgs search region. 

Other important reducible backgrounds are the Wti, Zti, Wbb and Zbb production processes. 
While Vtt (V = W, Z) production is only known at LO, the associated vector boson production with bb 
pairs is known at NLO including a soft gluon resummation [M]. Thus Vbb production can be considered 
as reliable from the theoretical point of view, while a full NLO calculation for Vtt production is highly 
desirable, since top mass effects will play a significant role. In addition, the background from gb — ► 
tH~ , gb — > iH + has to be taken into account within the MSSM framework. The full LO matrix elements 
are included in the ISAJET Monte Carlo program, which can easily be used for experimental analyses. 

(Hi) H — > ZZ — > 4£ ± : A sharp Higgs peak can be observed in the four lepton invariant mass distribu- 
tion. Hence, the ZZ — > 4£ ± backgrounds are directly measurable in the sidebands and can safely be 
interpolated to the signal region. 



1.3 qq — > qqH 

In the SM the WW, ZZ fusion processes qq — ► qqV*V* — > qqH play a significant role at the LHC for 
the entire Higgs mass range up to 1 TeV. We refer to them as weak boson fusion (WBF). The WBF cross 
section becomes comparable to the gluon fusion cross section for Higgs masses beyond ~ 600 GeV [15] 
and is sizable, of order 20% of a(gg — ► H), also in the intermediate mass region. The energetic quark jets 
in the forward and backward directions allow for additional cuts to suppress the background processes to 
WBF. The NLO QCD corrections can be expressed in terms of the conventional corrections to the DIS 
structure functions, since there is no color exchange between the two quark lines at LO and NLO. NLO 
corrections increase the production cross section by about 10% and are thus small and under theoretical 
control [34, 35]. These small theory uncertainties make WBF a very promising tool for precise coupling 
measurements. However, additional studies are needed to assess the theoretical uncertainties associated 
with a central jet veto. This veto enhances the color singlet exchange of the signal over color octet 
exchange QCD backgrounds [|36|-|3"9|]. 

In the MSSM, first parton level analyses show that it should be possible to cover the full MSSM 
parameter range by looking for the light Higgs decay h — > t + t~ (for vtla ^ 150 GeV) and/or the heavy 
Higgs H — > r + T~resonance (for a relatively small vtla) in the vector-boson fusion processes pOj]. 
Although these two production processes are suppressed with respect to the SM cross section, their sum 
is of SM strength. 

For the extraction of Higgs couplings it is important to distinguish between WBF and gluon fusion 
processes which lead to H + jj final states. With typical WBF cuts, including a central jet veto, gluon 
fusion contributions are expected at order 10% of the WBF cross section, i.e. the contamination is 
modest [41, 42]. The gluon fusion processes are mediated by heavy top and bottom quark loops, in 
analogy to the LO gluon fusion diagram of Fig. |l[ The full massive cross section for H + jj production 
via gluon fusion has been obtained only recently [41], while former analyses were performed in the 
heavy top quark limit [43]. Since stop and sbottom loops yield a sizeable contribution to the inclusive 
gluon fusion cross section, a similar feature is expected for H + jj production. Thus, it is important to 
compute the effects of stop and sbottom loops in H + jj gluon fusion processes, which has not been 
done so far. 



(i) H — * 77: Parton level analyses show that H —> 77 decays in WBF Higgs production can be isolated 
with signal to background ratios of order one Q36| ] and with statistical errors of about 15%, for 200 ft- 1 of 
data (see Fig. |3|). Like for the inclusive H — > 77 search, background levels can be precisely determined 
from a sideband analysis of the data. Prior to data taking, however, full detector simulations are needed 



9 



to confirm the partem level results and improve on the search strategies. 

Improved background calculations are desirable as well. In particular, the pp — ► 77 jj background 




fVWW 7 



^VWW 7 



Fig. 4: Typical diagrams contributing to gg — > 77J j af lowest order. 



(ii) H — ► r^r : The observation of if — > rr decays in WBF will provide crucial information on Higgs 



couplings to fermions Ql 3p and this channel alone guarantees Higgs observation within the MSSM [ ]40| ] 
and may be an important discovery channel at low pseudoscalar mass, m^. Recent detector simula- 
tions [42] confirm partem level results [37] for the observability of this channel. (See Fig. || for parton 



level estimates of statistical errors.) 

The r + T~-invariant mass can be reconstructed at the LHC with a resolution of order 10%. This 
is possible in the qq — ► qqH mode because of the large transverse momentum of the Higgs. In turn 
this means a sideband analysis can be used, in principle, to directly measure backgrounds. The most 
important of these backgrounds is QCD Zjj production (from QCD corrections to Drell-Yan) or elec- 
troweak Zjj production via WBF [[37|]. The (virtual) Z (or photon) then decays into a t + t~ pair. These 
Zjj backgrounds, with their highly nontrivial shape around m TT m mz, can be precisely determined be 
observing Z — > e + e~, fi + 'yT events in identical phase space regions. Theoretically the QCD Zjj back- 



ground is under control also, after the recent calculation of the full NLO corrections [|44|]. For the r + r~ 
backgrounds the inclusion of r polarization effects is important in order to obtain reliable tau-decay dis- 
tributions which discriminate between signal processes (h, H — > r + r~) and backgrounds. This can be 
achieved by linking the TAUOLA program [^] to existing Monte Carlo programs. 

(in) H —¥ WW — > £ + £~^> T : The most challenging WBF channel is H — > decay which does 

not allow for direct Higgs mass reconstruction and, hence, precludes a simple sideband determination 
of backgrounds. The important backgrounds [^,|3^] involve (virtual) W pairs, namely top decays in 
it+jets production, and QCD and electroweak WW jj production. QCD and EW rrjj production are 



subdominant after cuts, they are known at NLO [44], and they can be determined directly, in phase space 



regions for jets which are identical to the signal region and with high statistics, by studying e + e or 
H + fi~ pairs instead of t + t~. 

Demands on QCD calculations can be estimated by comparing the effects of systematic back- 
ground errors on the measurement of the signal rate with statistical errors achievable at the LHC with 
200 ftr 1 of data. Results are shown in Fig. || for an assumed 10% error on a(tt+jets), a 50% error on 
the QCD WWjj rate, and a 30% error on the electroweak WWjj rate. The latter two should be achiev- 
able from a LO extrapolation from signal poor to signal rich regions of phase space. A 10% error of 
a(ti+jets), on the other hand, may require a NLO calculation, in particular of the on-shell ti+ 1 jet cross 
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Fig. 5: Contributions of background systematic errors Aa to a measurement of an = o~B(H — -» tWK) in WBF. Shown, from 
bottom to top, are the effects of a 10% uncertainty of the rrjj rate (dotted line), a 50% error on the QCD WWjj rate (blue 
dash-dotted), a 30% error on the electroweak WWjj rate (green dash-dotted), and a 10% error on a(ti+jets) (red dashes). 
The long-dashed line adds these errors in quadrature. For comparison, the solid line shows the expected statistical error for 
200 fb~ . The vertical line at 145 GeV separates analyses optimized for small /|?9|/ and large /pq/ Higgs masses. 



section which dominates the tt background. Off-shell effects have recently been studied at LO Q32| ] and 
a 0(20%) increase of the background was found, which, presumably, is small enough to permit the in- 
clusion of off-shell effects at LO only. However, a dedicated study is needed to devise optimal techniques 
for a reliable background determination for H — ► WW searches in WBF, for all major backgrounds. 

(iv) Jet veto and Jet Tagging: Background suppression in the WBF channels relies on double forward jet 
tagging to identify the scattered quark jets of the qq — ► qqH signal and it employs a veto of relatively 
soft central jets (typically of px > 20 GeV) to exploit the different gluon radiation patterns and QCD 
scales of i-channel color singlet versus color octet exchange. Transverse momenta of these tagging or 
veto jets are relatively small for fixed order perturbative calculations of hard processes at the LHC. Thus, 
dedicated studies will be needed to assess the applicability of NLO QCD for the modeling of tagging 
jets in WBF and for the efficiency of a central jet veto in the Higgs signal. First such studies have been 
performed in the past at LO, for Wjj or Zjj events [|46|]. While NLO Monte Carlos for QCD Vjj 



production are now available [}44|,|47p, the corresponding NLO determination of electroweak Vjj cross 
sections would be highly desirable. This would allow a comparison of calculated and measured veto 
efficiencies in a WBF process. These efficiencies must be known at the few percent level for the signal 
in order to extract Higgs couplings without loss of precision. 

At present, the veto efficiencies for signal and background processes are the most uncertain aspect 
of WBF Higgs production at the LHC. Any improvement in their understanding, from QCD calculations, 
from improved Monte Carlo tools, or from hadron collider data would be very valuable. 
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1.4 tt(j) production 

SM Higgs boson production in association with ti pairs plays a significant role at the LHC for Higgs 
masses below about 130 GeV, since this production mechanism makes the observation of H — > bb pos- 
sible [HJ, 12|, ^4, 48, 49]. The decay H — > 77 is potentially visible in this channel at high integrated 
luminosity. For Higgs masses above about 130 GeV, the decay H — > W + W~ can be observed [50]. itiJ 
production could conceivably be used to determine the top Yukawa coupling directly from the cross sec- 
tion, but this requires either assumptions about the branching ratio for H — > 66, which are not justified in 
extensions of the SM, or observability of decay to either 77 or W + W~ . Recently, the NLO QCD correc- 
tions have become available. They decrease the cross section at the Tevatron by about 20-30% [51, 52], 
while they increase the signal rate at the LHC by about 20-40% [51]. The scale dependence of the pro- 
duction cross section is significantly reduced, to a level of about 15%, which can be considered as an 
estimate of the theoretical uncertainty. Thus, the signal rate is under proper theoretical control now. In 
the MSSM, tth production with h — ► 77, bb is important at the LHC in the decoupling regime, where 
the light scalar h behaves as the SM Higgs boson [|II|,|T2, 24, 48, 4^J. Thus, the SM results can also be 
used for tth production in this regime. 



(i) ttcf) — ► tibb: The major backgrounds to the cj> — > bb signal in associated titfi production come from 
ttjj and tibb production, where in the first case the jets may be misidentified as 6 jets. A full LO calcula- 
tion is available for these backgrounds and will be included in the conventional Monte Carlo programs. 
However, an analysis of the theoretical uncertainties is still missing. A first step can be made by study- 
ing the scale dependence at LO in order to investigate the effects on the total normalization and the 
event shapes. But for a more sophisticated picture a full NLO calculation is highly desirable. A second 
question is whether these backgrounds can be measured in the experiments off the Higgs resonance and 
extrapolated to the signal region. 



(ii) ttxf) — ► ££77: The W77 final states develop a narrow resonance in the invariant 77 mass distribution, 
which enables a measurement of the ti'yj background directly from the sidebands. 



(Hi) tt(f> — ► ttW + W~: This channel does not allow reconstruction of the Higgs. Instead, it relies on a 
counting experiment of multiplepton final states where the background is of approximately the same size 
as the signal [^0|]. The principal backgrounds are tiWjj and ti£ + £~ (jj), with minor backgrounds of 
ttW + W~ and titi. For the 31 channel, the largest background is ti£ + £~ where one lepton is lost. It 
is possible that this rate could be measured directly for the lepton pair at the Z pole and the result ex- 
trapolated to the signal region of phase space. However, for tiVjj backgrounds the QCD uncertainties 
become large and unknown, due to the presence of two additional soft jets in the event. Further inves- 
tigation of these backgrounds is essential, and will probably require comparison with data, which is not 
expected to be trivial. 



1.5 bb<p production 

In supersymmetric theories bb<p production becomes the dominant Higgs boson production mechanism 
for large values of tan/3 [15], where the bottom Yukawa coupling is strongly enhanced. In contrast to 
titfi production, however, this process develops potentially large logarithms, logm^/m^, in the high- 
energy limit due to the smallness of the bottom quark mass, which are related to the development of b 
densities in the initial state. They can be resummed by evolving the b densities according to the Altarelli- 
Parisi equations and introducing them in the production process [53]. The introduction of conventional 
b densities requires an approximation of the kinematics of the hard process, i.e. the initial b quarks are 
assumed to be massless, have negligible transverse momentum and travel predominantly in forward and 
backward direction. These approximations can be tested in the full gg — » bbcj) process. At the Tevatron it 
turns out that they are not valid so that the effective cross section for bb —* cf> has to be considered as an 
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overestimate of the resummed result. An improvement of this resummation requires an approach which 
describes the kinematics of the hard process in a better way. Moreover, since the experimental analyses 
require 3 or 4 b tags [10-[l2|], the spectator b quarks need to have a sizeable transverse momentum of 
at least 15-20 GeV. Thus a resummation of a different type of potentially arising logarithms, namely 
log m 2 ^/ (m 2 + p 2 min ) is necessary. This can be achieved by the introduction of e.g. unintegrated partem 
densities [ p4| ] or an extension of the available resummation techniques. 

As a first step, however, we have to investigate if the energy of the Tevatron and LHC is sufficiently 
large to develop the factorization of bottom densities. This factorization requires that the transverse mo- 
mentum distribution of the (anti)bottom quark scales like da / dpxb <x PTb I (to 2 + Prb) f° r transverse 
momenta up to the factorization scale of the (anti)bottom density. The transverse momentum distribu- 
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Fig. 6: Transverse momentum distributions of the bottom quark in bbH production for two Higgs masses at the LHC. We 
have adopted CTEQ5M1 parton densities and a bottom mass of nib = 4.62 GeV. The solid lines show the full LO result from 
qq, gg bbH and the dashed lines the factorized collinear part, which is absorbed in the bottom parton density. The upper 
curves are divided by the factor prb / (mf +Prb) of the asymptotic behavior, which is required by factorizing bottom densities. 



tions at the LHC are shown in Fig. ^, for two different Higgs masses. The solid curves show the full 
distributions of the qq, gg — > bbcb processes, while the dashed lines exhibit the factorized collinear part, 
which is absorbed in the bottom density. For a proper factorization, these pairs of curves have to co- 
incide approximately up to transverse momenta of the order of the factorization scale, which is usually 
chosen to be ftp = C(mjf). It is clearly visible that there are sizeable differences between the full result 
and the factorized part, which originate from sizeable bottom mass and phase space effects, that are not 
accounted for by an active bottom parton density. Moreover, the full result falls quickly below the ap- 
proximate factorized part for transverse momenta of the order of run /10, which is much smaller than the 
usual factorization scale used for the bottom densities. We conclude from these plots that bb<p production 
at the LHC develops sizeable bottom mass effects, so that the use of bottom densities in the process 
bb — > <p ma y lead to an overestimate of the correct theoretical result due to too crude approximations in 
the kinematics of the hard process. The full NLO calculation of the gg — » bbcb wm yield much more 
insight into this problem, since the large logarithms related to the evolution of bottom densities have to 
appear in the NLO corrections, if the picture of active bottom quarks in the proton is correct. 
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1.6 ZH, WH production 

Higgsstrahlung in qq — ► WH, ZH plays a crucial role for the Higgs search at the Tevatron, while it 
is only marginal at the LHC. At the Tevatron it provides the relevant production mechanism for Higgs 



masses below about 130 GeV, where H — > bb decays are dominant [10]. The NLO QCD corrections have 
been analyzed in the past. They are identical to the QCD corrections to the Drell-Yan processes qq — ► 
W,Z, if the LO matrix elements are replaced accordingly. QCD corrections increase the production 



cross sections by about 30^10% [35, 55] 



The most important backgrounds at the Tevatron are Wjj and in particular Wbb production. Both 
are known at NLO and are contained in a NLO Monte Carlo program p7|]. The same applies also to the 



Zjj and in particular Zbb backgrounds [33, 44 1. In addition, the tt background is relevant 



1.7 Conclusions 

Considerable progress has been made recently in improving QCD calculations for Higgs signal and 
background cross sections at hadron colliders. Noteworthy examples are the NNLO corrections to the 



gluon fusion cross section [19], the QCD Zjj cross section at NLO [44] and the determination of full 



finite top and W width corrections to tt and ttj production at LO [32]. These improvements are crucial 
for precise coupling determinations of the Higgs boson. 

Much additional work is needed to match the statistical power of the LHC. Largely, QCD sys- 
tematic errors for coupling measurements have not been analyzed yet. Additional NLO tools need 
to be provided as well, and these include NLO corrections to tt production with finite width effects 
and ttj production at zero top width. A better understanding of central jet veto efficiencies is crucial 
for the study of WBF channels. These are a few examples where theoretical work is needed. Many 
more have been highlighted in this review. Higgs physics at the LHC remains a very rich field for 
phenomenology. 

2. Direct Higgs production and jet veto[] 

Direct Higgs production through gluon-gluon fusion, followed by the decay H — ► W*W* , Z*Z* is a 
relevant channel to discover a Higgs boson with mass 140 ^ Mjj 190 GeV both at the Tevatron and 
at the LHC. In particular, the decay mode T^*W* — > l + l~vv is quite important [1C-12| 31], since it is 



cleaner than W*W* — > Ivjj, and the decay rate H — > W*W* is higher than H — ► Z*Z* by about one 
order of magnitude. 

An important background for the direct Higgs signal H — > W*W* — > l^l^vv is tt production 
(tW production is also important at the LHC), where t — ► lub, thus leading to b jets with high pr 
in the final state. If the b quarks are not identified, a veto cut on the transverse momenta of the jets 
accompanying the final-state leptons can be applied to enhance the signal/background ratio. Imposing a 



jet veto turns out to be essential, both at the Tevatron [10,56] and at the LHC [11,12,31], to cut the hard 
b jets arising from this background process. 

Here we study the effect of a jet veto on direct Higgs production. More details can be found in 



Ref. [27]. The events that pass the veto selection are those with pj, < p^ to , where pj, is the transverse 
momentum of any final-state jets, defined by a cone algorithm. The cone size R of the jets will be fixed 
at the value R = 0.4. 

The vetoed cross section a vcto (s, Mfj;p™ to , R) can be written as 

a vcto (s,M 2 H ;p v T cto ,R) = a(s,M 2 H ) - Aa(s,M 2 H ;p v T eto ,R) , (1) 

where er(s, Mfj) is the inclusive cross section, and Act is the 'loss' in cross section due to the jet-veto 
procedure. The jet-vetoed cross section is evaluated by using the large-M top limit. At NLO (NNLO) 

3 5. Catani, D. de Florian and M. Grazzini 
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the calculation is performed by subtracting the LO (NLO) cross section for the production of Higgs plus 
jet(s) from the inclusive NLO (NNLO) result. 

The NLO calculation is exact: apart from using the large- M top limit, we do not perform any further 
approximations. At the NNLO, the contribution Act to Eq. ([]]) is again evaluated exactly, by using the 



numerical program of Ref. [57]. To evaluate the contribution of the inclusive cross section we use the 



recent result of Ref. [21, 22], and in particular, we rely on our approximate estimate NNLO-SVC [21]. 
In the following we present both NLO and NNLO numerical results for the vetoed cross section. The 



results are obtained by using the parton distributions of the MRST2000 set [58], with densities and 
coupling constant evaluated at each corresponding order. The MRST2000 set includes (approximate) 
NNLO parton densities. 
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Fig. 7: Vetoed cross section and K-factors: NLO results at the Tevatron Run II. 

We first present the vetoed cross section at the Tevatron Run II. In Fig. [7] we show the dependence 
of the NLO results on the Higgs mass for different values of pj? to (15, 20, 30 and 50 GeV). The vetoed 
cross sections cr veto (s, M^;p™ to , R) and the inclusive cross section a(s, M^) are given in the plot on 
the left4iand side. The inset plot gives an idea of the 'loss' in cross section once the veto is applied, by 
showing the ratio between the cross section difference Act in Eq. ((IJ) and the inclusive cross section at the 
same perturbative order. As can be observed, for large values of the cut, say p^ to = 50 GeV, less than 
10% of the inclusive cross section is vetoed. The veto effect increases by decreasing p^ to , but it is still 
smaller than 30% when p™ to = 15 GeV. On the right4iand side of Fig. [7|, we show the corresponding 
K-factors, i.e. the vetoed cross sections normalized to the LO result, which is independent of the value 
of the cut. Figure [8] shows the analogous results at NNLO. In Fig. |9] we show the LO, NLO and NNLO- 
SVC K-factor bands, computed by varing renormalization (//r) and factorization (fj,p) scales in the range 
1/2Mh < Hf, < 2Mh and normalizing to the LO contribution at fj,p = hr = Mh- The calculation 
is done with p^ to = 15 GeV. Comparing Fig. |9] with the inclusive case (see Ref. [27]), we see that the 
effect of the veto is to partially reduce the relative difference between the NLO and NNLO results; the 
increase of the corresponding K-factors can be estimated to about 25%. 

The results for the vetoed cross sections at the LHC are presented in Figs. 10 and 11 for p™ to = 20, 



30, 50 and 70 GeV. At fixed value of the cut, the impact of the jet veto, both in the 'loss' of cross section 
and in the reduction of the K-factors, is larger at the LHC than at the Tevatron Run II. This effect can 
also be appreciated by comparing Fig. [l2| and Fig. [| At the LHC, the value of p^ to = 30 GeV is already 
sufficient to reduce the difference between the NNLO and NLO results to less than 10%. 

The results presented above can be interpreted according to a simple physical picture. The dom- 



inant part of QCD corrections is due to soft and collinear radiation [21]. The characteristic scale of the 



highest transverse momentum pS? ax of the accompanying jets is p5? ax ~ (1 — z)Mjj, where the average 
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Fig. 8: Vetoed cross section and K-factors: NNLO results at the Tevatron Run II. 
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Fig. 12: The same as in Fig. ^, but at the LHC and with p™ 4 = 30 GeV. 

value (1 — z) = (1 — Mjj/s) of the distance from the partonic threshold is small. As a consequence 
the jet veto procedure is weakly effective unless the value of p™ to is substantially smaller than p™ ax . 
Decreasing p™ to , the enhancement of the inclusive cross section due to soft radiation at higher orders is 
reduced, and the jet veto procedure tends to improve the convergence of the perturbative series. At the 
LHC Higgs production is less close to threshold than at the Tevatron and, therefore, the accompanying 
jets are harder. This is the reason why, at fixed p™ to , the effect of the jet veto is stronger at the LHC than 
at the Tevatron. 

When p™ to is much smaller than the characteristic scale p™ ax ~ (1 — z)Mh, the perturbative 
expansion of the vetoed cross section contains large logarithmic contributions that can spoil the conver- 
gence of the fixed-order expansion in as. Since (1 — z)Mh is larger at the LHC than at the Tevatron, 
the value of p^ to at which these effects become visible is larger at the LHC. Whereas at the Tevatron the 
perturbative calculation for p™ to = 15 GeV seems still to be reliable, at the LHC, with the same value 



of p™ to , the perturbative result suggests that the effect of these logarithmic contributions is large [27]. 
Note added. After the completion of this work, the full NNLO QCD contribution to inclusive Higgs 



boson production has been computed Q19|]. These results influence those in the present paper through 
Eq. (|l|), since in our NNLO calculation the inclusive cross section a(s,M H ) is evaluated by using the 
approximate (soft-collinear) estimate (named NNLO-SVC) of Ref. We have considered the effect 
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of the additional hard corrections of Ref. [19] and, as expected [^T, 27], we find that they are relatively 
small. The inclusive cross section at full NNLO is smaller than its NNLO-SVC approximation by about 
5% (7%) at the LHC (Tevatron Run II). This correction can directly be applied to our results. For instance, 



the NNLO K-factors in Figs, |8|, g and |11| can be modified as K -> K - AK, where AK = 0.20-0.21 
at the Tevatron and AK = 0.11-0.13 at the LHC (the variations of AK correspond to variations of the 
Higgs mass in the range considered in the Figures). 

3. The high-energy limit of H + 2 jet production via gluon fusionF] 

At the Large Hadron Collider (LHC), the main production channels of a Higgs boson are gluon fusion and 



weak-boson fusion (WBF) [|12|,|59|]. The WBF process, qq — > qqH , occurs through the exchange of a W 
or a Z boson in the t channel, and is characterized by the production of two forward quark jets [pOfl. Even 
though it is smaller than the gluon fusion channel by about a factor of 5 for an intermediate mass Higgs 



boson, it is interesting because it is expected to provide information on Higgs boson couplings [13]. In 
this respect, H + 2 jet production via gluon-gluon fusion, which has a larger production rate before cuts, 
can be considered a background; it has the same final-state topology, and thus may hide the features of 
the WBF process. 

In Higgs production via gluon fusion, the Higgs boson is produced mostly via a top quark loop. 
The computation of H+2 jet production involves up to pentagon quark loops [[fj]. However, if the Higgs 
mass is smaller than the threshold for the creation of a top-quark pair, M H < 2M t , the coupling of the 
Higgs to the gluons via a top-quark loop can be replaced by an effective coupling [^|] : this is called the 
large-Mt limit. It simplifies the calculation, because it reduces the number of loops in a given diagram 
by one. In H + 2 jet production, the large-M 4 limit yields a good approximation to the exact calculation 
if, in addition to the condition M H < 2M t , we require that the jet transverse energies are smaller than the 
top-quark mass, p ± < M t plj]. However, the large M t approximation is quite insensitive to the value 
of the Higgs-jet and/or dijet invariant masses. The last issue is not academic, because Higgs production 
via WBF, to which we should like to compare, features typically two forward quark jets, and thus a large 
dijet invariant mass. 

In this contribution, we consider H+2 jet production when Higgs-jet and/or dijet invariant masses 
become much larger than the typical momentum transfers in the scattering. We term these conditions 
the high-energy limit. In this limit the scattering amplitude factorizes into impact factors connected 
by a gluon exchanged in the t channel. Assembling together different impact factors, the amplitudes 
for different sub-processes can be obtained. Thus the high-energy factorization constitutes a stringent 
consistency check on any amplitude for the production of a Higgs plus one or more jets. 

In the high-energy limit of H + 2 jet production, the relevant (squared) energy scales are the parton 
center-of-mass energy s, the Higgs mass M|, the dijet invariant mass s JU2 , and the jet-Higgs invariant 
masses s JlH and s J2H . At leading order they are related through momentum conservation, 

S = Sj 1 j 2 + Sj-lH "I" S j2H — -^H ■ (2) 

There are two possible high-energy limits to consider: Sj 1 j 2 » Sj lH ,Sj 2ii » and Sj 1 j 2 ,Sj 2ll » 
Sj lH , Mg. In the first case the Higgs boson is centrally located in rapidity between the two jets, and very 
far from either jet. In the second case the Higgs boson is close to one jet, say to jet j\, in rapidity, and 
both of these are very far from jet j2- hi both cases the amplitudes will factorize, and the relevant Higgs 
vertex in case 1 and the Higgs-gluon and Higgs-quark impact factors in case 2 can be obtained from the 
amplitudes for q Q — ► q Q H and q g —> q g H scattering. 

4 V. Del Duca, W.B. Kilgore, C. Oleari, C.R. Schmidt and D. Zeppenfeld 
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The high-energy limit Sj ± j 2 ^> s 



»M 



We consider the production of two partons of momenta p\ and p% and a Higgs boson of momentum p H , 
in the scattering between two partons of momenta p 2 and p^, where all momenta are taken as outgoing. 
We consider the limit in which the Higgs boson is produced centrally in rapidity, and very far from either 
jet, Sj 1 j 2 3> Sj lH , Sj 2H 3> Mg, which is equivalent to require that 



Pi > Ph » Pj 



(3) 



where we have introduced the light-cone coordinates jr 1 = po±p z , and complex transverse coordinates 
P± = P x + ip y • I n the limit (g), the amplitudes are dominated by gluon exchange in the t channel, with 
emission of the Higgs boson from the i-channel gluon. We can write the amplitude for q Q — > qQ H 
scattering in the high-energy limit as 
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(4) 



where q\ = — (pi + P2), <?2 = P3 + P4> — — ki±| 2 > * = lj 2, and the i/'s are the quark helicities. In 
Eq. (Q) we have made explicit the helicity conservation along the massless quark lines. The effective 
vertex C q ' q for the production of a quark jet, q g* — ► q, contributes a phase factor [63]: its square is 1. 
The effective vertex for Higgs production along the gluon ladder, g*g* — ► H, with and off-shell g* is 

C H (gi,p H ,g 2 ) = 2g 2 M?/v (mlAi(qi,q 2 )-2A 2 {qi,q2 



(5) 

The scalar coefficients of the triangle vertex with two off-shell gluons, Ai t2 , are defined in terms of the 
form factors Ft and Fl of Ref . [ pd] ] as 



A\ = zFr/(47r) 



A 2 = i (F T qi-q 2 + F L qfq 2 ) /(Air 



(6) 



We have checked analytically that the amplitude for qg —> qg H scattering can also be written as 
Eq. (Q), provided we perform on one of the two effective vertices C^' q the substitution (for the sake of 
illustration, we display it here for the lower vertex) 
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and use the effective vertices g* g — > g for the production of a gluon jet [ |63| ] (which contribute a phase 
factor as well). The same check on the (squared) amplitude for g g — ► g g H scattering has been per- 
formed numerically Thus, in the high-energy limit (||), the amplitudes forqQ — ► qQ H, qg — * qg H 
and g g — > g g H scattering only differ by the color strength in the jet-production vertex. Therefore, in a 
production rate it is enough to consider one of them and include the others through the effective parton 
distribution function [Q, feg(x,p-%) = G(x,[Ap) + (Cf/Ca) X]/ [Qf( x iP?F) + Qf( x iP-F)]> wnere x 
is the momentum fraction of the incoming parton, /j, F is the collinear factorization scale, and where the 
sum is over the quark flavors. 

In Fig. 13 we plot the cross section in H + 2 jet production at the LHC energy y/s = 14 TeV, 
as a function of Ay, which is defined as the smallest rapidity difference between the Higgs and the jets, 
Ay = mm(\yj 1 — y H \, \yj 2 — y H |), with the kinematical constraint > y H > y j2 . The solid line is 
the exact production rate, with the amplitudes evaluated in Ref. | |4l| ] ; the dashed line is the rate in the 
high-energy limit (||), with the amplitudes evaluated using Eqns. (0)-(0). It is apparent that the high- 
energy limit works very well over the whole Ay spectrum. However, in the evaluation of the effective 
vertex @, we used the exact value of the scalar coefficients Ai >2 . A more conservative statement is to 
say that when any kinematic quantity involved in the amplitude (Q) is evaluated in the limit (||), we expect 
the high-energy limit to represent a good approximation of the exact calculation when Ay > 2. 
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Fig. 13: Cross section in H + 2 jet production in pp collisions at the LHC energy ^fa — 14 TeV as a function of Ay, with 
Mb = 120 GeV and Mt = 175 GeV. The dijet invariant mass fulfills the constraint -JSj 1 j 2 > 600 GeV. The rapidity interval 
Ay is defined as Ay = min^j/jj — y H \, \yj 2 — yn\), with the kinematical constraint > y K > y^. The solid line is the 
exact production rate; the dashed line is the rate in the high-energy limit. 

The high-energy limit Sj ± j 2 , Sj 2li 

Next, we consider the limit in which the Higgs is produced forward in rapidity, and close to one of the 

Sj 2 h ^> Sj lH , M|. This limit implies that 



jets, say to jet j\, and both are very far from jet i.e. Sj 1J2 . - , 

pi - pt > pi ! Pi - ph < ^3 



(8) 



In this limit, the amplitudes are again dominated by gluon exchange in the t channel, and factorize into 
in effective vertex for the production of a jet and another for the production of a Higgs plus a jet. For 
example, in the limit (|8|) the amplitude for q g — > q g H scattering [41] with the incoming gluon (quark) 
of momentum p2 (pi), can be written as 



(9) 



where C 9 ' ,Hg (p 1 ^ 2 ; p" 1 , Ph) is the effective vertex for the production of a Higgs boson and a gluon jet, 
9*9 9H. It has two independent helicity configurations, which we can take to be C s;H9 (j?^;p^~,£> H ) 
and C 9 ' ,Hg (pi;pi ,p H ) []62[]. High-energy factorization also implies that the amplitude forgg^ggH 
scattering can be put in the form (Q), up to replacing the incoming quark with a gluon via the substitu- 
tion (|7]). Likewise, the amplitude forqQ—>qQH scattering can be written as 



iM^ HQ (p^;pl\pn\P^;p A 



2s[gT^ a2 C^^;pl\p H )] 



1 r 



9K- a4 Ctefa*: 



PT, 



(10) 



where C^' Hq (p^ 1 ; p" 1 ,p H ) is the effective vertex for the production of a Higgs and a quark jet, g*q — ► 
qH. There is only one independent helicity configuration, which we can take to be C^ Hq {jh 1 \P\-,Pn), 
and its expression is given in Ref. []62|], where an analysis of the limit (|h with the kinematic parameters 



of Fig. |13] can also be found. 

In conclusion, we have considered H + 2 jet production via gluon fusion, when either one of the 
Higgs-jet or the dijet invariant masses become much larger than the typical momentum transfers in the 
scattering. These limits also occur naturally in Higgs production via WBF We have shown that we can 
write the scattering amplitudes in accordance to high-energy factorization, Eqns. (Q), @ an d (^0[ ). The 
corresponding effective vertices, whose squares are the impact factors, can be found in Ref. [|62[]. 
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4. FeynHiggsl.2: Hybrid MS/on-shell Renormalization for the MSSM HiggsQ 
4.1 Introduction 



In this section we present an updated version of the Fortran code FeynHiggs [ p5[ ] that evaluates the neutral 
CP-even Higgs sector masses and mixing ang les It differs from the previous version as presented 



in Ref. [65] by a modification of the renormalization scheme concerning the treatment of subleading 
terms at the one-loop level; the two-loop corrections, for which the leading contributions of 0(ata s ) 
and 0{a 2 ) are implemented, are not affected. In particular, an MS renormalization for tan/3 and the 
field renormalization constants has been used (where the MS quantities are evaluated at the scale m t ). 
The renormalization in the new version of FeynHiggs does no longer involve the derivative of the A boson 
self-energy and the AZ mixing self-energy. This leads to a more stable behavior around thresholds, e.g. 
at Ma ~ 2m f , and avoids unphysically large contributions in certain regions of the MSSM parameter 
space. Thus, the new renormalization scheme stabilizes the prediction of the masses and mixing angles 
in the CP-even Higgs sector of the MSSM. 

4.2 Renormalization schemes 

At the tree-level, The MSSM Higgs boson masses rrih and mjj can be evaluated in terms of the SM 
gauge couplings and two MSSM parameters, conventionally chosen as Ma and tan 0. Beyond lowest 
order, the Feynman-diagrammatic (FD) approach allows to obtain in principle the most precise evaluation 
of the neutral CP-even Higgs boson sector, since in this way the effect of different mass scales of the 
supersymmetric particles and of the external momentum can consistently be included. The masses of 
the two CP-even Higgs bosons are obtained in this approach by determining the poles of the h — H- 
propagator matrix, which is equivalent to solving the equation 

[ Q 2 ~ ™h,tree + M<f) ][ <f ~ ""^tree + Z H (q 2 ) } - [ t hH (q 2 ) f = 0, (11) 

where E s , s = h, H, hH, denote the renormalized Higgs boson self-energies. For the renormalization 
within the FD approach usually the on-shell scheme is applied J6^|. This means in particular that all 
the masses in the FD result are the physical ones, i.e. they correspond to physical observables. Since 



eq. dllD is solved iteratively, the result for mh and m# contains a dependence on the field renormalization 
constants of h and H, which is formally of higher order. Accordingly, there is some freedom in choosing 
appropriate renormalization conditions for fixing the field renormalization constants (this can also be 
interpreted as affecting the renormalization of tan (3). Different renormalization conditions have been 
considered, e.g. (£' denotes the derivative with respect to the squared momentum): 



1. on-shell renormalization for T,z, Ea> ^AZ, and 5vi/v± = SV2/V2 1 67] 



2. on-shell renormalization for T*z, S^, T^az, and 5vi = 5vi di v > i = 1,2 [68] 



3. on-shell renormalization for T,z, [p7|], MS renormalization for 6Zh, SZjj, tan (3 
The previous version of FeynHiggs is based on renormalization 1 , involving the derivative of the A boson 
self-energy. The new version of FeynHiggs, see www. f eynhiggs . de, is based on renormalization 3 



(a detailed discussion can be found in Ref. [59]) 



4.3 Numerical comparison 

In this section we numerically compare the output of the previous version (based on renormalization 1) 
and the new version (based on renormalization 3) of FeynHiggs. We also show results for the recently 
obtained non-logarithmic 0(a 2 ) corrections [|8],^] that are also included in the new version of FeynHiggs. 



The comparison is performed for the parameters of the three LEP benchmark scenarios [{70J. In this way, 
the effect of the new renormalization and the non-logarithmic 0(a 2 ) corrections on the analysis of the 
LEP Higgs-boson searches can easily be read off. 

5 M. Frank, S. Heinemeyer, W. Hollik and G. Weiglein 
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In Figs. 1 6 we show the results in the "m™ ax ", "no-mixing" and "large fi" scenario as a func- 
tion of Ma (left column) and of tan (3 (right column) for two values of tan f3 (tan (3 = 3, 50) and Ma 
(M a = 100, 1000 GeV for the m^ ax and the no-mixing scenario, M A = 100, 400 GeV for the large \x 
scenario), respectively. The solid lines correspond to the new result while the dashed lines show the old 
results. The dotted lines correspond to the new result including the non-logarithmic 0(a%) contribu- 
tions. Concerning the new renormalization scheme, in the m™ ax (Fig. |14]) and the no-mixing scenario 
(Fig. |l5|) the new result is larger by «l-2 GeV for not too small M A and tan (3. For small tan/? and 
large M A the enhancement can be even larger. In the large /i scenario (Fig. |l^) the largest deviations 
appear for small tan (3 for both large and small Ma- While the previous prescription for the field renor- 
malization constants leads to unphysically large threshold effects in some regions of the parameter space, 
which arise from the AZ mixing self-energy and the derivative of the A boson self-energy, no threshold 
kinks are visible for the result based on the new renormalization. The shift in of « 1-2 GeV related 
to the modification of the renormalization prescription lies in the range of the anticipated theoretical 



uncertainty from unknown non-leading electroweak two-loop corrections WW. The new 0(af) correc 



tions can further increase by up to 
elsewhere [|72|]). 



^3 GeV for large t mixing (a detailed analysis will be presented 



The new version of FeynHiggs can be obtained from www . f eynhiggs . de . 




Ma [GeV] 



tan P 



130 



120 



110 

> 

o 

^ 100 



90 



400 



600 



800 



1000 



80 



70 



.1111 












i\ \ 

:\ \ \ 


■•' 7/ \ \ 

// \ ! 
/ // \ ! 
/ // \ ! 




- \ \ ' 
: \ \ 


■ --' // \ ! 

// 











0.5 



10 



50 



Fig. 14: The new renormalization (3, solid) and the old scheme (1, dashed) are compared in the m™ ax scenario. The dotted 
line shows the inclusion of the non-logarithmic 0(af ) corrections. The lower curves are for tan/3 = 3 (left plot) or Ma = 
100 GeV (right). The upper curves are for tan (3 = 50 (left) or Ma = 1000 GeV (right). 
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Fig. 15: The new renormalization (3, solid) and the old scheme (1, dashed) are compared in the no-mixing scenario. The 
dotted line shows the inclusion of the non-logarithmic C(af ) corrections. The lower curves are for tan f3 = 3 (left plot) or 
M A = 100 GeV (right). The upper curves are for tan/3 = 50 (left) or M A = 1000 GeV (right). 
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Fig. 16: The new renormalization (3, solid) and the old scheme (1, dashed) are compared in the large fi scenario. The dotted 
line shows the inclusion of the non-logarithmic ©(a?) corrections. The lower curves are for tan/3 = 50 (left plot) or Ma = 
100 GeV (right). The upper curves are for tan j3 = 3 (left) or M A = 400 GeV (right). 

5. Suggestions for MSSM Benchmark Scenarios for Higgs Boson Searches at Hadron Colliders^ 
5.1 Introduction and theoretical basis 

Within the MSSM the masses of the CP-even neutral Higgs bosons are calculable in terms of the other 
MSSM parameters. The lightest Higgs boson has been of particular interest, since its mass, m^, is 
bounded from above according to < Mz at the tree level. The radiative corrections at one-loop order 
have been supplemented in the last years with the leading two-loop corrections, performed by renormal- 



tial calculation [| 
(FD) approach 




ization group (RG) methods [J73fJ , by renormalization group improvement of the one-loop effective poten- 
, by two-loop effective potential calculations []|,|9|], and in the Feynman-diagrammatic 
56]. These calculations predict an upper bound for rrih of about rrih < 135 GeV.Q 

After the termination of LEP, the Higgs boson search has now shifted to the Tevatron and will 
later be continued to the LHC. Due to the large number of free parameters, a complete scan of the 
MSSM parameter space is too involved. Therefore at LEP the search has been performed in three bench- 



mark scenarios [70]. Besides the m™ ax scenario, which has been used to obtain conservative bounds 
on tan/3 [74], and the no-mixing scenario, the large-/! scenario had been designed to encourage the in- 
vestigation of flavor and decay-mode independent decay channels (instead of focusing on the h — ► bb 
channel). The investigation of these channels has lead to exclusion bounds [Q] that finally completely 
ruled out the large-/x scenario. 

The different environment at hadron colliders implies different Higgs boson production channels 
and also different relevant decay channels as compared to LEP. The main production modes at the Teva- 
tron will be V* — ► V(j> (V = W, Z,(j) = h, H, A) and also bb — > btxft, while the relevant decay modes 
will be <f) — > bb and (f> — ► t + t~ [10]. At the LHC, on the other hand, the most relevant process for a 
Higgs boson with rrih < 135 GeV will be gg — > h — > 77, supplemented by tt — ► tth — » ttbb. In order 
to investigate these different modes, we propose new benchmark scenarios for the Higgs boson searches 



at hadron colliders. Contrary to the new "SPS" benchmark scenarios proposed in Ref. [75] for general 
SUSY searches, the scenarios proposed here are designed specifically to study the MSSM Higgs sector 
without assuming any particular soft SUSY-breaking scenario and taking into account constraints only 
from the Higgs boson sector itself. 

The tree-level value for m/j within the MSSM is determined by tan j3, the CP-odd Higgs-boson 
mass Ma, and the Z-boson mass Mz- Beyond the tree-level, the main correction to stems from 



the t-t-sector, and for large values of tan (5 also from the 6-6-sector (see Ref. J70| ] for our notations.) 
Accordingly, the most important parameters for the corrections to are mt, Msusy (in this work we 



6 M. Carena, S. Heinemeyer, C.E.M. Wagner and G. Weiglein 

7 This value holds for nit = 175 GeV and Msusy = 1 TeV. If mt is raised by 5 GeV then the rrih limit is increased by 
about 5 GeV; using Msusy = 2 TeV increases the limit by about 2 GeV. 
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assume that the soft SUSY-breaking parameters for sfermions are equal: MsuSY '■= = M^ R = 
= M^ r ), X t (= A t — ///tan/3), and X b (= A b — //tan/3) (A tj b are the trilinear Higgs sfermion 
couplings, // is the Higgs mixing parameter.) rrih depends furthermore on the SU(2) gaugino mass 
parameter, M 2 (the U(l) gaugino mass parameter is given by M\ = 5/3 s^/cyy M 2 .) At the two-loop 
level also the gluino mass, trig, enters the prediction for mu- 
lt should be noted in this context that the FD result has been obtained in the on-shell (OS) renor- 
malization scheme (the corresponding Fortran code, that has been used for the studies by the LEP col- 
laborations, is FeynHiggs []65[ |6^]), whereas the RG result has been calculated using the MS scheme; 



see Ref. J76| ] for details (the corresponding Fortran code, also used by the LEP collaborations, is subh- 
pole [73, 76]). While the corresponding shift in the parameter Msjjsy turns out to be relatively small 
in general, sizable differences can occur between the numerical values of X t in the two schemes; see 
Refs. [|| [7(J . For this reason we specify below different values for X t within the two approaches. 

5.2 The benchmark scenarios 

In this section we define four benchmark scenarios suitable for the MSSM Higgs boson search at hadron 
colliders^. In these scenarios the values of the t and b sector as well as the gaugino masses will be fixed, 
while tan (5 and Ma are the parameters that are varied.]] It has been checked that the scenarios evade 
the LEP2 bounds [|[] over a wide range of the M^-tan /3-plane, where the variation should be chosen 
according to: 0.5 < tan/3 < 50, M A < 1000 GeV. 



(i) The m™ 5 ™ scenario: This scenario is kept as presented in Ref. [J70p , since it allows for conservative 
tan/3 exclusion bounds [Q] (only the sign of // is switched to a positive value.) The parameters are 
chosen such that the maximum possible Higgs-boson mass as a function of tan /3 is obtained (for fixed 
Msusy, and Ma set to its maximal value, Ma = 1 TeV). The parameters are|^| 

m t = 174.3 GeV, Msusy = 1 TeV, // = 200 GeV, M 2 = 200 GeV, m g = 0.8 M SU sy, 
X? s = 2 Msusy (FD calculation), V t MS = Msusy (RG calculation), A b = A t . (12) 

(ii) The no-mixing scenario: This benchmark scenario is the same as the m™ ax scenario, but with van- 
ishing mixing in the t sector and with a higher SUSY mass scale to avoid the LEP Higgs bounds: 

m t = 174.3 GeV, M SU sy = 2 TeV, // = 200 GeV, M 2 = 200 GeV, m- g = 0.8M SU sy, 
X t = (FD/RG calculation), A b = A t . (13) 

(Hi) The gluophobic Higgs scenario: In this scenario the main production cross section for the light 
Higgs boson at the LHC, gg — > h, is strongly suppressed (see Ref. [|84]]). The parameters are: 

mt = 174.3 GeV, M SU sy = 350 GeV, // = 300 GeV, M 2 = 300 GeV, m~ g = 500 GeV, 
X? s = -750 GeV (FD calculation), V f MS = -770 GeV (RG calculation), A b = A t . (14) 



(iv) The small a e g scenario: Besides the channel gg — > h — ► 77 at the LHC, the other channels for 
light Higgs searches at the Tevatron and at the LHC mostly rely on the decays h — > bb and h — ► r + r _ . 

8 Here, we will comment only on the phenomenology of the lightest h boson. The couplings of the H, A and H bosons 
are also subject to important radiative corrections in the large tanbeta regime, see for instance ref. [[7^, p^]. It is customary 
to define the searches for these particles in terms of their tree-level couplings. Since the corrections to these couplings are 
strongly dependent on the value of the supersymmetry breaking parameters, a precise interpretation of these search analyses 
within the MSSM will demand a knowledge of the characteristic supersymmetry breaking parameters. In the case of the h 
boson, see [^,^| for more details. 

9 Plots that show the behavior of different Higgs production and decay channels in the four scenarios can be found at 
www . f eynhiggs . de . The numerical evaluation is based on Refs. [|80|— 183[| . 

'"Better agreement with BR(6 — > sj) constraints is obtained for the other sign of Xt (called the "constrained m™ 3 ^" 
scenario) |[78||. However, this lowers the maximum values by ~ 5 GeV. 
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In comparison to the Standard Model, both hff couplings have an additional factor of sina e g/ cos (3, 
where a c ff is the mixing angle of the neutral CP-even Higgs sector, including radiative corrections (see 



e.g. Refs. [82, 83]). If a c s is small, these two decay channels can be heavily suppressed (h — > bb can 



receive also large corrections from b-g loops [|8l|,|82j]). This case is realized for large tan (5 and not too 



large Ma (in a similar way as in the large-/i scenario [70]) for the following parameters: 



m t = 174.3 GeV, M SU sy = 800 GeV, \i = 2.5M SUSY , M 2 = 500 GeV, m- g = 500 GeV 
X° s = -1100 GeV (FD calculation), xf^ = -1200 GeV (RG calculation), A b = A t . (15) 

5.3 Conclusions 

We have presented four benchmark scenarios for the MSSM Higgs boson search at hadron colliders, 
evading the exclusion bounds obtained at LEP2. These scenarios exemplify different features of the 
MSSM parameter space, such as large nih values and significant gg — > h or h — > bb, h — > t + t~ 
suppression. In analyzing the new benchmark scenarios, it will be helpful to make use of the comple- 



mentarity of different channels accessible at the Tevatron and the LHC (see e.g. Ref. [82] for details) 



6. The invisible SUSY Higgs and Dark Matter^ 
6.1 Introduction 



Current limits [ |85[ ] on both the Higgs and the neutralino in a general SUSY model are such that it is 
kinematically possible for the light Higgs to decay into the lightest neutralino. If the decay rate is sub- 
stantial the Higgs will be mainly invisible, while its usual branching ratios will be dramatically reduced 
preventing a detection in the much studied channels at the LHC and the Tevatron. Some theoretical stud- 
ies [|86|-88] have addressed the issue of how to hunt an invisibly decaying Higgs at a hadronic machine. 



For the LHC it has been suggested to use WH/ZH production , tth [|87|] or more recently the the W 



fusion process [88]. The results for the latter are quite promising since for a luminosity of 100/6 -1 a 
branching ratio into invisibles as low as 5% is enough for Higgs discovery. The aim of the present report 
is to summarize our findings on the size of the branching ratio of the Higgs into neutralinos, taking into 
account the latest data from colliders as well as from cosmology. 

6.2 MSSM parameters and h — ► XX 

For a substantial branching fraction of the Higgs into invisible to occur one needs both enough phase 
space for the decay as well as a large enough coupling of the Higgs to neutralinos. Considering that the 
present experimental and theoretical limits on the lightest MSSM Higgs implies that its mass lies in the 
interval 113 — 135 GeV, the maximum LSP mass must be below 55 — 65GeV. In models with gaugino 
unification where M\ m A'h/2 , the lower limit on the chargino mass (which depends essentially on 
M2 and n) turns into a lower limit on the neutralino mass, leaving only a small window for the Higgs 



into neutralinos. In fact in this type of models we found that the branching is never above 20% [89] 



(see also ref. [90]). For this reason we will relax the relation between M\ and Mi and consider these as 



independent parameters. In order that the coupling of the LSP to the Higgs be large, it can be shown that 



the LSP has to be a mixture of gaugino and Higgsino [|89p . However a light LSP, which corresponds to 
Mi small, is mostly a Bino. To have a non negligible Higgsino component implies that fi should be small 
as well. However /j, is bounded below by the chargino mass constraint. From these arguments, we can 
already expect that if the Higgs invisible decay is large then the chargino and next to lightest neutralino 



should not be far above the present LEP limit. One also finds [ |89| | that positive /x values lead to larger 
couplings. Large tan j3 values also lead to large Higgs mass and more phase space for the invisible 
decays, however the LSP mass increases even faster with tan j3, and we found that their coupling to 
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the Higgs get smaller with increasing tan (3 J89|]. Therefore the largest effect for the Higgs occurs for 
moderate tan (3 and we will consider tan /3 = 5. 

We take a model with a common scalar mass ttiq (defined at the GUT scale) for the SUSY breaking 
sfermion mass terms of both left and right sleptons of all three generations. As for the gaugino masses, 
we take Mi = rMi at the weak scale. For r < 1/3 or so, this scheme leads to almost no running of the 
right slepton mass, since the contribution from the running is of order M 2 , while left sleptons have an 
added Mf contribution and would be "much heavier". Indeed, neglecting Yukawa couplings one has 

m\ R = ml + .88 r 2 M 2 2 - sin 2 6 W M 2 Z cos(2/3) 

m\ L = ml + (0.72 + .22 r 2 )M 2 2 - (.5 - sin 2 W )M 2 Z cos(2/3) (16) 

Even with a common scalar mass squarks are much heavier than sleptons, since they receive a large 
contribution from the SU(3) gaugino mass. For simplicity, we then assume all squarks to be heavy 
(ITeV). In any case heavy squarks especially stops would be required in order to get a heavy enough 
light Higgs. Of course, to allow for a low ji in this scenario one needs to appropriately choose the soft 
SUSY Higgs scalar masses at high scale. It is important to stress that the kind of models we investigate in 
this report are quite plausible. The GUT-scale relation which equates all the gaugino masses at high scale 
need not be valid in a more general scheme of SUSY breaking. SUGRA models with general kinetic 
terms [94,^5 ], superstring models with moduli-dominated or with a mixture of moduli and dilaton fields, 



as well as anomaly-mediated SUSY breaking mechanisms, all lead to non-universality of the gaugino 
masses 



6.3 Constraints 

Our scenario requires as large a Higgs mass as possible without a too large value for tan (3 . We will 
then only consider the MSSM in the decoupling limit with Ma ~ ITeV and choose large enough stop 
masses (mj =lTeV) and large mixing (A t = 2.4TeV). With these parameters we have rrih = 125GeV 
for tan (3 = 5 and we are never in conflict with the lower limit on the Higgs mass mu > 113GeV. 

The limits on Mi, M2, fi, the key ingredients for this analysis, are set from the chargino mass limit 
at LEP2, m ± > 103GeV [|85l]. This bound can be slightly relaxed depending on tan (3 and the sneutrino 
mass, however we prefer to take the strongest constraint so that our results are more robust. In addition 
to this, one must include the limits from LEP2 on pair production of neutralinos, as well as the limit on 
the invisible width of the Z. For the parameters we have studied these two constraints are weaker than 
the chargino mass constraint. We will also take mr > 96GeV, for all sleptons I, even though the limit on 
the lightest stau is slightly lower [|85|]. 

Apart from the chargino mass limit, the most important constraint comes from the relic density 
of the LSP. In the models we are considering the LSP is mainly (but not totally) a bino. Since it is 
rather light the main annihilation channels are into the light fermions. The largest contributions are 
from processes involving "right-handed" sleptons since they have the largest hypercharge. In this case 
the relic density may be approximated as S7/i 2 ~ 10~ 3 m~ /m 2 . (all masses in GeV) which imposes a 
strong constraint on . However this approximation does not hold if the neutralino mass is such that 
annihilation through the Z pole, X1X1 ~~ * Z, occurs. In this case the contribution of this channel alone is 
enough to bring the relic density in the relevant range irrespective of the slepton mass. 

We use a new code [ |9l| ] for the calculation of the relic density that tackles all s-channels poles, 
threshold effects and includes all co-annihilations channels. The program extracts all exact matrix ele- 
ments from CompHEP [92] and is linked to FeynHiggs [0,^5|] for the Higgs mass. Radiative correc- 
tions to Higgs partial widths are extracted from HDECAY [|80|]. 

Fig. la shows the allowed parameter space in the M2, /i plane with tan [3 = 5 and Mi = M2/5 
for a light slepton, too = lOOGeV. The chargino mass limit from LEP2 is delimited by a line. It does not 
depend on mo. The direct LEP2 limits, expectedly, cut on the lowest /z, M2 region. This is in contrast 
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to the relic density requirement which depends sensitively on mo- We delineate three regions set by the 
relic density: a) the overclosure region Mi 2 > .3 which we consider as being definitely ruled out 



b) .1 < fill < .3 which is the preferred region and c) fth~ < .1 where one needs other form of Dark 
Matter than the SUSY Dark Matter considered here. As mo increases the allowed region for the relic 
density shrinks. However there always remain allowed regions that correspond essentially to the pole 

M2/IO (Fig. lb), the effect of the Z-pole would be seen only 
94GeV, the relic density constraint leaves a sizeable allowed 



annihilation X1X1 



Z [89]. For Mi 



at much larger values of M^. For mo 
region, however as soon as mo increases the region allowed is restricted to the regions of parameter space 



where m 



x 



M z /2. 



In view of the latest theoretical calculations of the muon anomalous magnetic moment, showing 



consistency between the experimental limit and the SM within l.Qa J99j], all constraints that were previ- 
ously thought to play an important role, (in particular the preference for light smuons) disappear in the 
range of parameters considered here. Finally, we note that b — > 57 is irrelevant since the squarks and 
gluinos are assumed heavy and that we are choosing \i > anyway. 



6.4 Results 

The branching ratio into invisible due to neutralinos will be denoted by B xx . The opening up of this 
channel will not have any effect on any of the Higgs production mechanisms. This is in contrast to other 
SUSY effects on the production and decay of the Higgs, like those due to a light stop [p7|]. Thus the 
Higgs discovery significances of the different channels at the LHC (and the Tevatron) are only affected 
by the reduction in the branching ratio into bb and 77. We define (i? 77 )as the reduction factor of the 
branching ratio of h — ► bb (h — ► 77) due to invisible compared to the same branching ratio of a standard 
model Higgs with the same Higgs mass. Since in the absence of light neutralinos the width of the Higgs 
is dominated by that into bb, one has roughly R^ ~ i? 77 ~ 1 — B xx . This is well supported by our full 
analysis and therefore we will only show the behaviour of the branching into invisible. 
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Fig. 17: a) Contours of constant Br xx from .2 (far right) to .65 (far left) for M1/M1 = 1/5, mo — lOOGeV We have also 
superimposed the various constraints. The black area is excluded by the chargino mass at LEP. The other shadings refer to the 
relic density, with the allowed region (white), the overclosure region (light grey) and the region with fl < .1 (medium grey). 
The dotted lines are constant a M lines in units o/10 -9 . b)Similarly for M1/M2 = 1/10, mo ~ 94GeU 



Fig. 17 a shows the different contours in the M2 — /i plane of B xx for M\ = M2/5 and for a 
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light slepton. We see that, even after taking all constraints, we still find large branching ratio of the 
lightest SUSY Higgs into neutralinos. The largest branchings correspond to the smallest \x values, which 
as argued before maximises the Higgsino content. It is also worth stressing that even in these general 
models, the branching ratio into invisible is never larger than 70%. For a lower ratio M\ = M2/IO, 
the invisible branching ratio can reach over 60% (Fig. [L7|b). Even in the case of heavier sleptons large 
branching ratio into invisible are possible although the allowed region of in the A/2 — /i plane corresponds 
to a narrow region around the Z pole. 

We have also searched, by making a large scan over Mi, M2, fi and mo, but for fixed tan J3 = 5, 
which minimum value of Mi one can entertain. The parameters were varied in the range 10 < Mi < 
lOOGeV, 100 < Af 2 ,/i < 500GeV, 70 < m < 300GeV as given by Eq. 1. We find that, in order 
not to have too large a relic density, Mi must be above 20GeV independently of M2 and [i, as seen in 
Fig. [TJ. However, this is not a value that gives the largest branching into invisibles, largest values are 
in the range 40 < Mi < 60 GeV where one has both a significant Higgsino-gaugino mixing and a LSP 
light enough for the Higgs to decay into it. Note that this lower bound on Mi is more or less independent 



on tan j3 Q98|]. We also show the relic density as a function of M\. Note that one hits both the Z pole and 
the Higgs pole. However for the latter configurations, B xx is negligible. 

To conclude we have found that there are regions of parameter space that give a substantial branch- 
ing fraction of the lightest SUSY Higgs into invisibles that can account for the dark matter in the universe. 
We also find that these scenarios do not always require a very light slepton since we can obtain an ac- 
ceptable amount of LSP relic density through an efficient annihilation at the Z pole. However scenarios 
with the largest branching ratio into LSP do entail that the lightest chargino and at the least the next LSP 
are light enough that they could be produced at the Tevatron. The phenomenology at the Tevatron should 
somehow be similar to the Sugra SU(5) based "24-model" which was studied in [^5|]. Among other 
things, due to the fact that one has a larger splitting between the LSP and the NLSP, as compared to the 
usual unified scenario, one expects an excess of events containing many isolated leptons originating, for 
example, from a real Z coming from the decay of the NLSP. However to make definite statements about 
observability of these states at the Tevatron requires a thorough simulation. 
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Fig. 18: Large scan over Mi, M2, p., mo for tan/3 = 5. The first panel shows the branching ratio into invisibles vs Mi. The 
second panel shows the relic density as a function of Mi. 
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7. Search for the invisible Higgs in the Wh/Zh channel at the LHCQ 

One of the various methods to search for an invisible Higgs at the LHC is the associated production 
process pp — > Wh(Zh) followed by the invisible decay of h is one of them. The signatures being 
a single lepton with large transverse momentum and missing Pt, from the higgs as well as from the 
neutrino from the W decay, and lepton pair whose mass constructs to Z and missing Pt, for the Wh 
and Zh production respectively. A parton level study, which took into account the dominant irreducible 
background caused by the WZ(ZZ) production followed by the invisible Z — > vv decay, had been 



made [87Q. This had shown that that it is possible to have S/y B ~ 5.9 for the process for rrih = 120 
GeV and 100% B.R. into the invisible channel. The issue of the reach of these channels in terms of 
and the B(h — > invisibles), is being revisited here. 

In the Wh process: The backgrounds relevant for this signal are: (1) The irreducible background 
due to the WZ production followed by Z — > vv. (2) WW production followed by leptonic decays of 
both the Ws, one lepton being lost, due to low Pt or too large a rapidity. (3) The large QCD backgrounds 
caused by the production of W with jets which are lost. The lost jets can add on to the missing Pt of the 
decay v from the W and thus possibly give substantial missing Pp. (4) Another source of background 
would be the tt production with their decay producing a W pair with two b jets. This can cause a 
background if the b jets are lost along with one of the decay leptons. (5) Z + jets production will also 
give a background if the jet(s) are missidentified as a lepton. 

Before discussing the separation of the signal from the mentionned QCD backgrounds, we report 
the result of a repetition of the of the calculation of the signal and the irreducible background, with more 
modern parton densities. A calculation of the Wh signal, using the LO formulae, with \rji\ < 2.5, P T > 
100 GeV, gives 28 fb. The irreducible WZ background for the same cuts is 40 fb. The higher order 
corrections to the signal give rise to a moderate K factor and might just compensate for lepton detection 
effeciencies of 70%. Thus the above numbers, though calculated from a LO formula can be considered 
representative. All the numbers are for = 120 GeV and 100% B.R. into the invisible channel. 

The WW background mentioned in point (2) above has also been evaluated and is 9 fb. The 
canonical values for LHC used here for the soft and forward leptons that lost, giving rise to a final state 
similar to the signal are, i) a P| < 10 GeV, ii) P^ < 5 GeV or |r/ M | > 2.5, iii) P£ < 20 GeV. Incidentally 
the background also has been evaluated only at the LO. 

The big discriminant between the QCD background and the signal is the hadronically quiet nature 
of the signal. Hence one has to tune cuts and jet vetos such that we get rid of the background at low 
cost to the signal. Towards this end we have to first look at the effect of the initial state radiations on 
the kinematical distributions in the following variables: 1) Missing Pt which no longer is just the lepton 
Pt- 2)rapidity and Pt of the jet with the largest Pt in the event. 3) the same for the jet with the 
second largest Pt- The knowledge of 2 and 3 can help determine the vetos for the two b jets that will be 
produced in tt process. Results of a priliminary calculation are very encouraging and suggest that with 
a cut on the transverse mass Mt of the lepton and the missing Pt, as well as appropriate jet vetos, it 
should be possible to reduce the QCD background substantially without any harm to the signal. The last 
background due to Z+jets has been evaluated. With a rejection factor of 10 -5 against a misidentification 
of the jet as a lepton, this background can at the most be 1.2 fb, even for a missing Pt cut of 100 GeV. 

The NLO corrections to a pair of gauge bosons production [88] shows that a veto on the jet with 
Pt > 50 GeV and r] < 3 reduces them to within 20% of the Born cross-section. This also gives 
an indication that these kinds of jet vetos will work well to reduce the background. We will have to 
optimise these cuts once we after taking into account the effect of the initial state radiation on the signal. 

In the Zh process: This signal in this case is of course much smaller. The possible backgrounds in 
this case are: (1) ZZ production followed by the invisible Z — > vv decay of the Z. (2) WZ production 
with leptonic decay of the W lepton getting lost as outlined in the consideration of the background due 

12 5. Balatenychev, G. Belanger; F. Boudjema, R.M. Godbole, V.A. Ilyin and D.P. Roy 
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to the WW production in the earlier caser. This can give rise to a final state with coming from 
the Z decay and missing Pt. (3) Production of Z + 1 jet where the jet gets lost can also cause a small 
background. 

In the Zh case one would want to use a cut on the missing Pt unlike the case of the Wh signal 
where demanding a large Pt lepton automatically guranteed a large missing Pt. The signal and ZZ 
background as well as the WZ background mentioned in (2) has been calculated at the parton level. The 
numbers are for a missing Pt cut of 100 GeV, with P T < 20GeV, and \r) \ < 2.5. The size of the signal, 
ZZ and WZ background are about 8 fb, 21 fb and 3 fb respectively. The cut on the missing Pt can 
be increased without harming the signal but can bring down the QCD backgrounds effectively. Again to 
decide on the cuts to optimise the signal, one needs to know the effect of the initial state radiation on the 
Zh signal.These calculations are in progress. 



8. Simulation of neutral Higgs Pair Production in PYTHIA using HPAIR Matrix Elements^ 
8.1 Introduction 

Scalar Higgs boson pair production at LHC allows to study the trilinear Higgs self couplings in the 
Minimal Supersymmetric Standard Model (MSSM) scheme. In pp collisions the dominating process is 
the gluon fusion gg — > HH, where H can be any of A, H or h. In PYTHIA those processes can be 
generated through the resonance channel H — ► hh or via fa fa — ► AH, Ah. But for Standard Model (SM) 
and MSSM scenarios with high values of tan (3 the contribution of s-channels becomes negligible, thus 
PYTHIA alone can not be used to explore this region of the parameters space where the cross section can 
rise above 1 pb for values of rriA, up to 150 GeV. At LHC measuring trilinear self couplings in production 
modes dominated by the gluon fusion would require a huge amount of data. Nonetheless the MSSM high 
tan (3 values cases are also interesting as a discovery channel, as their cross section is large enough. 

This note presents the implementation of LO matrix elements calculated with HPAIRF]into PYTHIA 



6.1 niOOQ - This implementation will allow a more complete simulation of hh production in resonance 
region as well as all neutral Higgs pairs in continuum production. After a brief description of the im- 
plementation steps, we will present a comparison between cross sections values obtained with PYTHIA 
Monte Carlo generation and those computed with VEGAS. For a more complete version of this work, 
refer to the ATLAS note about to be published. 



8.2 Neutral Higgs pair production in PYTHIA 

The already implemented processes into PYTHIA are fa fa — ► Ah and fa fa — > AH . Although those 
processes are dominant in most MSSM scenarios, above tan/3 = 30 they contributes for only 10% 
of the total cross section. In those cases, the dominant process is gg — ► AA and all others processes 
contributes for 50% of the total cross section, as can be seen on figure 19. Thus for a more complete 
study of those scenarios we have added into PYTHIA, the following processes [101]: gg — > h°h°, 
gg — ► H°h°, gg — > H°H°, gg — > Ah , gg — ► AH and gg — > AA. All processes described above were 
implemented in a private version of PYTHIA, as standard 2^2 PYTHIA processes. We tried to keep 
the modifications into PYTHIA to the minimum. 

For each event tried, the kinematic is chosen by PYTHIA, taking into account resonance and 
continuum production in the shape of the phase space. The PYSIGH routine then call the PYHPAIR 
function to get the differential cross section. This differential cross section is computed with HPAIR 
Matrix Elements using PYTHIA parameters. Finally, the event is selected or rejected according to the 
PYTHIA standard Monte-Carlo procedure. In SM and MSSM large tan f3 scenarios, the Higgs width 
might be non negligible and Higgs particles might be produced off-shell. Thus in order to reproduce 
correctly the kinematic of the events, this should be taken into account. 

13 R. Lafaye 

14 HPAIR is a program written by M. Spira. 
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Fig. 19: Neutral Higgs pair LO cross sections as a function of ttla in maximal mixing for tan /3 — 30. Below the transition 
limit, the cross section is dominated by AA, Ah and hh production, while above AA, AH and HH dominate. The transition 
region is special as A, H and h have all very similar masses and all possible Higgs pairs can be produced. Higgs pair 
production studies can not disentangle different Higgs, as dominating processes are made of degenerated Higgs with similar 
masses and branching ratios. 

When MSTP{42) = 1 Higgs are produced off-shell by PYTHIA. Those off-shell values are then 
used to compute weight factors for the cross section like in WW production processes. The generation 
of events for those cases is much slower and can sometimes leads to unphysical and very high differential 
cross sections. Thus, except in the purpose of studying the differences in kinematics, one should rely on 
on-shell bosons production. 

For each off-shell Higgs the cross section is weighted by a factor B* H / Bh from the propagators 
of the virtual Higgs, defined as follow: 
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where H can be any of A, H and h. The total cross section suppression factor is then (i 
the two identical Higgs in the final state) 
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8.3 Monte-Carlo results 

The stand alone HPAIR program uses VEGAS to integrate the total cross section of Higgs pair production 
processes. It can also gives NLO order results with QCD corrections. Although, those corrections are not 
valid for tua larger than 200 GeV. As a check, one can compare cross section obtained using PYTHIA 
and VEGAS with the same matrix elements. Comparisons were made with 1000 VEGAS iterations and 
100 PYTHIA events per bin. Most of the cross sections obtained are very similar. But unfortunately, 
the stand alone HPAIR code and PYTHIA have some small differences in their Higgs mass spectrum 
algorithms. Thus the cross section obtained might be slightly different, especially in cases where it is 
sensitive to the H, h mass difference like for the resonance gg — > H — > hh process at low tan (3. 

Kinematic distributions have been investigated for the 4 b final state, that is when both Higgs 
disintegrate to bb. The most important for this final state is the pr distribution of the b, as current 
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selection algorithms ask for b with a px above 40 GeV. Furthermore, as there is no 6-trigger in ATLAS, 
high pt jets will be more likely to pass the jet trigger thresholds. 

First of all the distributions in 7] and p? of the b for the existing gg — > H — > hh PYTHIA 
process 152 and our implemented process gg — > /i/i process in the resonance region (tan /3 = 3 and 
= 300 GeV) are very similar. We then take a look at the px distributions of the b and mass resolution 
of the b pairs for tan j3 = 50 and mjj = 150 GeV for on-shell and off-shell Higgs bosons in the 
maximal mixing scenario. All Higgs pair production processes were turned on, including the already 
existing // — ► AH and // — ► Ah. For this set of MSSM parameters, Higgs pair production are largely 
dominated by A and H pairs. The following table summarize the MSSM parameters and cross sections 
obtained for those processes: 



tan/? 


m A 


m H 
(GeV) 


T H 


a 
(pb) 


Number of events 
generated 


Number of events 
tried 


50 


148.9 


150.0 


9.1 GeV 


1.46 


20000 


113186 


50 


148.9 


150.0 


9.1 GeV 


1.45 


20000 


334236 



Events were then analyzed with the ATLAS fast simulation, ATLFAST, to take into account the 
detector resolution. The b pairs invariant mass resolution of 15 GeV is convoluted by the Higgs width of 
9.1 GeV, when bosons are produced off-shell, as expected. The distribution of the pt of the jets is very 
similar and makes very little difference. 



8.4 Conclusion and prospects 

Higgs pair production processes have been successfully implemented in PYTHIA 6.1. Although in some 
rare points of the parameters space when using off-shell bosons the behavior of the matrix elements 
should be investigated. Neutral Higgs bosons production has two main interests [102] [103]. The first 
is Higgs discovery and this could be achieved through the study of the resonance production where the 
cross section according to NLO predictions is of the order of 2 pb or for high values of tan /9 where the 
cross section for tan f3 = 50 can reach 30 pb. The first case as already been studied in [ 104 ] while the 
second is under analysis. An other interest is the reconstruction of the Higgs potential which implies the 
measurement of the Higgs self couplings like Xnhh- Whether this coupling can be measured for double 
Higgs resonance production with sufficient integrated luminosity at the LHC is under investigation. 



9. Multiple NMSSM Higgs boson signals at the LHC0 



In Ref. [J105J] a no-loose theorem was established for the NMSSM [jl0(J|, guaranteeing that the LHC will 
discover at least one neutral NMSSM Higgs boson (unless there are large branching ratios for decays to 
SUSY particles and/or to other Higgs bosons). Here, we try to establish the plausibility of a NMSSM 
scenario in which multiple neutral Higgs boson detection is possible at the LHC, with a number of 
available Higgs states in excess of those pertaining to the MSSM. Similarly to what done there, we only 
consider the 'direct' production channels (at the accuracy described in [ ]107| ]), namely (V = W^,Z, 
Q = b,t and q( ) refers to any possible quark flavour): 

gg — ► Higgs (gluon — gluon fusion), qq^ ^ — > V Higgs (Higgs — strahlung), (17) 
qq( ) — * qq( ^ Higgs (VV — fusion), gg, qq — > QQ Higgs (quark associated production). 



We neglect 'indirect' Higgs production via decays/bremsstrahlung off SUSY particles [ 108| ] and Higgs 
production in association with squarks pl09Q. 



The parameter scan performed here is somewhat different though. By using the program described 
of Ref. [ |1 lfjp , we have first constrained the soft terms of the NMSSM by requiring universality at the 

15 C. Hugonie and S. Moretti 
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Fig. 20: Total number of events produced through processes ( pjp at the LHC after 300 fb _1 in the NMSSM, for the CP-even 
singlet S r (left plot) and the CP-odd singlet Si (right plot), at LHC after 300 fb _1 . (For an explanation of the colour coding, 
see the text.) 

GUT scale. The independent parameters of the model are then: a universal gaugino mass M1/2, a 
universal mass for the scalars tuq, a universal trilinear coupling Aq, the Yukawa coupling A and the 



singlet self-coupling k: see eqs. (2. 1)— (2.2) of [ |105[ ]. The (well-known) value of the Z-boson mass fixes 
one of these parameters with respect to the others, so that we end up with four free parameters at the 
GUT scale. As independent inputs characterising the NMSSM, we adopt here: mo/Mi/ 2 , Aq/Mi/ 2 , A 
and k. We then integrate numerically the Renormalisation Group Equations (RGEs) between the GUT 
and the weak scale and minimised the two-loop effective potential. Furthermore, we impose the current 
experimental bounds on (s)particle masses and couplings, especially the LEP limits on the Higgs mass 
vs. its coupling to gauge bosons, see Qlll[ ]. Finally, we assume the existence of one neutral CP-even 



Higgs boson with mass 115 GeV and sufficient coupling to gauge bosons, as hinted by LEP [112]. 

The main result of this numerical analysis, as already pointed out in Ref. [ 1 XOQ , is that the addi- 
tional couplings appearing in the Superpotential are always small: X(k) < 10~ 2 . The mixing angles 
of the additional singlet states to the non-singlet sector, being proportional to these couplings, are also 
small and the singlet sector of the universal NMSSM is then quasi decoupled. (In the non-universal sce- 



nario of the previous section, the outcome was quite different: see also Ref. [105]). Hence, the neutral 
Higgs sector consists of a quasi pure (qp) CP-even Higgs singlet state, S r , a qp CP-odd singlet, Si, and 
the doublet sector is basically MSSM-like, apart from small perturbations of order ~ A 2 , so that results 
known for the Higgs sector of the MSSM are also valid in our case. 

Fixing the mass of the lightest visible (non-singlet) CP-even Higgs at 1 15 GeV puts further con- 
straints on the parameter space of the model: we find that tan/3 is always larger than 4, the CP-odd 
doublet Higgs mass Ma is larger than 160 GeV and Mstjsy is larger than 350 GeV. In this limit, the 
CP-even doublet states are the qp interaction eigenstates. The Higgs state with mass 115 GeV is a qp H u , 
and the qp is heavy (with mass larger than 300 GeV). On the other hand, the masses of the singlet 
Higgs states, S r and Si, can vary from a few GeV to 1 TeV For each of the five neutral Higgs bosons of 
the NMSSM, we have computed the total number of events obtained by summing the rates of all produc- 



tion processes in (17), assuming 300 fb 1 as integrated luminosity, at the LHC. We have plotted these 



rates versus the mass of the given Higgs states in Fig. |20] If, as tentative threshold of detectability of a 
signal, we assume 100 events, the conclusions are quite encouraging]^]. 

At the LHC, all three non-singlet Higgs states, H u , and A, might be visible at the same time, 

16 However, we emphasise that this is not intended to be a definite claim of visibility, as the evaluation of such thresholds 
would require hadron-level simulations and detector-dependent considerations which are beyond the scope of this preliminary 
study. 
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as they are MSSM-like. In the singlet sector, S r could be visible if its mass is < 600 GeV and A is 
not too small. In the NMSSM, this covers most of the parameter space. Moreover, the CP-odd singlet, 
Si, might be visible too, for an appreciable part of the parameter space (when its mass is below 200 
GeV or so). To render this manifest, we have plotted in Fig. the total number of events produced at 
the LHC with S r in the final state, Ns r (left plot) in green (light) when the corresponding Si state is 
also visible (Ns t > 100) and in red (dark) when it is not (Ns 4 < 100). Similarly, we did for Si (right 
plot), with green (light) when the corresponding S r is visible (Ns r > 100) and red (dark) when it is not 
(N Sr < 100). 



Notice that the discovery areas of multiple Higgs boson states identified in Fig. |20j are indeed 
associated to the same regions of parameter space. In fact, a first glance at the total number of CP-odd 
singlet Si produced at the LHC might indicate that nearly all the parameter space of the model is already 
covered by the CP-even singlet S r search, as all the plotted points are in green (light). This is however 
not the case, as one can check from the left-hand plot (Ng r vs. Ms r ), where a lot of points are still 
under the 100 events threshold. The fact that one sees only green (light) points on the right-hand plot is 
due to the very high density of points considered, green (light) points being plotted after red (dark) ones. 
Hence, there are red (dark) areas, uncovered by the S r searches, behind green (light), covered, ones. 

The conclusions of this preliminary study are that, although the singlet sector of the NMSSM tends 
to decouple from the rest of the neutral Higgs spectrum in the universal case, quasi pure singlet states 
could still be found at the LHC. In fact, one has to remember that a very light CP-even Higgs state is not 
excluded by LEP searches if its coupling to gauge bosons is small enough. Such a Higgs state could be 
visible at the LHC in the form of a CP-even singlet Higgs state of the NMSSM (even with rather heavy 
masses), alongside a (light) singlet CP-odd state. Besides, often this scenario occurs where the MSSM- 
like non-singlet Higgses (H u , and A) should also be visible, hence making the whole neutral Higgs 
spectrum of the NMSSM in principle accessible at the CERN machine. (Rather similar conclusions also 
apply the the next-to-MSSM (nMSSM) of Ref. JTT31 : see JITl.) 

The caveat of our analysis is that we have not performed a full Higgs decay analysis in the 
NMSSM. One may question whether the additional Higgs states would actually be visible. For example, 
they would certainly couple to singlinos - S is always the Lightest Supersymmetric Particles (LSP) in 
our context - hence decay into the latter and thus remain undetected. This should however not be the 
case. In fact, the coupling of the singlet states to ordinary matter are generally stronger in comparison 
(of order A, whereas those to two singlinos are ~ A 3 ). So that, in the end, the main decay channels of 
singlet Higgs states should be those into detectable fermions and gauge bosons. 
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B. Higgs Searches at the Tevatron 

A. Bocci, J. Hobbs and W.-M. Yao 



Abstract 

Studies of the discovery reach for the SM and supersymmetric Higgs in Run II 
have been summarized. Combining the results from all possible decay chan- 
nels, and combining the data from both experiments, with 15 fb _1 the Tevatron 
experiments can exclude a SM Higgs at the masses up to about 190 GeV at 95 
% C.L. or discover it up to 120 GeV at the 5<r level. A great deal of effort 
remains in order to raise the performance of the accelerator and bring the de- 
tectors on line and fully operational to the level demanded by the Higgs search. 



1. Introduction 

The search for the Higgs boson and the dynamics responsible for electroweak symmetry breaking is 
the central goal of high energy physics today. The Tevatron luminosity increase provided by the Main 
Injector and Recycler, along with the upgrades of both CDF and DO detectors, will provide unique 
opportunities to search for the Higgs both in the Standard Model (SM) and in supersymmetry model 
(SUSY). The Tevatron is expected to deliver an integrated luminosity of 2 fb _1 in the first two years 
(Run Ha) and 13 fb _1 in the subsequent years (Run lib) until the LHC starts. Most of CDF and DO Run 



Ha upgrades [115] were successfully installed in spring of 2001 and are now collecting data from pp 



collisions. Since the design of Run Ha upgrades is for the initial goal of 2 fb 1 and will not survive the 



course of Run lib, it is now anticipated that the Run 2b upgrades [116], and in particular the replacement 



for the Run 2a silicon vertex detector, are necessary to carrying out this exciting program. 



2. Tevatron Run II SUSY/Higgs Workshop 

A year long workshop on the Tevatron Run II Higgs physics was held at Fermilab during 1998, a joint 
venture between CDF, DO and theory group at FermiLab. The aim is to explore the discovery sensitivities 
for the Standard Model and MSSM Higgs bosons in Run II at the Tevatron. The results is ultimately 
expressed in terms of the integrated luminosity required to either exclude the Higgs with 95% confidence 
level, or discover it with 3-a or 5-sigma statistical significance at a given mass. The details can be found 
in the report of the Higgs Working Group of the Tevatron Run II SUSY/Higgs Workshop [JTo|]. 

At the time of the Workshop, neither CDF nor DO has had full Run II detector simulation package 
available. Two complementary approaches were adopted. The first approach was based on a CDF Run I 
detector simulation with the geometrical acceptance extended to correspond to the Run Ha CDF detector. 
The second approach was based purely on SHW, a simple simulation package that uses an average of 
the expected CDF and DO detector performances as a set of parametrized resolutions and acceptances to 
perform simple reconstruction of tracking, jets, vertices and trigger objects. In addition, a multivariate 
analysis using neural network [117] has been pursued and leads to a potential gain of Higgs sensitivity 
above the conventional analysis. 

Further challenges must be met in bringing the detectors online and fully operational, and in devel- 
oping the techniques and understanding, particularly in bb jet-jet mass reconstruction and b jet tagging, 
necessary to extract the small signal of the Higgs boson from the larger Standard Model background. 
Here, we will quote the results with reasonably optimistic projections for what we might achieve after a 
great deal of hard work in the coming years. 
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Fig. 1: Jet energy resolution as a function of jet Et, comparing standard corrections based on calorimeter only with energy 
determination combining information from tracking, calorimetry, and shower max. 

2.1 bb mass resolution 

The bb mass resolution assumed in making these estimates is 10% in the central part of the distribu- 
tion. This represents a significant improvement over the 14-15% resolution obtained in Run 1. One can 
improve upon the jet energy corrections by making the best possible use of all detector information, in- 
cluding tracking, shower max, calorimeter, and muon chambers. Figure|I| shows the improvement of jet 
energy resolution possible by determining jet energy from an optimum combination of all jet informa- 
tion. A great deal of effort, presently underway, is needed to understand the jet energy corrections to the 
level required to attain 10% resolution. The required integrated luminosity for Higgs discovery scales 
linearly with this resolution. 

2.2 b tagging efficiency 

The estimates of required integrated luminosity assume that the b tagging efficiency and purity are essen- 
tially the same as in Run 1 in CDF, per taggable jet, shown in Figure ||. The better geometric coverage 
of the Run 2a and 2b silicon systems, however, is taken into account and leads to a much larger taggable 
jet efficiency. Unlike the Run 1 detector, the CDF Run II detector has a silicon vertex detector covering 
the entire luminous region, and has a 3D vertexing capability. Since the required integrated luminosity 
scales inversely with the square of the tagging efficiency (assuming constant mistagging rates), there is 
a potentially great improvement for developing high-efficiency algorithms for 6-tagging. [htpb] 

2.3 Background Systematic Studies 

Most particle searches are designed to have a small background, typically only a few events. The effective 
fractional statistical precision of the background is of the same order as the background size, implying 
that (fractional) systematic errors can be rather large. Systematic uncertainties of 30% of the background 
are common in searches. However, the Run II Higgs search will have hundreds of background events, 
which come predominantly from the direct production of vector bosons plus heavy flavors, top and single 
top, di -boson, and QCD jets production. In some cases, the magnitude or the shape of the Standard Model 
backgrounds are not known at the required level of accuracy. For example, in the vvbb channel analysis 
no estimates were made of the generic bb dijet background. This process has a very large cross section 
but tiny acceptance, and is thus not modeled reliably. In the CDF Run 1 analysis, this background was 
about half of the total, and estimated from the data. Uncertainties which affect the dijet mass spectrum 
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Fig. 3: Production cross section for Standard Model Higgs at the Tevatron as a function of Higgs mass. 
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Table 1: Fractional change in backgrounds to the WH signal as a function of q 2 . The result shows the fractional change of 
each component of the background as q 2 changes from the nominal to q 2 /2 and to 2q 2 . In the average, the absolute values of 
percent changes for the two q 2 choices were used. 



Source 


Percent change 


W + bb 
tt 

Single top 

WZ 
WH signal 


6±1% 
1±1% 
3.8 ± 0.2% 
2.4 ±0.1% 
0.4 ±0.1% 



are particularly important to control because the components have different spectra. With collider data 
in hand, the understanding will be improved using dedicated studies and by tuning the event simulations 
to match collider data control samples. 

The remainder of this section is an estimate, in advance of significant data taking, of the uncer- 
tainty in background estimates arising from systematic uncertainties in simulations. One of the possible 
sources of uncertainty has been chosen for this study. The effect from q 2 dependence on the background 
inside a mass window of ±2cr around a target Higgs mass is explored. Signal and background events 
for the channel pp — ► WH — » evbb were generated using the Pythia generator and passed through the 
detector simulation used during the Run II SUSY/Higgs workshop. The basic event selection outlined 
in the SUSY/Higgs report was applied, and the acceptance within the mass window was computed. This 
process was repeated with the q 2 scale in Pythia changed to 2q 2 and q 2 /2, ranges commonly chosen 
when assessing systematic uncertainties from simulation. The results, expressed as fractional change 
in acceptance, are shown in table [I]. Systematic effects from varying q 2 scale become important when 
the overall (fractional) uncertainty approaches 5%. As expected, the largest uncertainty comes from the 
steeply falling W + bb background. This is particularly important, because W + bb is the dominant 
component of the background. 



3. Standard Model Higgs 

The dominant SM Higgs production at the Tevatron is gluon-gluon fusion via a heavy quark loop, giving 
a single Higgs produced. The Higgs can also be produced in association with a W or Z boson via its 
couplings to the vector bosons. Figure [3] shows the production cross section for various modes as a 
function of Higgs mass. In the range below 135 GeV Higgs mass, the decay of bb dominates, and for 
larger masses the decay to W pairs dominates. 

In the gluon fusion case, since the Higgs decays predominantly to bb (for Higgs masses below 
135 GeV), there is an overwhelming background from QCD production of bb pairs. The WH and ZH 
modes, however, have been extensively studied [10] and lead to several distinct signatures in which a 
Higgs signal can be observed with sufficient integrated luminosity. 



3.1 Low-mass Higgs 

For low mass (< 135GeV) Higgs, the most sensitive signatures arise from the leptonic decays of the W 
and Z, and are denoted Ivbb, vvbb, and l + l~bb. Hadronic decays of the W and Z lead to the qqbb final 
state which suffers from large backgrounds from QCD multijet production. 

In Run 1 in CDF, all four of these channels were studied, and led to limits on the Higgs cross 
section times branching ratio to bb as depicted in Figure |] As the plots shown, the Run 1 limits are 
more than an order of magnitude above the expected Standard Model cross section. Improvements to the 
detector, coupled with much higher luminosity in Run II lead to the greatly enhanced sensitivity in the 
Standard Model search. 

Maximizing the sensitivity of the search for the Higgs in these channels depends most critically 
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Fig. 4: Limits on SM Higgs cross section times branching ratio to bb from CDF runl. 

on three things as mentioned above: attaining the best possible bb mass resolution, attaining the best 
possible b jet tagging efficiency and purity, and attaining as large a data sample as possible. 

Figure |5| shows the two 6-tagged dijet mass distribution and Figure ^ shows the background- 
subtracted signal in the Ivbb case, for a 120 GeV SM Higgs, combining data from both CDF and D0 
representing 15 fb _1 integrated luminosity, which clearly illustrates that even with the best resolution 
attainable, discovering the Higgs at Tevatron remains a major challenge. 



3.2 High-mass Higgs 

For larger Higgs masses (> 135 GeV), the Higgs decays predominantly to WW^*\ Two modes have 
been shown to be sensitive in this mass range: Ivtv (from gluon fusion production of single Higgs) and 



I I 33 (from tri-vector-boson final states) [56]. The critical issues in these search modes are accurate 



estimation of the WW background in the former channel and estimation of the W j Z+jets background 
in the latter. 



3.3 Standard Model Higgs Reach in Run II 

The integrated luminosity required to discover or exclude the Standard Model Higgs, combining all 
search channels and combining the data from CDF and DO , is shown in figure ^. The lower edge of the 
bands is the nominal estimate of the Run 2 study, and the bands extend upward with a width of about 
30%, indicating the systematic uncertainty in attainable mass resolution, b tagging efficiency, and other 
parameters. 

The figure clearly shows that discovering a SM (or SM-like) Higgs at the 5-sigma level requires a 
very large data sample: even with 15 fb -1 , the mass reach is about 120 GeV at best. A 95% CL exclusion 
can, however, be attained over the entire mass range 1 15-190 GeV with the integrated luminosity foreseen 
in Run 2b. 

4. SUSY Higgs 

In the context of the minimal supersymmetric standard model (MSSM) the Higgs sector has two doublets, 
one coupling to up-type quarks and the other to down-type quarks and leptons. There are five physical 



39 



CDF/DO Combined in WH -» Irbb 



800 



600 



400 



200 



■ Data in 1 5/fb 

: Total Background 

O Wbb/cc 

V Top+single top 

..... wz 

1 WH(120) 



+ 



+ 



^ 1 .. J., 



I 



6 * 'V:± 



BIO Blg'flT'-t 



100 , . 200 
Dijet Mass (GeV) 



300 



Fig. 5: Distribution of bb mass distribution in the Ivbb channel, showing expected signal from 120 GeV SM Higgs, combining 
15 fb" 1 of data from CDF and D0 . 
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Fig. 6: Background subtracted bb mass distribution in the ivbb channel, showing expected signal from 120 GeV SM Higgs, 
combining 15 fb _1 of data from CDF and D0 . 



40 




Higgs mass (GeV/ c 2 ) 



Fig. 7: The integrated luminosity required per experiment, to either exclude a SM Higgs boson at 95% CL or discover it at the 
3a or 5a level, as a function of the Higgs mass. These results are based on the combined statistical power of both CDF and D0 
and combining all search channels. 
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Fig. 8: CDF limits on MSSM Higgs using bbbb final state. 

Higgs boson states, denoted h, A, H, and H^. The masses and couplings of the Higgses are determined 
by two parameters, usually taken to be rriA and tan (5 (the ratio of the vacuum expectation value of the 
two Higgs doublets), with corrections from the scalar top mixing parameters. 

The light scalar h can appear very Standard-Model-like or nearly so over a larger range of MSSM 
parameter space. In this scenario the results of the search for the SM Higgs produced in the WH and 
ZH modes are directly interpretable. 

More interesting is the case of large tan (3. Since the coupling of the neutral Higgses (h/A/H) to 
down-type quarks is proportional to tan f3, there is an enhancement factor of tan 2 (3 for the production 
of bbcf), 4> = h, A, H relative to the SM rate appearing in figure |5[ This leads to distinct final states with 
four b jets; if we demand that at least three of the jets be tagged, the background from QCD multijet 
processes is relatively small. In Run 1, CDF searched for this process, and from the null result excluded 
a large swath of MSSM parameter space inaccessible to LEP, as shown in figure ||. 

Based on the Run 1 analysis, and taking into account the improved 6-tagging efficiency, figure || 
shows the regions of vtia versus tan that the Tevatron can cover for different integrated luminosities. 
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Fig. 9: Regions of MSSM Higgs parameter space where 5a discovery is possible, using SM Higgs search results. 

5. Conclusion 

Studies of the discovery reach for the SM and supersymmetric Higgs in Run II have been summarized. 
Combining the results from all possible decay channels, and combining the data from both experiments, 
with 15 fb _1 the Tevatron experiments can exclude a SM Higgs at the masses up to about 190 GeV at 
95% C.L. or discover it up to 120 GeV at the 5a level. A great deal of effort remains in order to raise the 
performance of the accelerator and bring the detectors on line and fully operational to the level demanded 
by the Higgs search. 
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C. Experimental Observation of an invisible Higgs Boson at LHC 

B. Di Girolamo, L. Neukermans, K. Mazumdar, A. Nikitenko and D. Zeppenfeld 
Abstract 

We present ATLAS and CMS simulation studies on the observability of an 
invisible Higgs boson produced via weak boson fusion at the LHC. With full 
and fast detector simulations we have checked the selection efficiency of the 
basic cuts proposed to search for such an object. The Level 1 and High Level 
trigger strategies for this purely jet and missing Et final state are discussed. 



1. Introduction 

Some extensions of the Standard Model (SM) exhibit Higgs bosons which can decay into stable neutral 
weakly interacting particles, therefore giving rise to invisible final states. In supersymmetric models, the 
Higgs bosons can decay with a large branching ratio into the lightest neutralinos or gravitinos in some 
region of parameter space Q118Q , leading to an invisible final state if R parity is conserved. Invisible Higgs 
decay also happens in models with an enlarged symmetry breaking sector, e.g. in Majoron models [[TH 



[120], where the Higgs disintegrates into light weakly interacting scalars. Another possibility arises in 



models with large extra dimensions []213|]. In Randall-Sundrum type models [121], the Higgs boson can 



mix with the scalar radion field which then predominantly decays, invisibly, to graviton states [122]. The 



invisible decay of the Higgs boson is a possibility which needs to be addressed in collider searches Q123J . 



Presently, the LEP II collaborations exclude invisible Higgs masses up to 114.4 GeV [124]. The 
presence of invisible Higgs decays makes much more difficult the Higgs boson search at hadron collid- 
ers. Phenomenological studies have been done on observability of the invisible Higgs in ZH and WH 
associated production [ |87| ] and tiH production [ |125| ]. Assuming that the Higgs boson is produced with 
SM strength, while decaying with an invisible branching fraction of « 100 %, associated ZH production 



was estimated to be sensitive to Higgs masses < 150 GeV [ |87| ] at the LHC, while ttH production might 



extend the Higgs mass range to 250 GeV Q125| ]. In recent work [ |88| ] it was shown that the LHC potential 
for the search of an invisibly decaying Higgs boson can be considerably extended by studying Higgs 
production via weak boson fusion. According to these parton level studies, 10 fb^ 1 of data should allow 
to discover these particles with masses up to 480 GeV, at the 5a level, provided their invisible branching 
ratio is 1. A method for background estimation directly from the data has also been proposed. 



The search strategy for an invisible Higgs described in [88] heavily relies on the performance 



of the ATLAS/CMS calorimetry for jets and fa reconstruction as well as on a dedicated calorimeter 



trigger. In the study presented here we basically repeat the analysis done in [88], but with a more 
dedicated simulation of the detectors. We also discuss possible Level 1 and High Level triggers for the 
most efficient on-line selection of invisible Higgs events. The efficiency of the basic selections proposed 
in J88l ] has been checked. Below we will refer to the following cuts on the Et of tagging jets, a rapidity 
gap between two tagging jets, an effective mass of tagging jets (Mjj), missing transverse momentum 
(fa) and the azimuthal angle between two jets in the transverse plane (4>jj) : 

E T l ' j2 > 40 GeV, b&l < 5.0, l^i-^l > 4.4, Wi2 < 0, (1) 

fa > 100 GeV, (2) 

M j:j > 1200 GeV, (3) 

<t>ii < 1 (4) 

A mini-jet veto (no jet with Et > 20 GeV in the -q gap between two tagging jets) and lepton veto (no 
lepton with pt > p^ 1 ) have also been used. The full set of these cuts we shall refer to, hereafter, as 
WBF cuts. 
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2. Trigger on invisible Higgs 

A purely multi-jets plus fix final state in the invisible Higgs search requires a dedicated calorimeter 
trigger both at Level 1 and at High Level. Off line analysis exploits a specific feature of the two tagging 
jets accompanying Higgs production via weak boson fusion, in particular a big gap in rapidity between 
them (1). Such a requirement on the two jet topology could be applied already in on-line selections - 
at Level 1 and in High Level trigger (HLT). Together with a cut on calorimeter fa, these cuts allow to 
suppress the QCD background rate down to an acceptable level of a few Hz, as will be shown later. 

The Forward Calorimeters of the ATLAS/CMS detectors will play a crucial role in the on-line 
and off-line selections of invisible Higgs due to the presence of two forward-backward tagging jets. The 
acceptance of the CMS Hadron Forward (HF) calorimeter ( 3.0 < \rj\ < 5.0) for these jets with Et > 
30 GeV is shown in Tab. [I], before and after cut on the rapidity gap between jets. One can see that with 
the rapidity gap constraint, almost 80 % of the Higgs events will have at least one tagging jet in the 
pseudorapidity region covered by HF. 

Table 1 : Acceptance (in %) of the CMS Forward Hadron Calorimeter for tagging jets (Et > 30 GeV) in qq —> qqH. 





no jets in HF 


one jet in HF 


2 jets in HF 


no cut on rjji — 77,-2 


49 


45 


6 


1 ~ Vj2 |>4.4 


22 


65 


13 



2.1 Level 1 trigger 

ATLAS: The implementation of a specific trigger for Weak Boson Fusion processes is still under 
discussion in ATLAS. Up to now the strategy at LVL1 is to not include a trigger on jets with \rj\ > 3.2. 
However all the information on jet energy and fa is potentially available at LVL1 for all the covered 
regions (\rj\ < 4.9). 

A discussion on the implementation of a dedicated trigger for tagging jets for the studies of the 
WBF channels is in progress. The invisible Higgs channel has been used to demonstrate how such a 
trigger is fundamental for such a search and, once it is implemented, how all the WBF channels benefit of 
the trigger redundancy that will be important for precise cross section measurements for these channels. 

Using the results obtained with the off line analysis of the invisible Higgs channel, the significance 
S/y/B at an integrated luminosity L = 10 fb _1 has been evaluated for two different regions of the jet 
rapidity acceptance, finding: (a) S/y/B « 10, if \<q\ < 4.9 and (b) S/y/B « 4, if \rj\ < 3.2. 

The LVL1 hardware in ATLAS offers the capability to recognise the jet hemisphere, therefore the 
data produced for the off line analysis have been analysed with the cuts enumerated in the following to 
evaluate the expected rate at LVL1 for the background processes: (a) 2 tagging jets; (b) \r]j\ < 4.9; (c) 
p£ > 40 GeV; (d) r] h ■ r, h < 0; (e) f T > 85 GeV. 

The signal efficiency has been evaluated to be about 95 % with these cuts for the events which 
passed off-line selections (l)-(4). Table || gives the rates for all the background channels and for the 
signal for a luminosity of L = 10 33 cm _2 sec _1 and the number of events for an integrated luminosity 
f Ldt = 10 fb _1 (corresponding to the first year of running at LHC). These preliminary results look very 
promising. Further studies are going on to evaluate the errors on these numbers with a more accurate 
simulation of the forward region. 

The trigger on the tagging jets has to be extended as much as possible in 77 to preserve a good signal 
to noise ratio, as demonstrated by the following study. The trigger efficiency versus 77 has been evaluated 
by measuring the ratio between N v , the number of events at a given LVL1 acceptance \rj\ < rjLVLi, and 
A^4.9, the number of events when \rj\ < 4.9: 

e( V ) = NJN 4 . 9 
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Table 2: Expected rate and number of events for signal and backgrounds during the first year of LHC operation. 



Process 


Rate (Hz) 


Nr. of events 


H -> inv (120 GeV) 


6.1 • 10~ 4 


6.1 • 10 3 


Wjj 


4.0- 10~ 2 


4.0 • 10 5 


Zjj 


2.1 • 10~ 2 


2.1 • 10 5 


QCD (100-150 GeV) 


1.9- 10~ 2 


1.9 • 10 b 


QCD (150-200 GeV) 


3.3- 10~ 2 


3.3 • 10 b 


QCD (200-250 GeV) 


3.5 • 10~ 2 


3.5 • 10 b 


QCD (250-300 GeV) 


2.4- 10" 2 


2.4 • 10 5 


QCD (300-x GeV) 


4.5 • 10~ 2 


4.5 • 10 5 



In figure p] the behaviour of e versus 77 is shown. It is clear from the plot that, reducing the accep- 
tance region, the signal is strongly reduced, while a lower impact is obtained on the QCD background as 
well as on the Wjj and Zjj backgrounds. 
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Fig. 1 : Level 1 efficiency versus 77 coverage of Level 1 calorimeter trigger. 
CMS: CMS will have Jet triggers at Level 1 operating over the entire calorimeter acceptance, includ- 



ing Hadron Forward calorimeter [126]. The present Level 1 calorimeter trigger table includes a set of 
inclusive Jet triggers, missing E t (MET) trigger as well as combined Jet+MET trigger. We have found 
that the combined Jet+MET trigger (with Jet trigger also implemented in the HF) is the most effective 
for the invisible Higgs selection. At low luminosity (L = 2 x 10 33 cm _2 sec -1 ), with thresholds of 70 
GeV on both Jet and MET, it provides ~ 96 % efficiency for the events selected with off-line cuts (l)-(4) 
and at an acceptable background rate of 0.6 kHz [ ]127[ - |129| ]. Such a high efficiency can be understood 
from Fig. || and Fig. || where transverse energy of the highest E t jet and calorimeter fa reconstructed in 
off-line and at Level 1 is shown for the events which passed the off-line selections (l)-(4). In the off-line 
reconstruction both E t of jet and fa are corrected for the effects of calorimeter non-linearity. Jet energy 
collections are also applied at Level 1, while it is not foreseen to correct fa at Level 1. Due to this the 
Level 1 fa spectrum shown in Fig. || is shifted in comparison with off line fa. 
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Fig. 3: Calorimeter j>T reconstructed at Level 1 (dashed his- 
togram) and off line (solid histogram) for the Higgs events 
passing cuts (l)-(4). 



2.2 High Level Trigger 

CMS: At High Level Trigger the off-line requirement (1) on the pseudorapidity gap between the two 
highest E t jets can be exploited together with the cut on fa. Full granularity calorimeter information 
is available at HLT and computer farms will perform jet and fa reconstruction like in off-line analysis. 
In Fig. |] the rate of QCD multi-jet events for L = 2 x 10 33 cm~ 2 s~ 1 after cuts (1) is shown as a 
function of the cutoff on fa. This plot has been obtained with full detector [130] and reconstruction [131] 
simulations processing QCD multi-jet events withp t bins from 15-20 GeV up to 2600-3000 GeV (about 
1M events). One can see that a cut on fa above 80 GeV will reduce the rate below 1 Hz. 
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Fig. 4: QCD multi-jet background rate after cuts (1) as a function of the threshold on fa- 
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3. Results of off-line analysis with detector simulation 

We have performed off-line analysis, partially with full and partially with fast ATLAS/CMS detector 
simulation programs, using mainly the selection criteria (l)-(4) together with mini-jet veto and lepton 
veto. As backgrounds, QCD multi-jet production and QCD and electroweak (EW) W and Z production 
with more than 1 reconstructed jet have been considered. Below, we separately provide two independent 
ATLAS and CMS analyses. 



3.1 CMS analysis 

Kinematics simulation: Signal events have been generated with PYTHIA6.158 (with CTEQ5L struc- 



ture functions) JlOOfJ for a SM Higgs of mass M#=120 GeV, produced via weak boson fusion. All 
backgrounds (except QCD multi-jet production) have been simulated with LO matrix elements inter- 
faced with PYTHIA for hadronization and additional initial and final state radiation. Colour and flavour 
information at the partem level is passed to PYTHIA. EW Zjj and Wjj events were generated with COM- 



PHEP [132] (with CTEQ5L); QCD Zjj and Wjj events (generated with CTEQ4L) were generated with 



MadCUP programs [133] based on the work of Ref. [134]. Loose selection criteria have been used to 



produce events at the parton level with tree level matrix elements: 

p> T > 20GeV, \rjj\ < 5.0, \rjji - t] j2 \ > 4.2, ^ir^ < 0, (5) 
Mjj > 900GeV, (6) 

These events are further processed through PYTHIA. Initial cross sections (in pb) given by the matrix 
element calculations with these cuts are presented in Tab. |[ Cross sections include Br(Z — > vv) and 
Br(W —y Iv) for three lepton generations. 



Table 3: Cross sections (in pb) for backgrounds as given by LO matrix element calculations with cuts (5), (6). Br(Z — ► vv) 
and Br(W — > iv) is included. 



QCDW+jj 


QCDZ+jj 


EW W+jj 


EWZ+jj 


76.0 


15.7 


4.7 


0.644 



Detector simulation: Full detector simulation has been performed for the Higgs and QCD multi-jet 
events at L = 2 x 10 33 cm _2 s _1 (on average 3.4 minimum bias events of PYTHIA MSEL=1 have been 
superimposed). One of the crucial questions of this study is a proper simulation of the tails in the px 
distribution of the QCD multi-jet background. Such tails could be due to real ]/>t from heavy quarks 
decays, but also due to a number of detector effects. To make confident estimates of such effects we 



used about 1 million QCD events, fully simulated [13C] and reconstructed [131]. As mentioned already, 
events have been generated in different p t bins from 10-15 GeV (a = 8.868 x 10 12 fb) up to 2600-3000 
GeV (a = 11.25/6). However, as will be shown later, this statistics is still not enough to directly prove 
that the QCD background could be suppressed to an acceptable level after all cuts are applied. 

The other backgrounds, QCD and EW production of Wjj and Zjj have been simulated with CM- 



SJET [135] fast simulation with no minimum bias events superimposed. 

Another key point of all searches for a light Higgs produced via weak boson fusion is the use of a 
mini-jet veto, namely a veto of events with additional soft (Et > 20 GeV) jet(s) inside the rapidity gap 
between two tagging jets. The efficiency of the mini-jet veto is expected to be sensitive to detector effects 
like calibration, electronic noise and readout thresholds, interaction of soft particles in the tracker in front 
of the calorimeter, magnetic field, or pile up activity. Since we did not expect that the fast CMSJET 
simulation can properly reproduce some of these effects, we did not evaluate mini-jet veto efficiency 
from CMSJET simulation. Instead we multiply the background efficiency by P SU rv as estimated in [ ]88| ] . 
P SU rv calculated in []136|] is a probability to radiate a jet (parton) in the rapidity gap between two tagging 
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jets. In the partem level study of [ pq ] it has been assumed that such jets will be reconstructed with 
100 % efficiency. CMS full simulation study on soft jet reconstruction [ |137j ] shows that with a dedicated 
window algorithm it is possible to reconstruct 20 GeV jets at low luminosity with reasonably good purity 
and about 100 % efficiency. The question of the false jet suppression is still under investigation. 

Results on QCD background: Fig. |5] shows the fa disuibution of the QCD jet background (blue empty 
histogram) and of the Higgs signal, for M#=120 GeV, (red full histogram) after cuts (1),(3). With an 
additional cut (4), fa for the signal events is shown as the light green histogram in Fig. ||. One observes 
that the tail in the background distribution goes well beyond 100 GeV. In Fig. |5] QCD events in the tail 
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Fig. 5: j>T for Higgs of Mh=120 GeV and QCD multi-jet background after selections (1),(3). 

come from p t bins between 300 and 600 GeV. Once the cut (4) on 4>jj is applied, no background event 
with fa > 100 GeV is left. With the statistics used in the analysis, this leads to an upper limit of ~ 
1 pb on the QCD background contribution which is about of 10 times higher than the signal expected 
after the same selections (l)-(4). The ATLAS fast simulation study (see below), which uses much higher 
statistics, shows that the QCD background can be suppressed to a negligible level with cuts (l)-(3) plus 
a cut on the minimal angle in the transverse plane between fa and a jet, or with cuts (l)-(4). 

Results on Higgs signal and QCD and EW Z+jj, W+jj backgrounds: 

Estimated cross-sections (in fb) for the Higgs and backgrounds at different steps of the event 
selection are shown in Tab. ||. Numbers in parentheses are the estimates obtained in [88]. Standard 
Model production cross-sections and Br(if — ► invisible) = 1 are assumed. 

The first row of Tab. |] presents cross sections after cuts (l)-(3) and a veto on identified electrons 
(muons) for the Wjj backgrounds (including e(/x) from r decay in W — > r + v) with pfp^ >10 (5) 



GeV and |7y e ^| <2.5. The lepton veto in [88] includes a veto on r leptons with p T T >20 GeV and 
|^e,/x,T| <2.5. Here we discuss the veto on taus separately as a lepton veto or a jet veto, dependent on 
whether the tau decays leptonically or hadronically. Fig. |6] and Fig .[7] show fa distributions for the signal 
and background events after cuts (1) and (3) and e, /j, veto for the Wjj backgrounds. 

The second row of Tab. ^ presents cross-sections after mini-jet veto. As has been mentioned, the 
efficiency of the mini-jet veto for the backgrounds is taken from [jgq]. For the Wjj background with 
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Table 4: Cross sections in fb for the background and Higgs of Mh=120 GeV assuming Br(H — > invisible)=l and Standard 
Model production cross-section for the Higgs. Numbers in parentheses are results from [|88||. 



cross section, fb 


Higgs 


QCD Zjj 


QCD Wjj 


EW Zjj 


EW Wjj 


after cuts (l)-(3) and e(/i) veto for Wjj 


238(274) 


857(1254) 


1165 (1284) 


141.5(151) 


145.1(101) 


+ mini-jet veto 


180(238) 


240(351) 


237 (360) 


116(124) 


84.5 (83) 


+ < 1 


74.7(96.7) 


48.0(71.8) 


40.0 (70.2) 


12.8(14.8) 


8.7(9.9) 




50 100 150 200 250 300 
missing E t , GeV 

Fig. 6: i> T for Higgs of M H =120 GeV after cuts (1), (3) 
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Fig. 7: j>T for backgrounds after cuts (1), (3) and e(pi) veto 
for Wjj background 



W — > tv and hadronic tau-decay, the veto on the r jet is included in this second row. The efficiency of 
the r jet veto is estimated with fast CMS JET simulation, counting events with a reconstructed energy of 
the r jet greater than 20 GeV. A mini-jet veto efficiency of 0.76 for the signal is obtained from the full 
simulation. It is lower than the efficiency of 0.87 used in JB8||. This may be due to the reconstruction of 
additional soft jets from minimum bias events. Since we plan to use tracker information to suppress such 
contributions, our estimate of the mini-jet veto efficiency for the Higgs events is conservative. 

The last row of the Tab. |] presents cross-sections after all selection cuts. Fig. || and Fig. |9] show 
distributions for the Higgs and background events after cuts (l)-(3), e(/i) veto for Wjj background 
and mini-jet veto. After all selections are applied, Tab. || shows good agreement between our estimates 
and the ones obtained in [88] (except for a 40 % lower Wjj background). Our simulations therefore 
confirm the conclusion reached in [ |88| ] : the LHC potential in the search for an invisibly decaying Higgs 
boson can be considerably extended by studying the weak boson fusion channel. 



3.2 ATLAS analysis 

Kinematics and detector simulation: Higgs production, QCD multi-jet and QCD Zjj and Wjj produc- 
tion have been generated with PYTHIA6. 158 [ ]100| ]. QCD multi-jet events were generated with PYTHIA 
MSEL= 1 subprocesses and for p t in the interval 50-300 GeV, divided in bins of 50 GeV, and for p t > 
300 GeV. QCD Zjj (Wjj) backgrounds have been produced switching on the processes 15, 30 (16, 31) 
with p t > 30 GeV. These processes generate Z(W)+1 partem only and additional jets are produced due 
to initial and final state radiation (ISR, FSR) and fragmentation. Processes 123 and 124 have been used 
to generate Higgs production via weak boson fusion. 
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Fast detector simulation with the ATLFast [138] package was performed both for the signal and 
all backgrounds. 

QCD background rejection: After cut (1) on the tagging jet topology, the QCD multi-jet background is 
about of factor 10 4 larger than the other backgrounds. Cuts (l)-(3), mini-jet and lepton veto suppress this 
background to the level of QCD Wjj and Zjj backgrounds. Considerable contributions of QCD multi- 
jet events ; ~ ] — ' ail of the fir 

3 




80 100 120 140 160 180 200 

Prmiss (GeV) 

Fig. 10: Missing energy distribution after cuts (1), (3), mini-jet and lepton veto. 

distribution in QCD multi-jet events is mainly due to semi-leptonic decays within jets (> 70% of b-jets). 
Therefore the missing energy is carried by the jet whereas it is carried by the Higgs for the signal. 

An additional cut is used to ensure that the QCD multi-jet background is well under control. We 
define an isolation variable / as the minimal angle in the transverse plane between fir and the tagging 
jets: / = min{\4){j)T) — 4>{ji,2)\)- The scatter plots in Fig. 11 show the correlation between (pjj and / for 
the signal and the different backgrounds. A cut on / > 1 removes 97% of the QCD multi-jet background 
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Fig. 11: Event distribution in 4>jj versus Isolation variable for backgrounds and signal. 

at the price of a 25% loss of the signal, but it does not affect the <f>jj region which is important for the 
final counting. It may affect however, the accuracy in the prediction of the Zjj and Wjj backgrounds from 
the experimental data since it reduces the number of useful Zjj (Z — > 2e, 2/x) and Wjj (W — > e(/x) + v) 
events by about of 40-50 % as one can see in Fig. [H]. After selection of events with I > 1 the QCD 
mini -jet contribution in the region of fa > 100 GeV becomes negligible as shown in Fig. |12|. 




Tmiss ( GeV ) 

Fig. 12: Missing energy distribution after cuts (1), (3), mini-jet and lepton veto and Isolation cut. 

Higgs signal and Wjj and Zjj background estimates: Tab. || presents the estimated cross-sections in 
fb for a Higgs of M#=130 GeV and QCD Wjj and Zjj backgrounds for the different selection criteria. 
Standard Model Higgs production cross-section and Br(H — > invisible) = 1 is assumed. Lepton veto in 
the first row of Tab. § includes veto of Wjj events with lepton (p,, e, r) of p e ^^ ,r > 5, 6, 20 GeV and 
\r] e ' M ' T | <2.5. One can see that QCD backgrounds after all cuts are about of factor 1.8 larger than in [88]. 
It is mainly due to the fact that the rejection factor due to the mini-jet veto obtained with this simulation 
is about a factor 2 smaller than the one used in [Bq]. The discrepancy in the signal is understood as FSR 
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Table 5: Cross sections in fb for the background and Higgs of Mh=130 GeV assuming Br(H — > invisible)=l and Standard 
Model production cross-section for the Higgs. Numbers in parentheses are results from mm. 



cross section, fb 


Higgs (130 GeV) 


QCD Zjj 


QCD Wjj 


after cuts (l)-(3) and lepton veto 


187 (266) 


817 (1254) 


899(1284) 


+ mini-jet veto 


146 (232) 


457 (351) 


451 (360) 


+ fa < 1 


60.1 (94.3) 


132.3 (71.8) 


125.6 (70.2) 



in PYTHIA. For the QCD backgrounds it is known that PYTHIA predicts a smaller V+3jet cross section 
than matrix element calculations, when the hard process is simulated as V+lparton events. This might 
lead to an underestimate of mini-jet activity in PYTHIA. This point requires further study. 

Since EW Wjj and Zjj backgrounds have not been simulated, we use the background fractions 
of [88] and assume that EW Wjj and Zjj events contribute ~ 20% to the total background. Tab. ^ 
presents cross sections in fb for the total background and for a Higgs of M#=130 GeV, after all cuts 
including the Isolation cut. Fig. [13| shows (fijj distributions for the signal and backgrounds after cuts 
(l)-(3), mini-jet and lepton veto and Isolation cut. 



Table 6: Cross section in fb for background and Higgs of Mh= 130 GeV after all cuts including the Isolation cut. The EW 
Vjj background has been estimated to contribute as 20 % of the total background and has been added. 



cross section, fb 


H(130 GeV) 


Zjj 


Wjj 


all cuts not including cut on <pjj < 1 


130 


446 


428 


with cut on <pjj < 1 


60.1 


158. 


150.7 
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Fig. 13: 4>jj a ft er cuts (l)-(3), mini-jet and lepton veto and Isolation cut 

Discovery potential: The observation of the invisible Higgs is fully determined by the knowledge of 
background cross sections in the search area. At present, the leading order calculations for QCD Zjj and 
Wjj backgrounds lead to uncertainties of a factor of 3 to 4, depending on the renormalization scale [ ]88| ] . 
However, these backgrounds could be directly predicted at LHC using Z — > U or W — > Iv data samples 
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Fig. 14: 4>jj distribution for the QCD Zjj background (solid line) and predicted from Z — > e + e events (dashed line). Cross 
points and the right y-axis are for the residue of the two distributions, defined as NBk! L tiI Fred ■ 



as proposed in [|88|]. Fig. |14| shows the predicted cj)jj distribution for the QCD Zjj background using 
Z — ► ee events. Residuals shown in this plot are defined as {Nskg — Np rec i) / Nyf . Residuals indicate 
some pattern which is understood as the different acceptances in n for Z — ► ee and Z — ► vv events. 
Moreover, the minimal lepton pr threshold to insure its observability induces a bias in the Zjj prediction. 
Then the systematic error on the normalisation factor is ~ 2.4%. At NLO these uncertainties should be 
negligible and they are not considered in the following. 

We define the sensitivity to invisible Higgs as 1.96 standard deviations (95 % CL) from the back- 
ground error which includes statistical error and the accuracy of the background prediction from the 
(Z — ► 2/)jj and (W — > I + z^)jj data. The systematic error on the background prediction is still under 
investigation, therefore we use the predicted accuracy of 3 % as evaluated in J88|]. One should keep in 
mind however, that if the PYTHIA estimates of QCD Vjj cross sections with mini-jet veto are correct, 
then the statistics for 11 jj and Ivjj events, which are used for the background measurement, would also 
be a factor 2 higher. Hence the 3% error would go down by a factor y/2. Tab. [7] and Fig. 15 show the 
parameter 



£ 2 = Br(H — > invisible) x 



a 



(qq^qqH) 



a 



(qq->qqH)SM 



(V) 



that can be probed at 95 % CL as a function of Mh- 



Table 7: Sensitivity to the H — > invisible signal for different Higgs masses. The first line is the cross section after all cuts. 
The two last lines give values of f 2 which can be probed at 95% CL for an integrated luminosity of 10 fb^ 1 and an expected 
background of 310 fb without and with a 3% uncertainty on the total background. 



M H (GeV) 


110 


120 


130 


140 


150 


200 


250 


300 


400 


e SU rvv(4>jj < 1) (fb) 


56.2 


61.1 


60.1 


64.6 


64.2 


58.8 


51.2 


42.5 


31.0 


¥ (%) (only stat.) 


20.7 


17.8 


18.2 


16.9 


17.0 


18.6 


21.4 


25.7 


35.1 


¥ (%) stat.+ 3% 


40.4 


34.8 


35.4 


32.9 


33.1 


36.2 


41.7 


50.1 


68.5 



We have investigated the possibility of the invisible Higgs observation in the MSSM sector. In this 
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Fig. 15: 95% CL sensitivity to H — ► invisible signal for an integrated luminosity of 10 fb' 1 . 
case, Higgs (h and H) production via weak boson fusion is suppressed by the factor 

V(qq^qqH)MSSM _ , h ,H,2 

\9y ) 



\qq-*qqH)SM 



(8) 



with (g v ) 2 = sin' 2 (a — (3) and {gy) 2 = cos 2 {a — (5). It means a Standard Model like production for h 



-,H\2 



and a strongly suppressed production for H for large Ma and tan((3). Fig. 16 shows Higgs production 
cross sections via weak boson fusion for different tan(f3) values. The black line is the Standard Model 
cross section. The marked line is the ATLAS sensitivity assuming that Br(H — > invisible) = 1. The 
lines of different colours are production cross sections for MSSM h and H for different values of tan(f3). 



In the case that gaugino mass unification is ruled out [89] and Ml/M2=5 we have for medium 
tan/3 = 5 a large Br(H — > XX) which could cause a dangerous situation for Higgs discovery at the 
LHC. At the decoupling limit (M^=lTeV) and for large stop mixing A t =2A TeV and large stop mass (1 
TeV) the lightest Higgs mass is 120 GeV and it is produced with Standard Model cross section. Fig. |l7| 
shows the region where the Higgs does not escape detection. For the region which is not covered by the 
invisible Higgs search, the Higgs boson will be detected by other decay channels. 



4. Summary 

We have presented results of two independent studies of ATLAS and CMS on the experimental observ- 
ability of an invisible Higgs produced via weak boson fusion at the LHC. Background estimates in the 
two analyses differ substantially even though the performance of the two detectors is similar for this par- 
ticular study. One of the reasons is that different MC samples are used in the analysis: purely PYTHIA 
(2^2 processes) generation of W,Z plus multi-jet backgrounds vs. generation with full LO matrix 
elements implemented as an external process into PYTHIA. The most important reason however, which 
makes a difference of more than a factor of 2 between CMS and ATLAS estimates of QCD Wjj and 
Zjj backgrounds, is different rejection factors for the mini-jet (central jet) veto. In the CMS study, the 
survival probability of the mini-jet veto has been taken from analytical calculations as a multiplicative 
factor, while in the ATLAS study it comes from PYTHIA generation of additional soft central jets be- 
tween two tagging jets and the further reconstruction of these jets with the calorimeter using the fast 
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Fig. 16: Higgs production cross sections via the WBF process for different tan(/3) values. The black line is the Standard 
Model cross section, the marked line is the ATLAS sensitivity for Br(H — » invisible) = 1. 
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detector simulation. The expected performance of mini-jet veto in weak boson fusion events requires 
further study, for all Higgs decay modes, not just the invisible decay considered here. 

Even with the more pessimistic background estimates of the ATLAS simulations it has been shown 
that an invisible branching ratio of ~ 30-40 % can be probed at 95 % CL up to Higgs mass of 250 GeV 
with the first 10 fb _1 of data. Fortunately, for the invisible Higgs search in the real experiment there 
will be the possibility to estimate Wjj and Zjj backgrounds directly from the data. In addition one can 
directly measure the mini-jet veto efficiency with these events. 

A detailed study of the possible trigger strategy at both Level 1 and High Level trigger has now 
been performed. It was shown that very high trigger efficiency (~ 95 %) for invisible Higgs can be 
achieved with an acceptable background rate, by making use of topological selections of the tagging jets 
in addition to a missing Er cutoff. 
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D. Search for the Standard Model Higgs Boson using Vector Boson Fusion 
at the LHC 

G. Azuelos, C. Buttar, V. Cavasinni, D. Costanzo, T. Figy, R. Harper, K. Jakobs, M. Klute, R. Mazini, A. Nikitenko, 

E. Richter-Was, I. Vivarelli and D. Zeppenfeld 

Abstract 

The weak boson fusion process has been suggested and discussed recently as a 
discovery process for a Standard Model Higgs boson in the intermediate mass 
range rnn < 2m z at the LHC. The additional jets in the forward region of 
the detector and the requirement of low jet activity in the central region al- 
low for a significant background rejection. In the present paper the analyses 
for the H — ► WW^*) and the H — > rr decay modes have been performed 
using a more realistic simulation of the expected performance of the LHC de- 
tectors. The results obtained confirm both the large discovery potential in the 
H — * WW^ decay channel and the sensitivity to Higgs boson decays into 
r-pairs, which is important for the determination of the Higgs boson coupling 
to fermions. 

1. Introduction 

The search for the Higgs boson is one of the primary tasks of the experiments at the Large Hadron 



Collider (LHC). It has been established in many studies over the past years [ ]11 , |12| ] that a Standard 
Model Higgs boson can be identified with a high significance over the full mass range of interest, from 
the lower limit set by the LEP experiments of 114.1 GeV ^ up to about 1 TeV. 

At the LHC the production cross section for a Standard Model Higgs boson is dominated by 
gluon-gluon fusion. The fusion of vector bosons radiated from initial state quarks represents the second 
most important contribution to the production cross section. The relative contribution depends on the 
Higgs boson mass. In the intermediate mass range vector boson fusion amounts to about 20% of the 
total production cross section and becomes more important with increasing mass. However, for this 
production mode additional event characteristics can be exploited to suppress the large backgrounds. In 
these events the Higgs boson is accompanied by two jets in the forward region of the detector originating 
from the initial quarks from which the vector bosons are emitted. Another feature of the vector boson 
fusion process is the lack of color exchange between the initial state quarks, which leads to suppressed 
jet production in the central region. This is in contrast to most background processes, where color flow 
in the t-channel appears. Jet tagging in the forward region of the detector together with a veto of jet 
activity in the central region are therefore useful tools to enhance the signal to background ratio. These 



techniques have so far been applied in the search for heavy Higgs bosons [ ]12[ ] 



The observation of the Standard model Higgs boson at the LHC in the vector boson fusion channels 



in the intermediate mass range has first been discussed in Refs. [ )36| ] and [ ]38[ ] for the H — > 77 and 
H — ► WW^ decay modes and in Ref. [37] for the H — > rr decay mode. In the framework of the Les 
Houches workshop the analyses for the WW^*' and rr decay modes have been repeated using more 
realistic simulations of the performance of the LHC detectors, including forward jet tagging and jet veto 
efficiencies. In the present study the performance at low LHC luminosity, i.e. C = 10 33 cm~ 2 sec -1 , is 
addressed, and the discovery potential is evaluated for integrated luminosity values up to 30 fb~ , which 
are expected to be reached during the first years of operation. 

2. Signal and Backgrounds 

The cross sections for the vector boson fusion process have been calculated using the programme VV2H 



[139]. Although next-to-leading order calculations are available [34], leading order cross sections have 
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been used. The size of the QCD corrections amounts to about 10% and is thus small. Another reason for 
this approach is the consistency with the background estimates, for which NLO cross section calculations 
are not available for all relevant processes. The Higgs branching ratios have been calculated using the 
programme HDECAY [80]. The values for the total cross section for the vector boson fusion process 
as well as the cross sections times branching ratios for the H — ► WW(*) and H — > tt decay mode are 
given in Table [j] as a function of the Higgs boson mass. They have been computed using the CTEQ5L 
structure function parametrization []140|]. 



Table 1: Total vector boson fusion production cross sections o(qqH) and o ■ BR(H — > VKVK'*') and o ■ BR(H —> tt) as a 
function of the Higgs boson mass. 



rriH (GeV) 


120 


130 


140 


150 


160 


170 


180 


190 


o(qqH) (pb) 


4.36 


4.04 


3.72 


3.46 


3.22 


3.06 


2.82 


2.64 


o ■ BR(H — > WW { * ) ) (fb) 
a ■ BR(H -> tt) (fb) 


531 
304 


1127 

223 


1785 
135 


2370 
64.4 


2955 
11.9 


2959 
2.8 


2620 
1.6 


2054 
1.0 



The following background processes are common to all channels considered, as described in more 



detail in Ref. [38]: 



tt production: due to the appearence of two b-jets, tt events contribute already at leading order to 
the background, if the two b-jets fulfill the identification criteria of the two tag jets. 
QCD WW background: the continuum production of W-pairs, where two or more tag jet candidates 
arise from parton emission. 

Electroweak WW background: pair production of W bosons via t-channel vector boson exchange. 

Due to the similarity to the signal process the rejection of this particular background is expected 

to be much harder than for the QCD type backgrounds. 

QCD Drell-Yan Z/j* + jet production, with Z/j* — > ee, and tt. 

Electroweak tt production: tau pair production via a t-channel weak boson exchange. 



Table 2: Cross sections times leptonic branching ratios (W — > Iv, I = e, fl, and T)for the major background processes. 



process 


Pt -cutoff 


cross-section 


tt 




55.0 pb 


QCD WW + jets 




16.7 pb 


Z/'y* + jets, Z/'y* — > tt 


> 10 GeV 


1742.0 pb 


EW WW + jets 




81.6 fb 


EW tt + jets 




170.8 fb 


Z/'y* + jets, Z/"/* — > ee/n/j, 


> 10 GeV 


3485.0 pb 


ZZ 




37.8 pb 


H — > ZZ 




0.26 - 2.5 pb 



The signal processes and all background processes except the electroweak WW and tt back- 
ground have been generated using the PYTHIA 6.1 Monte Carlo event generator [10C]. The Drell-Yan 
Z/l* +jet background has been generated using matrix element calculations for qq — > Zg and qg — > Zq 
with a Pt cutoff of the outgoing quark or gluon of 10 GeV. A summary of the major background pro- 
cesses and the relevant cross sections multiplied by the branching ratio BR(W — ► lu), where I = e,p, 
and r are listed in Table ^[ In the PYTHIA simulation initial and final state radiation (ISR and FSR) and 
fragmentation have been switched on, thereby allowing for a study of the jet activity in the central detec- 
tor region due to radiation. The CTEQ5L parametrization [|14C|] of the parton distribution functions has 
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been used in the generation of all signal and background processes. To take the spin correlations in tau 



decays properly into account, tau decays have been modelled using the TAUOLA r decay library []45|]. 



The two electroweak processes, which are not included in PYTHIA, have been generated by interfac- 



ing the matrix element calculation of Ref. [ |38| ] to PYTHIA, which was then used to perform the parton 



showering, including initial and final state radiation [142]. The W + jet background which is relevant 



for the H -> WW^ lu jj decay channel has been generated using the matrix elements from the 



VECBOS Monte Carlo [141], interfaced to PYTHIA. The fast simulation packages ATLFAST [138] and 



CMS JET [135] of the ATLAS and CMS detectors have been used to perform the detector simulation. 



3. Experimental Issues in the Search for the Vector Boson Fusion Process 

Trigger aspects: all channels considered in the following have leptons (e or /j) in the final state and can 
be triggered by either the single or the di-lepton trigger. It has been assumed that full trigger efficiency 
can be reached for a single electron or muon for Pt values above 25 GeV or 20 GeV respectively. The 
corresponding threshold values for the lepton pair triggers are 15 GeV (for e) and 10 GeV (for (j). 

Lepton Identification: it has been assumed that leptons (e and /i) can be identified in the pseudorapidity 
range, \r/\ < 2.5, with an efficiency of 90%. Hadronically decaying taus can be identified over the same 
range of pseudorapidity. The tau reconstruction efficiency is correlated with the rejection against QCD 
jets and the results obtained in detailed simulation studies [11 ,|l2] have been used. 



Jet Tagging: from the signal production process it is expected that the two tag jets are reconstructed 
with a sizeable Pt in opposite hemispheres and have a large separation in pseudorapidty. In case where 
there is no further hard initial or final state radiation the transverse momentum of the tagging jets should 
be balanced by the transverse momentum of the Higgs boson. 

In the present study the two tag jets are searched over the full calorimeter coverage of the detectors 
(| 77 1 < 4.9). For all jets a calibration has been applied which corrects the jet energy on average back to the 
original parton energy. After calibration the jet with the highest Pt in the positive and negative region 



of pseudorapidity is considered as the tag jet candidate. Detailed studies have shown [143] that this 
choice of the tag jets has a high efficiency for a correct tag jet identifiction. Since the tag jets originate 
from quarks in the incoming proton it is unlikely that they are b-jets. Consequently a b-jet veto has been 
applied in the pseudorapidity range of the detectors, where b-jet tagging is available, i.e. \rj\ < 2.5. In 
this procedure a b-tagging efficiency of 60% has been assumed with a corresponding efficiency of about 



99% for a light quark or gluon jet not to be b-tagged 012Q. 

An important question is how well the tag jets can be identified at the LHC in the presence of 
pile-up. To answer this question a full GEANT simulation of the performance of the ATLAS detector 



in which also pile-up effects have been considered, has been performed J144 ]. In this study it has been 



demonstrated that tag jets can be reliably reconstructed in the ATLAS detector and that the fast simula- 
tion package of the ATLAS detector provides a sufficiently good description of the tagging efficiency. 
Differences between the fast and full simulation have been found in the transition regions between differ- 
ent calorimeters and at very forward rapidities. The ratio between the efficiency for reconstructing a jet 
with Pt above 20 GeV as determined in the full and fast simulation has been parametrized as a function 



of Pt and r\ and has been used to correct the fast simulation results accordingly [144]. 



Jet Veto Efficiencies: at the LHC, jets in the central region can also be produced by pile-up events. In 



the full simulation study Q144| ] it has been found that after applying a threshold cut on the calorimeter 
cell energies of 0.2 GeV at low and 1.0 GeV at high luminosity, that fake jets from pile-up events can 
be kept at a low level, provided that Pt thresholds of 20 GeV at low and 30 GeV at high luminosity are 
used for the jet definition. 
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4. The H -> WW& decay mode 

In this Section the analyses of the H — * WW^ channels is briefly described. The acceptance cuts 



proposed in Ref. [38] have been used as a starting point. Finally a multi-variate optimisation has been 



performed to find the best combination of values for the cuts [ ]143[ ] for Higgs boson masses in the range 
between 150 and 170 GeV. The cuts found in this optimization have also been used to get a first estimate 
of the discovery significance outside this mass range. The signal significance may still be improved, if 
the cut optimization is done as a function of mass, in particular for lower Higgs boson masses. 

4.1 Di-lepton final states: H -> WW W -> lulu 

As discussed already previously [12, 3^] a large rejection against the tt and the WW backgrounds is 



obtained by exploiting the anti-correlation of the W spins from the decay of the scalar Higgs boson []31|]. 



Background from real taus from Z + jet production with Z — » tt can be rejected if the tau momenta 



and thereby the tt invariant mass can be reconstruced in the collinear approximation [38]. Due to the 
high Pt of the Z boson in Z + jet events it can be assumed that the neutrinos in the tau decays are 
emitted in the direction of the visible charged leptons. From the lepton momenta and the p™ lss vector 
the fractions x Tl and x T2 of the r energy carried by each lepton and thereby the tt invariant mass m TT 
can be reconstructed. For decays of real t's values of x Tl 2 in the range < x Tl 2 < 1 are expected. 
The background from Z/y* Drell-Yan production in association with jets can be efficiently rejected by 
a cut on the reconstructed transverse mass mr(Hu) of the di-lepton and neutrino system, defined as 

mr(llv) = J 2P^P™ %SS ■ (1 — cosA<p), where A(p is the angle between the di-lepton vector and the 
pmiss vec |; 0r m ^ transverse plane. 

In the event selection the following cuts have been applied: 

• Two isolated leptons with Pt > 20 GeV and \r)\ < 2.5; 

• Two tag jets with P T > 40 GeV, P| > 20 GeV and Ar] tags = \r}\ ag - rf ag \ > 3.8; 

in addition it has been required that the leptons are reconstructed within the pseudorapidity gap 
spanned by the two tag jets: rj^ g n <rji 12 < rf^ g x ; 

• Lepton Angular Cuts: A<j> u < 1.05, AR a < 1.8, cos 6^ > 0.2 

Mu < 85 GeV, P T {h, 2 ) < 120 GeV, 
where A4>u is the azimuthal separation between the leptons, cos On is the cosine of the polar 
opening angle, ARu is the separation in rj — 4> space, and Mu is the invariant mass of the di-lepton 
system. 

• Real tau rejection: events are rejected, if x Tl , x T2 > 0.0 and 
M z - 25 GeV < M TT < M z + 25 GeV; 

• Invariant mass of the two tag jets: Mjj > 550 GeV; 

• Transverse momentum balance: |P|?*| < 30 GeV. 

If no hard initial or final state gluons are radiated, it is expected that the transverse momentum of 
the Higgs boson is balanced by the transverse momentum of the two tag jets, such that an upper 
cut on the modulus of the vector 

ptot _ . pi, 2 . fymiss _,_ f)j,l . pj,2 

r T — T T T ' T T 

can be used to reject background. 

• Jet veto: no jets with Pt > 20 GeV in the pseudorapidity range \rj\ < 3.2; 

• Z/-y*, Z/j* -> tt rejection: m T {Uu) > 30 GeV. 

The additional background contributions for the signal from same-flavour leptons, of which the 
the ee- and /x^-Drell-Yan backgrounds are the dominant ones, can be efficiently rejected by tightening 
the di-lepton mass cut and by introducing a cut: 

• Mu < 75 GeV and JK? iss > 30 GeV. 
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Table 3: Accepted signal (for mn = 160 GeV) and background cross sections infbfor the H — > WW — > ep, channel after 
the application of successive cuts. For the signal the contributions via the vector boson fusion and the gluon fusion channel are 
given separately. The last two lines give the final numbers if the contributions from W — > TV — > Ivv v are added for both the 
ep and the ee/ fifi final states. 





signal (fb) 






background (fb) 






VV 


gg 


tt 


WW + jets 


z h* + jets 


total 




ni H = 


160 GeV 




EW 


QCD 


EW 


QCD 




Lepton acceptance 


25.3 


107.4 


5360 


12.9 


513.7 


3.56 


12589 


18479 


+ Forward Tagging 


10.7 


2.35 


186.4 


7.79 


1.37 


1.04 


125.8 


322.4 


+ Lepton angular cuts 


6.99 


1.46 


22.0 


0.47 


0.12 


0.40 


22.7 


45.7 


+ Real r rejection 


6.69 


1.44 


21.0 


0.42 


0.12 


0.06 


3.54 


25.1 


+ Inv. mass Mjj 


5.30 


0.89 


12.5 


0.42 


0.05 


0.06 


2.54 


15.6 


+ P* ot 


4.56 


0.63 


2.71 


0.33 


0.04 


0.05 


1.77 


4.90 


+ Jet veto 


3.82 


0.45 


0.72 


0.31 


0.03 


0.04 


1.16 


2.26 


+ M T -cut 


3.71 


0.42 


0.69 


0.30 


0.03 


0.01 


0.03 


1.06 


H -> WW ( * ] -^efi + X 


















incl. t — > e,p contribution 


4.14 


0.46 


0.71 


0.33 


0.03 


0.01 


0.03 


1.11 


H -» WW ( * ] -> ee/fip + X 


















incl. t — > e, fi contribution 


3.89 


0.43 


0.64 


0.33 


0.02 


0.01 


0.15 


1.15 



The acceptance for a Higgs boson with a mass of 160 GeV and for the backgrounds after the 
application of successive cuts is summarized for the e/x final state in detail in Table |5[ In addition to 
the signal from the vector boson fusion also contributions from the gluon gluon fusion process gg — > 
H — > WW^*) where the two tag jets are produced from initial and final state radiation, have been found 
to contribute to the final signal rate. 

All numbers given in the upper part of Table || come from direct decays into electron and muon 
final states. Di-leptons can, however, also be produced via cascade decays of tau leptons, for example, 
~W — ► tv — ► Ivv v. These contributions have also been calculated and have been added to the accepted 
signal and background cross sections. An increase of about 10% for the cross sections has been found. 
Due to the softer Pt spectra of leptons from tau decays this contribution is smaller than the one expected 
from a scaling of branching ratios. The final acceptance including the contributions from r cascade 
decays, is also given for the sum of the ee and ftfi final states. Due to the additional cuts the signal 
acceptance is slightly lower than in the e/x case. 

It has to be pointed out that the numbers for the dominant it background given in Table ||| have 
been obtained from the PYTHIA parton shower simulation. An independent estimate of that background 
has been made by using tree level matrix element calculations for it + 0, 1, and 2 — jets. In order 
to avoid double counting when adding the three contributions the procedure proposed in Ref. [38], to 
define three distinct final state jet topologies, has been adopted. For it + jets, only the two b-jets are 
considered as tag jet candidates. Initial and final state radiation in these events may lead to a rejection of 
the event due to the jet veto. A distinctively different class is defined by those ti+ 1 jet events where the 
final state light quark or gluon gives rise to one tag jet and one of the two b-jets is identified as the other 
tag jet. Finally, a third class is defined where in it + 2 jet events the final state light quarks or gluons are 
identified as tag jets. 

Using this procedure the total it background in the e\i channel has been estimated to be 1.65 fb, 
which is about a factor of 2.3 higher than the PYTHIA prediction. The largest contribution has been 
found to arise from events where one tag jet originates from a b-jet and the second one from an emitted 
parton. In the following estimate of the signal significance, a conservative approach is taken and this 
number is assumed for the it background. 

After all cuts a large signal to background ratio can be reached, which leads to an impressive 
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discovery potential for a Higgs boson with a mass around 
out that even if the larger tt background estimate is taken, 
than in the gg — > WW^*' channel considered so far [|l2, 
is much less affected by systematic uncertainties on the 
gg — ► WW (*' channel no mass peak can be reconstructed, 
an excess of events above the sum of the backgrounds, for 



160 GeV in this channel. It has to be pointed 
the signal to background ratio is much better 



145]]. Therefore, the final signal significance 
Similar to the situation in the 



background 

Evidence for a signal has to be extracted from 
example, in the transverse mass spectrum. 

Following the discussion in Ref. [38] the transverse mass of the Higgs boson has been calculated 



as 



M T = J(E!£ + E!f) 2 - (pU +p™ iss ) 2 . 



where 



£*S = J(P%) 2 + ml E V T V = J(P™ iss ) 2 + m\. 



The corresponding distribution is shown in Fig. [jjfor Higgs boson signals of 140 GeV and 160 
GeV above the total background. 
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Fig. 1: Distributions of the transverse mass Mt for Higgs boson signals of 140 GeV (left) and 160 GeV (right) above the total 
background after all cuts are applied. The accepted cross sections a a cc (infb/5 GeV) including all efficiency and acceptance 
factors are shown in both cases. The number of events observed in the detector is obtained by multiplying with the integrated 
luminosity. 



It should be noted that in the present study numbers for signal and background have been found 
which are somewhat different from the numbers quoted in the original parton level study of Ref. [38]. A 
detailed comparison between both simulations has been performed and the main differences have been 
understood. One reason for a reduced signal efficiency observed in the present study is a lower lepton 
acceptance. In addition, the efficiency for reconstructing the tag jets is found to be lower. Both are 
related to effects from initial and final state gluon radiation. They lead to a degraded lepton isolation 
as well as to non-Gaussian tails in the jet response which can not be fully corrected in jet calibration 
procedures. However, the main conclusions of Ref. [38], that the search for vector boson fusion in the 
intermediate mass range at the LHC has a large discovery potential for a Standard Model Higgs boson in 
the H — > WW^*> decay channel are confirmed. 
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4.2 The Zz^-jet-jet channel 

It has also been investigated whether the larger branching ratio of the W-bosons into quark pairs can 
be used and the process qq — * qqH — * qqWW^ — > qq Iv jj can be identified above the larger 
backgrounds, in particular the W + jet background. This process has already been established as a 
discovery channel for a heavy Higgs boson [12] in the vector boson fusion process, but has so far not 
been considered in the intermediate mass region. 

The final signal rate in this channel is expected to be much lower than the corresponding numbers 
in the di-lepton channel [146]. However, a possible observation of a Higgs boson with a mass around 
160 GeV can be confirmed in this channel for higher integrated luminosities around 30 fb -1 . It must be 
stressed that very hard cuts on the Pt and on the invariant mass of the forward tag jets, as well as on 
the separation AR between the lepton and the jets from the W-decay are necessary to extract the signal 
above the large backgrounds. These extreme cuts might also lead to larger systematic uncertainties on 
the background prediction. 



4.3 Discovery potential as a function of mass 

The analyses outlined above have been performed in the full range of Higgs boson masses from 1 10 
to 190 GeV. The expected numbers of signal and background events expected in the transverse mass 
interval 50 < My < mjj + 40 GeV are given in Table ^ for integrated luminosities of 5 and 30 fb , 
respectively, for the three l^lfW channels considered. The interval of transverse mass has been chosen 
to maximize the signal to background ratio. For signal events about 98% of all events are contained in 
that interval. For the estimate of the signal significance the more conservative matrix element estimate 
of the tt background has been used. 



Table 4: Expected signal and background rates for the three WW^*' decay channels as a function of mn assuming an inte- 
grated luminosity of 5 fb -1 and 30 fb -1 (for the Ivjj channel). In addition, the signal significances are given for an integrated 
luminosity of 5 and 30 fb -1 . They have been computed using Poisson statistics and assuming a systematic uncertainty of 5% 
on the background. 



run 


(GeV) 


110 


120 


130 


140 


150 


160 


170 


180 


190 


H -> WW { -* > ^ep + X 




















Signal 


(Sfb- 1 ) 


0.4 


1.8 


4.8 


8.4 


13.9 


22.7 


21.6 


16.9 


12.7 


Background 


(5 ft)- 1 ) 


5.2 


6.0 


7.0 


8.0 


8.3 


9.1 


9.3 


9.4 


9.8 


Stat, significance 


(5 fb" 1 ) 






1.5 


2.5 


3.9 


5.8 


5.5 


4.4 


3.4 


Stat, significance 


(30 fb -1 ) 




1.5 


3.9 


6.9 


11.1 


17.3 


16.3 


12.6 


9.3 


H -> WW { *> -> 


ee/pp + X 




















Signal 


(5 ft- 1 ) 


0.3 


1.6 


4.4 


7.9 


13.2 


21.5 


20.4 


16.6 


11.6 


Background 


(5 fb" 1 ) 


5.2 


5.9 


6.9 


7.6 


8.0 


8.6 


8.8 


8.9 


9.2 


Stat, significance 


(5 fb" 1 ) 






1.4 


2.4 


3.8 


5.6 


5.3 


4.4 


3.2 


Stat, significance 


(30 fb" 1 ) 




1.6 


3.9 


6.7 


10.8 


16.9 


15.8 


12.8 


8.8 


H -> WW W - 


Iv jj + X 




















Signal 


(30 fb" 1 ) 






4.5 


7.5 


10.5 


24.0 


24.0 


18.0 


15.0 


Background 


(30 fb" 1 ) 






6.0 


6.0 


6.0 


18.0 


18.0 


18.0 


18.0 


Stat, significance 


(30 fb -1 ) 






1.5 


2.4 


3.3 


4.6 


4.6 


3.5 


3.0 



The signal significances, expressed in the equivalent number of Gaussian standard deviations, 
have been calculated using Poisson statistics and assuming an integrated luminosities of 5 and 30 fb -1 
and a systematic uncertainty of 5% on the background. A 5a discovery can be claimed for 5 fb _1 for 
run = 150 — 185 GeV if the e/x and ee/fi/j, channels are combined. For an integrated luminosity of 30 
fb -1 the discovery range increase to ran = 130 — 190 GeV. 
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Table 5: Accepted signal (for tuh = 120 GeV) and background cross sections infbfor the H — > tt — > e/i channel after the 
application of all cuts for both the ep and the sum of the ee and pp channels. For the signal the contributions via the vector 
boson fusion and the gluon fusion channel are given separately. 





signal (fb) 
VV gg 
m H =120 GeV 


tt 


WW 
EW 


backgrc 
+ jets 
QCD 


und (fb 

TT H 

EW 


) 

- jets 
QCD 


Total 


H — > tt — > ep 


0.23 0.01 


0.02 


0.01 


0.0 


0.02 


0.04 


0.09 


H — > tt — » ee/pp 


0.24 0.02 


0.05 


0.01 


0.0 


0.04 


0.08 


0.17 



5. The H — > tt decay mode 

In the following searches for H — > tt decays using the double leptonic decay mode, qqH — > qq tt — > 
qq l + vv l~vv and the lepton-hadron decay mode qqH —* qq tt — > Z^z/ had v, are described. Due 
to the tt final state the Z + jet, Z — ► tt background contributes to the irreducible background and 
constitues the principal background for H — > tt decays at low Higgs boson masses. The main points 
of the analyses are briefly summarized in the following subsections. For details the reader is referred to 
Refs. pi and M]. 



5.1 Di-lepton final states: H -» tt -> Z+Z~ P™ ss + X 

The di-lepton final state is characterized by two tag jets in the forward regions of the detector, two leptons 
in the central region and missing transverse momentum. The following cuts have been applied to select 
e/i final states: 

• Two isolated leptons with Pr(e) > 15 GeV and \rj e \ < 2.5 and 
Pr(ji) > 10 GeV and |^| < 2.5; 

• Two tag jets with P^ > 50 GeV, P| > 20 GeV and Ar] tags = \r]} ag - rf ag \ > 4.4. 

In addition, it has been required that the leptons are reconstructed within the pseudorapidity gap 
spanned by the two tagging jets: rftag n <% 2 < Vt^'i 

• Missing transverse momentum: p™ ss > 50 GeV; 

• Invariant mass of the two tag jets: Mjj > 700 GeV; 

• Jet veto: no jets with Pt > 20 GeV in the pseudorapidity range defined by the two tag jets 

min ^ veto ^ max. 
'hag ^ 'Ij 'Itag > 

• Azimuthal separation Acj)jj between the tag jets: Afyj < 2.2. 

This cut is applied to reduce the electroweak Zjj background, for which back-to-back jets are 
preferred [|88|]. 

• Separation AR e ^ in r] — <p space between the two leptons: AR efM < 2.6; 

• Real tau reconstruction: x Tl , x T2 > and x 2 Tx + x^ 2 < 1; 

• Mass window around the Higgs boson mass: mg — 10 GeV < m TT < ma + 15 GeV. 

For ee and final states the additional background from Z decaying into ee or is efficiently 
rejected by requiring in addition: m« < rriz — 15 GeV. 

The results are summarized for both the e/i and the sum of the ee and /i/i channel in Table ^ 
where the accepted cross sections for the signal with ttih = 120 GeV and the background contributions 
are given after the application of all cuts. After t reconstruction the signal to background ratio is still 
much smaller than 1. This situation is drastically changed after the application of the mass cut around 
the Higgs boson mass. Due to the reconstructed Higgs boson mass the sidebands can be used for the 
determination of the absolute level of the background. 

The distribution of the reconstructed tt invariant mass is shown in Fig. |2| (left) for the sum of 
the e/i,ee and /i/i channels for a Higgs boson signal of 120 GeV above the background assuming an 
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integrated luminosity of 30 fb 
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Fig. 2: fe/f: 77ze reconstructed tt invariant mass for a Higgs boson signal of 120 GeV in the ll-channel above all backgrounds 
after application of all cuts except the mass window cut. right: The reconstructed tt invariant mass for a Higgs boson signal 
of 135 GeV in the (J, — had)-channel compared to the QCD plus electroweak Z + jj {Z —> tt —> I + had backgrounds 
after application of all cuts except the mass window cut. In both cases the signal and background numbers are shown for an 
integrated luminosity of 30 fb" 1 . 



The analysis has been performed for Higgs boson masses in the range from 110 to 150 GeV. The 
expected numbers of signal and background events and the statistical significance for a Higgs boson 
discovery expressed in terms of Gaussian standard deviations are given in Table ^ for an integrated 
luminosity of 30 fb -1 . 



Table 6: Expected signal and background rates and statistical significance for the three tt decay channels as a function of ma 
assuming an integrated luminosity of 30 fb -1 . 



m H (GeV) 


110 


120 


130 


140 


150 


H — ► tt — > ep, P^' ss 












Signal (30 fb" 1 ) 


7.7 


7.0 


5.1 


3.3 


1.5 


Background (30 fb" 1 ) 


7.0 


2.6 


2.3 


1.9 


1.5 


Stat, significance (30 fb" 1 ) 


2.4 


3.2 


2.5 


1.8 




H — > tt — > ee/p,p,P T alas 












Signal (30 fb" 1 ) 


9.2 


7.2 


5.7 


3.1 


1.5 


Background (30 fb" 1 ) 


10.5 


5.2 


3.8 


3.1 


2.3 


Stat, significance (30 fb -1 ) 


2.4 


2.6 


2.3 


1.4 





5.2 The lepton-hadron decay mode: H — > tt — > l^vv had v 

The I — had decay mode of the tt final state has also been studied using a fast simulation of the CMS 
detector [135]. The backgrounds considered in this study are QCD and electroweak production of Z + 
jj (Z — > tt — > I + r-jet) and W + 3j (W — > e{fj) + v). As in the previous studies, the bb + jj 
background it expected to be small [37]. QCD Z + jj production has been generated using the lowest 
order matrix element provided by the authors of Ref. [ |37| ] interfaced to PYTHIA. The electroweak Z+jj 
production has been simulated with COMPHEP [132], again interfaced to PYTHIA. The W + 3j events 
have been produced with PYTHIA which may lead to an underestimate of this background cross-section. 
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For the present study the jet veto efficiency has not yet been evaluated in a full detector simulation 



and the survival probability as determined in Ref. [ ]38| ] has been used to account for the acceptance of the 
jet veto cut. In Ref. [ [38| ] the jet veto efficiency was found to be 0.87 for signal events and background 
from electroweak production and 0.28 for QCD type backgrounds. In addition, it has been assumed that 
the jets are reconstructed with full efficiency. Results based on a full simulation of the CMS detector 
have shown that a reconstruction of low Pr jets around 20 GeV is possible with a high efficiency in a 
low luminosity scenario. For the identification of the hadronic tau with the calorimeter and the tracker 
information an efficiency of 0.32 has been used. Using these criteria a probability of 0.0019 is obtained 
for mis-identifying jets as hadronic r's. 

In the event selection the following cuts are applied: 

• One isolated lepton with P T > 15 GeV and \rj\ < 2.4; 

• One hadronic tau jet with Pt > 30 GeV and \t]\ < 2.4; 

• Two tag jets with Ar] tags = \r)} ag - rf tag \ > 4.4. 

In addition, it has been required that the lepton and the tau-jet are reconstructed within the pseu- 
dorapidity gap rfi™ 1 + 0.7 < rn,T-jet < rCg X - 0.7; 

> 1000 GeV; 



Invariant mass of the tag jets: Mjj 
Transverse mass m t (l, P™ ss ) := 



Tau reconstruction: 
Mass window: 



< x T , < 0.75, < x 



Th 



(1 - cosA<f>) 
< l; 



< 30 GeV; 



< 15 GeV. 



The number of signal events for m#=135 GeV and the number of different background events 
expected after all selections for an integrated luminosity 30 fb~ x are given in Table. [7| The errors quoted 
result from the statistical uncertainty of the Monte Carlo data samples. For comparison the number of 



events estimated in Ref. [ 37 ] are also shown in the second row of the table. 



Table 7: Number of signal (ma 
integrated luminosity 30 fb~ 



135 GeV) and background events in the I — had channel expected after all selections for an 



Higgs, Afe=135 GeV 


QCDZ+jj 


EWZ+jj 


W+3j 


6.7±0.3 


0.63±0.10 


0.74±0.08 


0.14±0.05 


6.2 (from Ref. {37]) 


total background from Ref. [37] is 1.1 



The number of signal and the total number of background events expected after all selections 
for different Higgs boson masses in the range between 115 and 145 GeV and assuming an integrated 
luminosity 30 fb^ 1 are shown in Tab. ||. 

Table 8: Expected signal and background rates and statistical significance for the l-had tt decay channels as a function of 
rriH assuming an integrated luminosity of 30 fb _1 . 



m H (GeV) 


115 


125 


135 


145 


H — > tt — > lhad P^" B 










Signal (30 fb _1 ) 


12.6 


9.9 


6.7 


3.6 


Background (30 fb -1 ) 


5.5 


2.3 


1.5 


1.1 


Stat, significance (30 fb _1 ) 


4.1 


4.5 


3.7 


2.4 



The reconstructed rr invariant mass for the QCD and EW Z + jj (Z — > tt — > I + T-jet) 
backgrounds and for a Higgs boson with m#=135 GeV is shown in Fig. |2| (right). The distributions are 
normalised to the expected number of events after all cuts, except the mass window cut, for an integrated 
luminosity of 30 fb~ l . 
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6. Conclusions 

The discovery potential for the Standard Model Higgs boson in the intermediate mass range has been 
studied using the vector boson fusion process. It has been demonstrated that the LHC experiments have 
a large discovery potential in the H — ► WW^ — ► l + l~P£ ass channel. The additional signature of tag 
jets in the forward and of a low jet activity in the central region of the detector allow for a significant 
background rejection, such that a better signal to background ratio than in the inclusive H — > WW^*\ 
which is dominated by gluon gluon fusion process, is obtained. As in the inclusive case, only the trans- 
verse mass of the Higgs boson can be reconstructed and a signal has to be claimed from an excess of 
events above the background. Due to the larger signal to background ratio in the search for the fusion 
process the signal sensitivity is less affected by systematic uncertainties on the background prediction. 
The present study shows that the ATLAS and CMS experiments at the LHC would be sensitive to a Stan- 
dard Model Higgs boson in this decay channel in the mass range between 150 and 185 GeV with data 
corresponding to an integrated luminosity of 5 fb _1 only. 

In addition, it has been shown that in the low mass region for mu < 140 GeV the LHC experi- 
ments are also sensitive to the rr decay mode of the Standard Model Higgs boson, if the characteristics of 
the vector boson fusion are exploited. However, a discovery in this final state would require an integrated 
luminosity of about 30 fb _1 and a combination of the I — I and I — had decay modes. The measurement of 
the r decay mode is particularly important for a measurement of the Higgs boson coupling to fermions. 

The present study confirms the results published earlier [|3^-38], that the search for vector boson 



fusion in the intermediate mass range at the LHC has a large discovery potential over the full range from 
the lower limit set by the LEP experiments up to 2 mz, where the high sensitivity H — > ZZ — ► 4 I 
channel takes over. 
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E. Study of the MSSM channel A/H — > rr at the LHC 

D. Cavalli, R. Kinnunen, G. Negri, A. Nikitenko and J. Thomas 
Abstract 

Sudies both from ATLAS and CMS with fast and full detector simulation have 
shown that the discovery potential of the A/H —> tt channel in the MSSM is 
large in the range from ~100 GeV to ~1 TeV already with 30 fb _1 col- 
lected at low LHC luminosity (10 33 cm _2 s _1 ). The results of these studies, in 
particular for the lepton-hadron and the hadron-hadron final decay channels, 
are presented here. The question of the trigger for the hadron-hadron final 
state that is a very important issue for this purely hadronic final state process 
is also discussed here. 



1. Introduction 

In the minimal supersymmetric extension of the Standard Model (MSSM), two Higgs doublets are re- 
quired, resulting in 5 physical states, referred to as H + , H~, h (neutral lighter scalar), H (neutral heavier 
scalar) and A (neutral pseudoscalar). At tree level their masses can be computed in terms of only two 
parameters, typically m\ and tan/3 (the ratio of the vacuum expectation values of the two doublets). 

The MSSM H — » tt and A — > tt rates are strongly enhanced with respect to the SM case over 
a large region of the parameter space. 

A/H can be produced via two different mechanisms. For low values of tan(3, the gg — > A/H (di- 
rect production mode - Fig. [[]) rates are dominant and significantly larger than in the SM case. For large 
values of tan/3, the production is dominated by gg,qq—>hh A/H, gg— >bb A/H is largely dominant 
between the two at the LHC, (associated bb production mode- Fig. ||]). 

For niA> 150 GeV, the H and A bosons are degenerated in mass, so the signal rates in the tt- 
channel can be added, whereas a more complicated procedure depending on the experimental resolution 
and on the mass difference m^-mA has to be applied for niA< 150 GeV [ ]149p . Higgs-boson masses 
below 120 GeV have not been considered in this channel because of the large resonant background from 
Z — » tt decays. 

For high tan [3 values A and H couple dominantly to the heaviest lepton and to the heaviest down- 
type quark; the branching ratio of A/H is ~ 90% into bb and ~ 10% into tt. 

Including the decay of the r leptons, the three possible final states are: 

• the lepton-lepton (e/i) channel, with a branching ratio (BR) of 6.3%; 

• the lepton-hadron channel, with BR = 46%; 

• the hadron-hadron channel, with BR = 41%. 

The lepton-lepton channel has been studied in the low tua region (<~500 GeV) both in ATLAS 



[149] and CMS [166]: compared to the lepton-hadron channel it turns out to provide a worse sensitivity 




g 

Fig. 1 : Direct A/H production mode. 
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to a possible signal, due to both its lower rate and the less favourable kinematics of the r-decay. CMS is 
studying the possibility to use the impact information to reduce the backgrounds in this channel. 



The lepton-hadron channel has been studied in the low to a region both in ATLAS [ |12| ] Q150| ] and 



CMS [169]. In this channel the application of r-jet identification strongly reduces the jet-background 



from various sources. In ATLAS the toa region studied has been recently extended to ~1 TeV [150] 
with promising results. 



The hadron-hadron channel has been studied in CMS [162] and recently in ATLAS [151] in the 
higher toa region (>~ 500 GeV). The leptonic decay channels include a trigger lepton which allows for 
an efficient background reduction; in the hadron-hadron channel the purely hadronic final states compete 
with QCD jets, so it is difficult to maintain the trigger rates acceptable and it is also difficult to find 
criteria to reduce the huge QCD background. To exploit fully the 2 r-jet final states, especially in the 
very low (~ 200 GeV) mass range, an efficient hadronic r trigger has been developed in CMS based on 
Level- 1 calorimeter selection, Level-2 electromagnetic calorimeter isolation [|163|] and a Level-3 tracking 



(isolation) [ ]164| ]. 

The search strategy for all channels is based on kinematical cuts, T-jet identification (for the chan- 
nels where at least one r decays hadronically) and the recontruction of the tt invariant mass m TT , 
so it relies on two very important detector performance requirements. One important feature of the 
A/H — > tt analysis is in fact the possibility to reconstruct the invariant tt mass. The energies of the 
two t's are evaluated from the energies of the r decay products, assuming that they have the same direc- 
tion of the T-parent; the neutrino energies are obtained solving a system containing the two p™ tss com- 
ponents. Therefore, it is crucial to have a very good p™ tss resolution. For the channels where at least 
one t decays to hadrons, a very good T-jet identication is also crucial, to have the possibility to reject the 
huge jet-background from different sources. 

2. ATLAS Results 

2.1 Event Generation, A/H Production Cross-Sections and Branching Ratios to tt 

The signal and background events were generated with the PYTHIA 6.152 Monte Carlo event generator. 
The CTEQ5L parametrisation of the structure functions was used. The fast ATLAS detector simulation 



was used [138] 



The direct A/H production (from gg — > A/H — > tt) cross-section is calculated using the program 



HIGLU [[155]], based on the results of []17|]. The associated bb A/H cross-section is calculated using 
the program HQQ, which calculates the production cross-section of Higgs bosons via gg, qq—>bb A/H 
according to the results presented in [ ]156| ]. The MSSM Higgs sector is implemented in the approximate 



two-loop RGE approach of [73]. The program HDECAY is used to calculate the total decay widths 



and the branching ratios [8C]. The cross-sections calculated at leading order (LO) are used here. The 



direct A+H production cross-sections and the associated bb A+bb H cross-sections, both multiplied by 
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m A (GeV) 

Fig. 3: Direct A/H production cross-section times 
BR( A/H—>tt) as a function of tua for different 
tanp values. 
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Fig. 4: Associated bb A/H production cross-section times 
BR(A/H — >tt) as a function of tua for different 
tanfl values. 



the BR(A/ H— ►tt), are shown in Fig. || and in Fig. § respectively, for three different tanf3 values. For 
large values of tan[3, the production is dominated by the associated production mode, moreover, for a 
fixed tan/3 value, the ratio between the associated and direct production increases as tua increases. 
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Fig. 5: Average pt(A/H) as a function of tha in direct and associated production in PYTHIA6. 1 

There are differences in the event topology and kinematics between the events from the two dif- 
ferent production processes: 

• there are b-jets in bb A/H events 

• the px distribution of the generated A/H is different. As can be seen from Fig. ||, the average pxA. 
is larger in direct production for masses larger than 150 GeV. 

It must be underlined that the theoretical uncertainty for both the computation of cross-section and 



simulation of events for the bb associated production is still large [157]. 
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2.2 r-jet identification, p r rp' tss resolution and the reconstructed m TT resolution in ATLAS 

An excellent r-jet identification performance to suppress the huge jet-background from various sources 
is necessary for the A/H — » tt study for the channels where at least one r decays hadronically. In 
ATLAS the T-jet identification/jet-rejection has been studied with full detector simulations of signal and 
background events [ 152| |. The criteria to identify a hadronic jet as a r-jet are based on both calorimeter 
and tracker information. 

In the px region 30-150 GeV, the requests are (here called TDR criteria because they were used to 



obtain the results reported in in the ATLAS Physics Performance Technical Design Report (TDR) [|12|]): 

• Rem < 0.07, where R em is the jet radius computed using only the e.m. cells contained in the jet; 

• AE^j? < 0.1, where AE)j? is the difference between the transverse energies contained in cones of 
size AR = 0.2 and 0.1, normalised to the total jet transverse energy Et\ 

• N tr = 1, where N tr is the number of reconstructed charged tracks with pr > 2 GeV pointing to the 
cluster. 

With these cuts the r-jet identification efficiency (e T ) is ~ 25% for r's from m,A=150 GeV (< 
p T T ~ jet >~ 50 GeV) and the jet - rejection goes from ~ 170 to ~ 1700 for jets in 30<pr<150 GeV, 
depending on the pr and on the jet type (light quark, gluon, b-jet). With the same criteria an e T ~ 40% 
for m A =800 GeV (< p T T ~ jet >~ 200 GeV) and a jet - rejection ~ 2500 against QCD jets with 
Pt> 150 GeV can be achieved. In Fig. || the T-jet identification efficiency as a function of niA is shown, 
while in Fig. [7] the jet-efficiency is shown as a function of the jet pr- The r-jet identification criteria 
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Fig. 6: Efficiency of the r-identification criteria for r's 
from A — > tt decays as a function of iba. 

were optimized for the high ttia and p^(>150 GeV) [|150| combining the Atlas TDR criteria and the 
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Fig. 7: Efficiency of the r-identification criteria in QCD jet 
events as a function of p^ 1 . 



criteria used by CMS in [ |162[ ]. Asking for: 

• R em < 0.12 

• CMS cuts 

- 1 isolated track with px>40 GeV within Ai?<0.1 from the jet axis 

- track isolation: no other track with pr>l GeV in a cone of AR=0.4 

an e T ~ 45% (e T ~ 55% - only CMS criteria) for r's from m-A=800 GeV can be achieved and the 
jet-rejection can be significantly improved to ~ 3500 ( ~ 2000 - only CMS criteria) for QCD jets with 
p T >150 GeV. 

However, the jet-efficiencies are still determined with large errors: a larger statistics of fully sim- 
ulated events containing jets is still needed to complete this study. 
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Moreover a very good E™ lss -resolution performance for the reconstruction of the rr mass is 



required for the A — > rr channel study. Crucial for a good p™ ss resolution Q153[ ] are the calorimeter 
coverage until |r/|<5, p™ lss and p™ tss have to be reconstructed from all calorimeters cells (in clusters 
and outside the clusters), a careful calorimeter calibration and intercalibration is necessary and finally a 
careful choice of electronic noise cutoff has to be made (only cells with Et> l.5a(noise) are kept). 

In this way, the p™ tss resolution in ATLAS is found to be: 

a(p™ iss ) = 0.46 * y/EEr 

where T>Et is the total transverse energy in the calorimeters expressed in GeV. 

This formula is valid at low luminosity and it takes into account both the coverage effect and the 
energy resolution. At high luminosity, there is a strong degradation of the p™ lss resolution (about a 
factor of 2 worse) due to the pile-up as described in the TDR [O]. 



The invariant mass of the r-pair in A/H — ► rr can be reconstructed in the collinear approxima- 
tion that the directions of the two neutrino systems from each r-decay coincide with the ones of the 
measured r-decay products and under the condition that the r-decay products are not back-to-back: 



v/2pi + E Ul )(E 2 + E U2 )(l - cos9) 



where E\, E<i are the energies of the measured r-decay products, E Vl , E V2 are the energies of the two 
neutrino systems and 6 is the angle between the directions of the measured r-decay products. E Ul 
and E U2 are obtained by solving a system containing the two p™ lss components. The measurement 
uncertainties on p™ lss , p™ lss combined with the assumption on the directions of the decay-products 
often result in unphysical negative solutions for the neutrino energies, in that case the A mass cannot be 
reconstructed. 

The mass resolution a{m TT ) is proportional to a(p™ lss ) /\sin(A(f)(pip2))\, therefore both the 
p-miss reso i u tion and the A</> separation between the charged r-decay products are important in the rr 
mass reconstruction [154]. 

The reconstructed m TT resolution has been compared for the three different final states (after 
having applied a cut on the lepton pt>24 GeV and on the r-jet p?>40 GeV - the dependence of the 
mass resolution on the px cutoffs is weak - and the cut Atfi < 165°) (see Fig. @) and it is found to 
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be comparable, as it is expected on the base of the p™ 
calorimeters, shown in Fig. 



resolution, shown in Fig. |9[ and the T,Et in 
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Fig. 9: Pt" s resolution in A events as a function of A Fig. 10: E £ T in calorimeters in A events as a function of 

mass. A mass. 

The reconstructed mass resolution is systematically worse for the associated production events. 
In the associated production events in fact the average px of the generated A is lower with respect to 
the direct A production and the difference increases with m A (see Fig. |5|). This implies that in bb A 
events the two r's from the A tend to be more back to back, with two important consequences that are a 
lower acceptance of the A(j) cut and a final worse solution of the system giving the neutrino energies and 
consequently a lower efficiency in the mass reconstruction and a worse resolution. At m A =450 GeV the 
relative a(m TT ) is ~ 10% and ~ 13% respectively for direct and associated production with an overall 
efficiency of the A0 cut and of the request of positive solutions for the neutrino energies of ~ 30% and 
~ 20%. 



2.3 The lepton-hadron channel analysis in ATLAS 

The irreducible backgrounds for this channel are tt —>bW + bW~—>bb lep r and Z/7* — > tt , the re- 
ducible backgrounds are tt -^bW + ~b~W~ ^hh lep had, W—>lep +jets and bb — >lep had. 

The analysis criteria (Standard A analysis) are: 
, pT ie P ton >24GeV(40GeVform A >500GeV)and|7 ? | /e P ton < 2.5; 

• Isolation of the trigger lepton (which rejects leptons from bb by a factor 100 for a 90% efficiency for 
isolated leptons); 

• mrilepton — p™ lss ) < 25 GeV (against the backgrounds containing W), where my(/epton — pjp tss 
is the transverse mass of the lepton-neutrino system; 

• p^ iss > 18 GeV (40 GeV for m A >500 GeV); 

• E T i et > 40 GeV (80 GeV for m A >500 GeV), |r/p e * < 2.5 (r-Candidate); 

• Acj)(jet - lepton) in 100° - 165°; 

• m rr in the window m A ± AM (AM = 1.5<7 mTT ). 



Each event is weighted using the r-jet identification factor for the r-jet candidate (see section |2.2| ). 

Due to the topological and kinematical differences in the direct and associated events, two different 
analyses, one optimized for the direct production process, the other one optimized for the associated 
production are performed with the following criteria: 
- Direct analysis: 
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• zero b-jet tagged (against tt and bb backgrounds); 

• all cuts of the Standard A analysis (r-jet identification, kinematic and mass cuts). 
- Associated analysis: 

• 1 b-jet tagged (against Z and W+jets backgrounds); 

• number of non b-jets < 3 (against tt backgrounds); 

• cuts of Standard A analysis (r-jet identification, kinematic and mass cuts) except the cut on A<j)(jet — 
lepton), to not reduce too much the signal acceptance. 

To choose the b-jets, a b-tagging efficiency of 60%, with a corresponding rejection of 100 against other 
jets and of 10 against c-jets has been randomly applied. 

The dominant background selected by the direct analysis arises from W+jets, which have the 
largest production cross-section and from the Z — » tt at the lower masses. The analysis optimized for 
the associated production rejects much better W+jets and Z — » tt backgrounds and the tt background 
becomes dominant. 

Having the opposite request to have or not to have a b-jet tagged, the two analyses are not corre- 
lated, so, after having applied them separately to both signal samples (direct and associated A production) 
and to background events, the significances can be combined. 

The results at lower masses have been compared to the results reported in the TDR and they 
have been found to be in reasonable agreement [ |150( ]. The differences observed are due to the different 
PYTHIA version, to the use of fast instead full simulation and to the use of different cross-sections values 
for signals and backgrounds. 



At the higher masses (ma > 500 GeV), the analysis is performed only on the events from the 
bb A production channel, due to the complete dominance of that production mode (see section 2.1). 

Despite the low production cross-section and the low acceptances of the analysis (~ 0.5% for 
mA=800 GeV) the backgrounds are strongly reduced (the total background for m,A=800 GeV is ~4.6 
events, dominated by tt background, in 30 fb _1 ). 



Figure 1 1 shows the distribution of m TT after the analysis cuts (except the mass cut) for associated 
signal events at ttia=800 GeV for tanf3=45 and for the main backgrounds normalized to the expected 
event number for an integrated luminosity of 30 fb _1 . 
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Fig. 11: Reconstructed m TT in the lepton-hadron channel after the analysis cuts for mA=800 GeV for tan/3=45 and for the 
tt and W+jets backgrounds (plotted separately) assuming an integrated luminosity of 30 fb~ . 
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2.4 The hadron-hadron channel in ATLAS: trigger study and analysis 

Due to the very low production cross-section of the A/H especially in the higher mass region between 
0.5 and 1 TeV, the trigger acceptance for the signal needs to be excellent; on the other hand, the signature 
of the signal with only hadronic decays is similar to QCD-background. 

In the first level trigger, the total rate of Jet and r triggers is dominated by 2-jet QCD-background, which 
has to be controlled by appropriately high settings of the Ex thresholds for jet triggers [ ]161[ |. However, 
these settings also reduce the acceptance of the trigger for the hadron/hadron channel of the A/H decay. 
A high input acceptance of the trigger for A — > rr — > hadron — hadron signal events can be achieved 
by using combined ]e,t+E™ lss and T+EJp iss triggers. 

The fast trigger simulation ATL1CT [158] [159] interfaced to the fast ATLAS simulation program 
has been used for the study. The following trigger types are of relevance for the hadron4iadron channel: 

• Jet + EJp iss (Threshold: E J T et > 50GeV, E™ iss > 50GeV) 

• t + EJp iss (Threshold: E T > 20GeV, E™ iss > 30GeV) 

• Single Jet (E J T et > 180GeV) 

• Three Jet (E J T et > 75GeV) 

• Four Jet (E J T et > 55GeV) 

The threshold settings on the cluster Et of the jets and r's are set to accept 90 % of the jets with the Et 
value given in the trigger menu, which is identified with the Et value in the reconstruction [159]. The 
isolation criteria of the r trigger are set to fixed values of 2 GeV for the electromagnetic and 4 GeV for 
hadronic trigger towers [160]. 

The total acceptance in the Level- 1 trigger is determined by the number of events accepted by at least 
one of the trigger types. This is equivalent to a logical OR. Many events are accepted by more than one 
trigger. 

The input acceptance of the Level- 1 trigger for signal events rriA/ H = 450,600 and 800 GeV 



for the individual trigger types listed above and their combination (OR) are shown in Fig. 12. For 
m A/H = 800 GeV, an input acceptance of 76.6 % is reached using fixed isolation thresholds in the 
r trigger. Using an alternative scheme of a dynamic isolation in the r trigger (electromagnetic trigger 
towers: 4 % of cluster energy, hadronic: 8 % of cluster energy), this value can be improved to 80.6 
% [151]. The total trigger rate was evaluated being ~ 1.4 kHz, which fits well within the limitations of 
the Level- 1 trigger menu [160]. 

The influence of the trigger acceptance on the discovery contour, however, is given by the combination 
of the trigger acceptance and the offline analysis, discussed after. 

Since the associated production is dominant for high masses, the expected signal event signature of 
the hadron-hadron channel consists of two high-p-r r's in hadronic decay with two b-jets. Backgrounds 
for this channel are 2-jet QCD, tt , W+jets and Z+jets events. The analysis uses event weighting, where 
the two jets with the highest px are considered as 'r candidates', while b-tagging is used for all other jets, 
where one of the two expected b-jets is requested to be tagged. Using this method, background events 
are not rejected by the requests on the number of r's and b-jets, but weighted accordingly, therefore 
background rate estimates can also be given for channels with extremely low acceptance, especially 2-jet 
background from QCD and also tt . The r identification described in Sec. 12 is used to derive the 
r acceptance factors for each r candidate. The r identification is assumed to be e T = 55% and the 
corresponding jet rejection is used for the other jets. The b-jet tagging efficiency is here assumed being 
n, = 70%. 

The following cuts have been applied: 

• Two jets in the event with pt > 100 GeV and \rj\ < 2.5 (r candidates); 

• No lepton (e, n) with p T > 10 GeV; 

• Not more than 4 jets in |jy| < 3.2 with p T > 20 GeV; 

• at least one b-jet tagged ; 
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Fig. 12: Level-1 trigger input acceptance for Signal m^/H = 450, 600 and 800 GeV: total acceptance by all Jet and r-triggers 
(OR), and acceptances of each separate trigger. 



and ?tit(t2,p™ ss )); 



• p^ lss > 65 GeV; 

• A<f> between r candidates in 145° — 175°; 

• transverse mass cut: mx < 50 GeV (minimum of mr(Ti,p^ 8 

• m rr in the window toa ± AM (AM = l.5a mTT ). 

Some cuts lower the acceptance of signal events significantly, especially the p™ lss cut and b-tagging, 
however, those cuts are necessary to suppress the background channels efficiently The acceptance of the 
analysis cuts for m,A=800 GeV is ~ 0.6% and the total background is ~ 5.4 events, with dominance of 
tt (2.2 events) followed by Z (0.8 events) and 2-jet QCD background in 30 ftr 1 JBlj ]. 



Figure |13] shows the distribution of m TT after the cuts (except the mass cut) for the events at 
?7ja=800 GeV for tan(3=50 and for the main backgrounds normalized to the expected event number for 
an integrated luminosity of 30 fb _1 . 

In Fig. 



14 the reconstructed m TT is shown for ttia=600 GeV in the A — > rr — > hadron — hadron full 



simulated events, after the pt J6< > 50 GeV and Ac/) < 165° cutoffs. 



The combined acceptance of the Level- 1 trigger and offline analysis for the hadron-hadron channel 
have been studied. The acceptance after all cuts for events passing the trigger conditions are evaluated 
to be 92.5%, 95.6% and 95.3% respectively for m A =450, 600 and 800 GeV. Especially for high m A , the 
acceptance is very good due to the similarity of the kinematic cuts used in the Level- 1 trigger and offline 
analysis; both require large values of E™ lss and jets/r's with high px- 



2.5 Combinations of the lepton-hadron and the hadron-hadron channels in ATLAS 

For higher mass values, the results from the lepton-hadron and from the hadron-hadron channels can be 
combined to improve the signal significance. 

In Table [IJ, the discovery tanf3 values (giving a 5a significance) are reported for the two channels 
separately and combined. The extended discovery contour plot is given in Fig. [15| in the tanj3 logarith- 
mic scale and with a linear tan(3 scale in Fig. [l6[ 



3. CMS Results 
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for the tt and Z backgrounds assuming an integrated lumi- simulation, 
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Table 1: Discovery tan/3 values in ATLAS (5<r confidence). 
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22.1 


19.0 


600 


32.9 


30.0 


25.2 


800 


50.0 


45.0 


41.4 




100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000 
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Fig. 15: Discovery contour curves for the Fig. 16: Discovery contour curves for the 
A/H —> tt channel in the (rriA,tanf3) plane in AT- A/H — > tt channel in the (niA,tanf3) plane in AT- 
LAS for an integrated luminosity of 30 fb -1 . LAS for an integrated luminosity of 30 fb -1 shown using a 

linear scale for tan/3. 
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Fig. 17: Efficiency of QCD background and Higgs signal at High Level trigger when the size of the pixel isolation cone is 
varied. 
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* IV 



3.1 ^4, — > tt — ► 2 t jets in CMS 

A, H — > rr with 2 r jet hadronic final states have been shown to extend significantly the SUSY Higgs 



discovery reach into the large mass (600 - 800 GeV) range [(162J]. To exploit fully the 2 r jet final states 
- especially in the low (~ 200 GeV) mass range - an efficient hadronic r trigger has been developed 
based on Level- 1 calorimeter selection, Level-2 electromagnetic calorimeter isolation [|163|] and a Level- 



3 tracking (isolation) using only the pixel detector information [ ]164| ]. 

Level- 1 calorimeter single or double Tau trigger with thresholds of 80 and 65 GeV for L = 2 x 
10 33 cm -2 ^ 1 selects A,H—>tt—> 2r-jet events useful for off-line analysis with an efficiency of 
about 0.9 while giving an output QCD background rate of about 6 kHz. A further reduction of the QCD 
background rate by a factor ~ 10 3 is possible at the High Level trigger path (Level-2 calorimeter and 
Level-3 Pixels) with an efficiency of ~ 40% for the signal at mu = 200 and 500 GeV [165] as one can 
see in Figure |l7|. Even better performance is expected using the regional tracking option of the CMS 
High Level trigger once the CPU performance is proven to be satisfactory. 

For the off-line r identification the tracker information is used. The fast simulation of the CMS 



detector [|135J is used to study the signal to background ratios. The track reconstruction efficiency eval- 
uated with full simulation of CMS tracker is included as a function of pt and rj for the track. The r jet 
candidate (E t > 60 GeV) is required to contain a hard (p t >40 GeV) charged track within AR < 0. 1 
around the calorimeter jet axis. Around this leading track in a cone of AR < 0.03 two other tracks with 
Pt > 1 GeV are accepted to include the 3-prong r decays. This narrow cone with one or three hard tracks 
is required to be isolated demanding that no track with p t > 1 GeV is found in the surrounding larger 
cone of AR < 0.4. The efficiency for this r selection is 7.2% for rriA = 200 GeV and 34% for itia = 
500 GeV. Accepting the 3-prong decays in the narrow cone of AR < 0.03 increases the event rate for 
A, H — > tt — ► 2 r jets in the high mass range ( tua = 500 GeV) by ~ 1.7 but also degrade significantly 



the QCD rejection factor for hard QCD jets. Figure |18| shows the rejection factor against the QCD jets for 
the 1/3 prong selection a function of E t jet compared to the one prong selection with one hard (p t >40 
GeV) charged track within AR(jet, track) < 0.1. Optimization is still needed for the low mass range 
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for more efficient selection mainly by increasing the size of the narrow cone. 

A further suppression can be obtained exploiting the r lifetime using a r vertex reconstruction 
or impact parameter measurement or a combination of them. A full simulation study indicates that an 
additional rejection factor of ~ 5 against the 3-prong QCD jets and an efficiency of ~ 70% for the r jets 



can be obtained with r vertex reconstruction [ |167| ]. Promising results are also obtained from the impact 
parameter method in the channel A, H —> tt —> £ + £~ + X using full simulation combining the impact 
parameter measurements for the two leptons from r decays to reduce the backgrounds with W — ► Iv 
and Z -> ^ decays [|66|]. 
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Fig. 18: Rejection factor against the QCD jets as a function pig. 19: Reconstructed Higgs mass for bb, H -> rr -> 
of £ t jet for the 1 and 1/3 prong r selection. 2 r jets with m H = 500 GeV. 

The resolution of the reconstructed Higgs mass and even more so the mass reconstruction effi- 
ciency in A, H — > rr events is very sensitive to the £™ ss measurement. The absolute value of EY 11813 
is relatively small in these events making the mass reconstruction and background reduction with a cut 

in E^ iss 



a difficult task. Figure 19 shows the Higgs mass reconstructed with full simulation for bbA, 



A — ► rr — > 2 r jets with rriA = 500 GeV and tan(3= 20 [ 168 ]. The resolution of the reconstructed Higgs 
mass is 14.5% for ttla = 200 GeV and 14.9% for tua = 500 GeV and the corresponding reconstruction 
efficiencies are 37% and 36%, respectively (including Acf) < 175° cut and requiring positive neutrino 



energies). This confirms the earlier results of the fast simulation study [ ]162| ]. 

The large Z, 7* — ► rr background can be reduced efficiently only with b-tagging in the associated 
production processes bbHsusY- The associated b-jets are soft and uniformly distributed over < 2.5. 
Nevertheless, a study with full simulation shows that a b-tagging efficiency of ~ 34% per jet can be 



obtained for the signal events with a mistagging rate less than 1% for Z + jets events [ |166| |. Requiring 
one tagged b-jet reduces efficiently also the QCD background thus improving significantly the signal 



visibility. Figure 20 shows the signal for m# = 500 GeV and tan(3 = 25 superimposed on the total 
background with b-tagging. The missing transverse energy plays a major role in the Higgs mass 

reconstruction as discussed above. However, a cut in E™ tss does not improve significantly the mass 



resolution and therefore, in order to retain the signal statistics, is not used in this study. Figure [21] shows 
the expected discovery reach for for 30 fb^ 1 assuming a maximal stop mixing scenario [70, 79]. The 



expectations for other important MSSM Higgs discovery channels in CMS [|170J and the exclusion region 
from LEPII [Q] are also shown in the figure. The A, H — > rr channels are found to be insensitive for 
stop mixing, the SUSY scale and for the sign of the Higgsino mass parameter \x for high enough tan(3 
( >, 10). A systematic study of the A,H — > rr with £ + £~, lepton + r jet and 2 r jet final states 
is presently in progress in CMS including full simulation of the hadronic r trigger, r identification, r 
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tagging with impact parameter and vertex reconstruction, Higgs mass reconstruction and b-tagging in 
the associated production channels. 




Fig. 20: Higgs mass for H — > tt — > 2 r jets with mn Fig. 21: Expected 5a discovery reach for the MSSM Higgs 
= 500 GeV and tan(3=25 superimposed on the total back- bosons in CMS in the maximal mixing scenario for 30/fe -1 
ground for 30 fb^ 1 . One tagged b-jet is required. as a function of mA and tanfi. The shaded area is excluded 

by LEP [fjj. 

4. Conclusions 

The LHC discovery potential for A/H — » tt has been studied in ATLAS and CMS in the three different 
final decay channels at low luminosity (10 33 cm s~ ) in the range from ~100 GeV until ~1 TeV. 

• at lower masses (<~500 GeV), the lepton-lepton and lepton-hadron channels have been studied in 
both experiments; the lepton-hadron gives the best sensitivity. To study the hadron-hadron channel 
also in the lower mass range, CMS is developing an efficient hadronic r trigger (Level- 1 + High 
Level trigger). 

• at higher masses, for large tan/3 values the A production is dominated by bb A. In this mass 
range the hadron-hadron channel can be studied because it is possible to reject the huge QCD 
background with kinematical cuts and T-jet identification. The trigger is a very important point for 
this purely hadronic channel and it has been studied in the two experiments. The lepton-hadron 
channel has been also studied in ATLAS in the high toa region giving promising results; the 
combination of the hadron-hadron and lepton-hadron channels improves the discovery potential. 

The expected 5<r-discovery contour curves for the combined Aj H — > tt signal show that a signal should 
be observed over a large region of the {m^tanfS) plane, with mA up to ~ 1 TeV assuming an integrated 
luminosity of 30 fb _1 . 

At high mA (~0.5-~l TeV) the A/H — ► tt channel should be observable for tanj3 values greater 
than ~ 25. This is an important result because the A/H — > tt channel is up to now the only one giving 
access to this high mA region. 

Studies are still in progress, both on the experimental and on the theoretical sides, to improve the 
results presented here. 
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F. Searching for Higgs Bosons in ttH Production 



V. Drollinger 



Abstract 

Higgs boson production in association with ti pairs with the subsequent de- 
cay into l^vqqbbbb is analyzed including all relevant background processes. 
Excellent 6-tagging performance and a good mass resolution turn out to be 
the most important components for a successful analysis. The top Yukawa 
coupling can be determined with an accuracy of about 17% in this process, 
provided the branching ratio of H — ► bb is known with a sufficient accuracy. 
Finally, a first estimate of the potential size of higher order corrections to the 
ttbb background is given. 



1. Introduction 

If the Higgs boson is lighter than 130 GeV/c 2 , it decays mainly to a bb pair. To observe the Higgs boson 
at the LHC, the ttH channel turns out to be the most promising channel among the Higgs production 
channels with H° — > bb decay [171]. In this study, wediscuss the channel ttH — ► l^vqqbbbb (Figure |]), 
where the Higgs Boson decays to bb, one top quark decays hadronically and the second one leptonically. 
The relevant signal and background cross sections at the LHC (y / s^ = 14 TeV) and particle masses 
used in the simulation are listed in Table |T[ This is the first set of signal and background processes 



Table 1 : CompHEP \92J cross sections for signal and background relevant for the tiH° — ► vqqbbbb channel, calculated with 
parton density function CTEQ41 ^172\ j. The branching ratio of the semileptonic decay mode (one W ± decays to quarks the 
other W decays leptonically, where only decays to electrons or muons are taken into account) is 29% (not included in the 
cross sections of this table) and m w ± = 80.3427 GeV/c 2 . 



LO cross sections 


masses 


a tm o x BR H o^ b - b = 1.09- 0.32 pb 
a tiz o = 0.65 pb 
a m = 3.28 pb 
vttjj = 507 pb 


m H o = 100 - 130 GeV/c 2 
m z o = 91.187 GeV/c 2 
m b = 4.62 GeV/c 2 
mt = 175 GeV/c 2 



completely calculated at LO for the ttH channel. The hard processes are generated with CompHEP and 



then interfaced to PYTHIA, where the fragmentation and hadronisation are performed [92]- [100]. The 
combined package CompHEP-interface-PYTHIA includes all features of a pure PYTHIA simulation, 
such as initial sate radiation, final state radiation, multiple interactions and underlying event. After the 
final state has been obtained, the CMS detector response is simulated, with track and jet reconstruction 



with parametrisations FATSIM [ 175] and CMSJET [135], obtaining in this way tracks, jets, leptons (the 
electron or muon reconstruction efficiency is assumed to be 90%; taus are not considered here) and 
missing transverse energy. These parametrisations have been obtained from detailed simulations based 
on GEANT. 



2. Reconstruction 

From Figure [l] we expect to find events with one isolated lepton, missing transverse energy E™ and six 
jets (four 6-jets and two non-fe-jets), but initial and final state radiation are sources of additional jets. So 
the number of jets per event is typically higher than six. On the other hand, not all six quarks of the 
hard process can be always recognised as individual jets in the detector, in which case it is impossible to 
reconstruct the event correctly - even if there are six or more jets. 
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Fig. 1: Example of a pp — > ttH° — > ^uqqbbbb signal event. The LO process is drawn in red. HO events include gluon 
radiation (light green) in addition. The expected final state consists of one isolated lepton, missing transverse energy, four b-jets 
and two (or more) non-b-jets. 

For the reconstruction of resonances it is necessary to assign the n jets of an event to the cor- 
responding quarks of the hard process. In general, and ignoring information on 6-jets, the number of 
possible combinations N is given in Table |2| as a function of the number of jets per event. We obtain 
N for the case, when the masses of the Higgs boson, both top quarks and the hadronically decaying W 
boson are reconstructed. The nominal mass of the leptonically decaying W boson, together with E™ 
and the lepton four momentum, is used to calculate two solutions of the longitudinal momentum of the 
neutrino pz{v) which is needed for the mass reconstruction of the leptonically decaying top. 

Table 2: Number of jets per event n and the corresponding number of possible combinations N. If there are more than a dozen 
jets, only the twelve with highest Et are considered. 





x 6! x 


2 x \ X 2 = (g) 


x 360 




n = 6 7 


8 


9 10 


11 


12 


N = 360 2520 


10080 


30240 75600 


166320 


332640 



Good mass resolution and the identification of 6-jets is essential to reduce the number of wrong 
combinations in the event reconstruction. A good mass resolution can be obtained when the energy and 
direction of each reconstructed jet agree as closely as possible with the quantities of the corresponding 
parent quark. This can be achieved with jet corrections as described in [176] and [177]. For 6-tagging 
we use the ^-probability functions which depend on impact parameters of tracks and leptons inside the 
jets. They are determined using ti six jet events, as described in [171]. The identification of 6-jets is even 
more important for efficient background suppression. 

Figure || shows the invariant mass distributions of the reconstructed resonances of ttH° -> l±vqqbbbb 
events in the case of an ideal reconstruction: after the "preselection" and the calculation of pz(y) (see 
later on) each quark of the hard process is matched with exactly one jet, the closest one in R = \J 4> 2 + rj 2 
if AR(q,j) < 0.3 and if the jet energy is closer than ± 30 % to the parent quark energy. The mean val- 
ues and widths of the top and W mass distributions are used to define likelihood functions used in the 
selection procedure described in the following. 

o Preselection 

Events are selected if there is an isolated lepton (e^ or fi^ with pr > 10 GeV/c within the tracker 
acceptance; no other track with px > 1 GeV/c in a cone of 0.2 around the lepton) and at least six jets 
(Et > 20 GeV , \r]\ < 2.5). 
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Fig. 2: Invariant resonance masses of the tiH° — > ^vqqbbbb signal: Higgs boson, leptonic top, hadronic top and hadronic 
W . The leptonic W is not reconstructed but its nominal mass is used to calculate pz(f)- The generated masses are: m H o = 
115GeV/c 2 , mt = 175GeV/c 2 andm w ± = 80.3427 GeV/c 2 . 



o Event Configuration 

In order to be able to reconstruct the Higgs mass, we have to rind the correct event configuration among 
all possible combinations listed in Table ^[ The best configuration is defined as the one which gives the 
highest value of an event likelihood function ([[]) which takes into account 6-tagging of four jets, anti-6- 
tagging of the two jets supposed to come from the hadronic W^, mass reconstruction of and the 
two top quarks, and sorting of the 6-jet energies. 

L.EVNT = J] P.ih) x H [1 - P b ( qi )\ x n e" - 5 *^^ x f[E b (t,F) - E b (H°)\ (1) 

i=l,4 i=l,2 i=W ± ,t,t 

o Jet Combinations 

Events with more than six jets can contain gluon jets from final state radiation, which are not yet used 
in the analysis. The combination of these jets with the correct quark jets can improve the event recon- 
struction further. The additional jets are combined with the decay products of both top quarks if they 
are closer than AR(j, j) < 1.7, if the corresponding mass is closer to the expected value of Figure ||| If 
there are still jets left, they are considered as Higgs decay products and are combined with the closest of 
the corresponding two fe-jets, if AR(j,j) < 0.4. 

o Event Selection 

Three likelihood functions: for resonances (LJIESO > 0.05), 6-tagging (L.BTAG > 0.50), and 
kinematics (L_KINE > 0.2) are used to reduce the fraction of background events. Finally, the events 
are counted in a mass window around the expected Higgs mass peak (mi nv (j, j) in fh ± 1.9 a ; fh and a 
are obtained from mass distributions as shown in Figure || with various generated Higgs masses). 

The overall efficiency for a triggered event to be finally selected is 1.3% for ttH° (m H o = 115 
GeV/c 2 ), 0.2% for tiZ°, 0.4% for ttbb and 0.003% for tijj events. This shows that the reducible 
background is reduced very effectively. In addition, there is little combinatorial background left (see 
Figure |3|) with this reconstruction method. 
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Fig. 3: Simulated invariant mass distribution of signal (dark shaded, m H o —115 GeV/c 2 ) plus background for L int — 30 
fb~ . The dashed curve is obtained from the fit of the background without signal, the solid line describes the fit of signal plus 
background. The small dashed line corresponds to the LO background contribution. 

3. SM Results 

After the whole reconstruction and event selection procedure, it turns out that the irreducible background 
(with four real 6-jets) is dominant. Even the tijj background, where only two 6-jets from the top decays 
are generated in the hard process, is dominated by events with four real 6-jets. This is possible after the 
fragmentation of PYTHIA: e.g. gg — ► tigg — ► l^vqqbbgbb with one bb pair coming from g — > bb (gluon 
splitting). In this case the final state consists of nine partons or leptons which is one more than expected 
at LO and is therefore considered as HO (in this case NLO) process. Together with the number of tibb 
events (considered as LO) we obtain an intrinsic k-factor k t j q q = 1.9 for all tiqq events as indicated in 
Figure [| (For the full m H o =115 GeV/c 2 selection we get 23 events from "ttbb" and 20 events from 
"tijj" with four real 6-jets plus 6 events with two real and two false 6-tags. The total number for the 
non resonant background amounts to 49 events, whereas the corresponding number from the PYTHIA ti 
process with additional jets from fragmentation is only 24 events.) In case of the tiH° signal additional b- 
jets from fragmentation cannot enhance the signal, but rather complicate the reconstruction of the correct 
invariant mass, the Higgs mass, in the end. 

The signal to background ratio S/B, the significance S/VB for Lj nf = 30 fb^ 1 , the integrated 
luminosity Li n t required for a significance of five or more and the precision on the top Higgs Yukawa 
coupling y t for Li nt = 30 fb^ 1 are shown in Figure |] as a function of the generated Higgs mass: S/B 
is around 50% and an relatively low integrated luminosity is sufficient to discover the Higgs boson in 
this channel with a significance above five. An integrated luminosity Li nt = 100 fb^ 1 would be enough 
to explore all points considered in Figure |] up to a Higgs mass of 130 GeV/c 2 . If we assume a known 
branching fraction of the decay H° — ► bb, it is possible to determine the precision of y t with accuracy 
of about 17%. Apart from these results, the Higgs mass can be determined from the Gaussian fit of the 
final mass distribution (see Figure ||) with a precision of better than 6% for Li nt = 30 fb~~ l . 

4. Conclusions 

From our present understanding, it is experimentally possible to observe the tiH° — ► l^vqqbbbb channel. 
Most important for a successful analysis are excellent 6-tagging performance and a good mass resolution. 
The effects of event pile up still have to be investigated. 

From theoretical point of view, the first complete LO simulation has been performed for signal and 
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Fig. 4: S/B, S/VB, L int (required for S / y/B = 5)andAy t /yt versus generated Higgs mass in the SM. All results are based 
on CompHEP cross sections calculated at LO; from the fragmentation with PYTHIA an intrinsic k-factor k t t q q = 1.9 for tiqq 
background events is included. 



background. After the K-factor for the signal has been calculated (see Ref. [51]), the main uncertainty 
in this channel is background cross section at HO, because our first estimate gives a factor of almost two. 
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A Comparison of ttbb Events 

Table 3: CompHEP (ISR and FSR included) PYTHIA (default) comparison of ttbb background: the cross sections do not 
include any branching fractions (all decays are allowed). All four fr-quarks are required to be within \q\ < 2.5 and additional 
the 6-quarks of the top decays have to satisfy pribtop) > 15 GeV/c. The pT(b g i u ) cut for both 6-quarks not coming from the 
top decay is varied. 



Pt Cuts 


CompHEP 


PYTHIA 


CompHEP / PYTHIA 


Pr(b g iu) > 15 GeV/c 


a = 2407 fb 


a = 2927 fb 


0.82 


Mbglu) >30 GeV/c 


a = 1123 fb 


a = 1189 fb 


0.94 


p T (b g iu) >50 GeV/c 


a = 512 fb 


a = 431 fb 


1.19 


PT(bgiu) > 100 GeV/c 


a = 116 fb 


a = 53 fb 


2.19 


p T {bgiu) > 200 GeV/c 


(7= 13/6 


a = 2/6 


6.73 


p T (bgiu) > 300 GeV/c 


(7 = 3/6 


£7 = 0/6 
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G. Studies of Charged Higgs Boson Signals for the Tevatron and the LHC 

K.A. Assamagan, M. Bisset, Y. Coadou, A.K. Datta, A. Deandrea, A. Djouadi, M. Guchait, Y. Mambrini, 
F. Moortgat and S. Moretti 



Abstract 

Two Higgs doublet models are a viable extension to the Standard Model (SM) 
and can be incorporated into supersymmetry (SUSY). In such models, elec- 
troweak symmetry breaking leads to five Higgs particles, three neutral and a 
charged pair. We discuss various analyzes of the charged Higgs boson, carried 
out in the context of the Minimal Supersymmetric extension of the Standard 
Model (MSSM) and also in models with singlet neutrinos in large extra dimen- 
sions. Specific studies for the Large Hadron Collider (LHC) and the Tevatron 
are presented. 



1. Introduction 

The Higgs sector of the MSSM contains five physical states, two of which are charged, H^, and the other 
three are neutral (h°, H°, and A ) [§, 178]. Searches for the charged Higgs boson have been carried out 
at LEP and at the Tevatron: at LEP 2, a lower bound of 78.6 GeV has been set on the charged Higgs 
boson mass independent of the —* t^v t branching ratio (BR) [179]. At the Tevatron, CDF and D0 
performed direct and indirect searches for the charged Higgs boson, and excluded the low and high tan (5 
regions up to ~ 160 GeV f|l80|l . 

The sensitivity of the ATLAS and CMS detectors at the LHC to the discovery of the charged 
Higgs boson has been investigated in detail [11, 12]. Some of these studies have been carried out as 
particle-level event generation in PYTHIA, HERWIG and ISAJET [[TOOKl8l|,|T8|], at y/s = 14 TeV, with 
the detector resolutions and efficiencies parameterized in ATLFAST [138] and in CMSJET [135] from 
the full detector simulations. Some of the LHC studies assume that the mass scale of supersymmetric 
partners of ordinary matter is above the charged Higgs boson mass so that charged Higgs boson decays 
into supersymmetric partners are forbidden. The main production processes considered in these studies 
are the gluon fusion mechanism, gg — ► tbH^ and the 2^2 process gb — > tH^ shown in Figure [T[ A 
central value of 175 GeV is used for the top-quark mass. The decay channel H ± — ► f^v T has been stud- 






gg^tH + b gb^tH + 

Fig. 1: The charged Higgs boson production at the LHC through the 2^3 process, gg — > tbH ± and the 2^2 
process, gb — > tH ± . The inclusive cross section is the sum of both conttibutions after the subttaction of the 
common terms. 
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ied extensively for ATLAS for m H ± < m t , and the signal appears as an excess of r leptons [183]. The 
channel — > Wh° is only relevant in a tiny range of MSSM parameter space although it constitutes a 
unique test for MSSM and may be sensitive to the singlet extension to MSSM, i.e., NMSSM fl!84| , |185p . 

r ± ^ r are the dominant decay channels of the charged Higgs boson in most of 



H± -> tb and ± - 
the parameter space. In the H ± 

top-quark mass in the low and high tan (3 regions up to 



tb channel, upwards of 5-a discovery can be achieved above the 



A00 GeV [186]. 



T^Vt 



extends the 



discovery reach to high Higgs boson masses and to lower tan values in the high tan (3 region as seen in 
Figure 0. However, in the low tan (3 region, the t v T channel offers no sensitivity for the charged Higgs 
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Fig. 2: The ATLAS 5-a discovery contour of the charged Higgs boson for an integrated luminosity of 300 fb -1 
(left plot); the CMS 5-er discovery contour of the MSSM Higgs bosons for 100 fb _1 of luminosity (right plot). 
Below the top-quark mass, the charged Higgs boson is produced from top decay and the r ± v r channel provides 
coverage for most tan j3 below ^160 GeV. Above the top-quark mass, the tb channel covers the low and the high 
tan (3 regions while the t^u t channel extends the discovery reach to high Higgs boson mass and to lower tan (3 in 
the high tan (3 region. 



boson discovery as the — > r v T branching vanishes [ |187[ ]. Further studies are needed to cover the 
remaining areas of the 5-a discovery contour of Figure § 

• The lack of sensitivity in the intermediate tan [3 region is due to the fact that the charged Higgs 
boson coupling to SM fermions is proportional to 

H + (mt cot l3ibL + mb tan pibji), (1) 



the square of which goes through a minimum at tan (3 = \Jmtjmb- The studies of charged Higgs 
boson production from SUSY cascades and charged Higgs boson decays to SUSY particles might 
help cover this region. 

The gap in the ttia axis around vtla = 160 GeV corresponds to the transition region where, for the 
correct description of the charged Higgs boson production and decay mechanisms, it is mandatory 
to use the production process gg — > tbH^ which includes not only gg — > ti with t — > bH^, but 
also the Higgs-strahlung mechanism and the relative interferences [ ]188j ]. 

The discovery reach could be extended to high Higgs boson masses by studying the process gg — > 



tbH^ with — > tb and tagging all the four b-jets in the spectrum [189]. 

Recent studies which attempt to cover these remaining regions of the parameter space are presented 
along with the observability of charged Higgs boson signals in models with singlet neutrinos in large 
extra dimensions and the prospects for the determination of the charged Higgs boson mass and tan (3 at 
the LHC. 
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2. Mass and tan f3 Determination at the LHC 

In this section, we discuss the expected precisions on the charged Higgs boson mass and tan j3 measure- 
ments at the LHC — above the top-quark mass — in the — > t^v t and — > tb channels. Details 
of this analysis can be found in [ ]190| ]. 



2.1 Motivation 

The detection of a charged Higgs boson signal would constitute an irrefutable proof for physics beyond 
the SM. The subsequent determination of the charged Higgs boson parameters such as the mass, the 
decay width, the spin, the rates in the various decay channels and the couplings to SM and SUSY particles 
will be necessary not only to establish that the observed particle is indeed consistent with a charged scalar 
boson but also to identify the actual scenario that is realized. The measurements of the charged Higgs 
boson mass and tan j3 will be essential to the determination of the charged Higgs boson properties. 



2.2 H ± Mass Determination in H ± 



This channel does not offer the possibility for the observation of a resonance peak above the background, 
only the transverse Higgs boson mass can be reconstructed because of the neutrino in the final state. The 
background comes from single top W t, and ti productions with one — > r v T . Thus, the transverse 
mass is kinematically constrained to be less than the W mass while in the signal the upper bound is the 



charged Higgs boson mass. Furthermore, the distributions of one-prong hadronic decays of t^s, 



TV V T 



(11.1%) 
(25.2%) 
(9.0%), 



(2) 



are sensitive to the polarization state of the r-lepton [ |191| , |192[ ]. In fact, it is to be noted that the spin 
state of t s coming from H^- and -boson decays are opposite (neglecting leptonic mass effects, 
as we did here). This is true for the case of one-prong decays into both tt^'s and longitudinal vector 
mesons, while the transverse component of the latter dilutes the effect and must be somehow eliminated 



by requiring that 80% of the r-jet (transverse) energy is carried away by the vr^s, i.e.: 



R 



p" 
Pt 



> 0.8. 



(3) 



Ultimately, the polarization effect leads to a significantly harder momentum distribution of charged pions 
from r-decays for the ff ± -signal compared to the W /± -background, which can then be exploited to 



increase the signal-to-background ratios and the signal significances [ |187[ , |193| ]. Indeed the background 
is relatively small as shown in Figure where the transverse mass 



mx = y 2pT r p ! T (\ — cos A(f>), 



(4) 



is reconstructed from the visible r-jet and the missing energy. As a result, although there is no resonance 
peak in this channel, the charged Higgs boson mass can be extracted from the transverse mass distribution 
with a relatively good precision. For the mass determination in this channel, we use the likelihood method 



described in [194], which we summarize as follows: 



Suppose we wish to estimate the expected precision 5mo on a Higgs boson reference mass rriQ. 
We generate samples of events with charged Higgs boson masses = uiq + k x 5m and for 
each rrik we calculate the probability density function Pk{m) from the reconstructed transverse 
mass distribution of a charged Higgs boson with mass m^. For example, for a charged Higgs 
boson reference mass tuq = 250 GeV, we generate signal events at charged Higgs boson masses 
m k = 230, 235, 240, 245, 250, 255, 260, 265 and 270 GeV. 
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Fig. 3: The reconstruction of the transverse charged Higgs boson mass in H ± — > r v r for mjj± = 250 and 
500 GeV. The background is relatively small in this channel. The discovery reach is limited to high tan/3 but 
extended to higher mass compared to the tb channel. 



• Assuming iVo is the expected number of events — signal and background — corresponding to 
the reference mass mo, we draw randomly N = Nq + 5 No masses rrij from each distribution 
Pk( m ) ($Nq is the statistical error on No). For each m^, we calculate the likelihood function 
Ck = ^jLi ^og(Pk{rrij)). The differences = Co — Ck show a minimum around mo, where 
a parabolic fit is performed to get the actual expected value of mo- This exercise can be repeated 
many times within the statistical error 5 No and the distribution of the expected values, so obtained, 
of mo would be a Gaussian whose mean is the reconstructed mass and whose standard deviation 
is the statistical precision on the reconstructed mass. 

• Three main sources of systematic uncertainties are included in the mass determination: the shape 
of the background, the background rate and the energy scale. The background shape becomes 
more significant at lower Higgs boson masses where there is more overlap between signal and 
background. To include this effect, we assumed a linear variation of the background shape, from 
— 10% to +10% between the minimum and the maximum of the transverse mass distribution. 
Another source of systematic uncertainty is the rate of the backgrounds. It is expected that the 
background rate (W^t and ti) could be known to 5% [194]. Therefore, to take this effect into 
account, we increase the background rate by 5% while at the same time we decrease the signal by 
5%. Finally, we also include the scale uncertainty: 1% for jets and 0.1% for photons, electrons 
and muons. The overall precisions on the charged Higgs boson mass determination, including the 
systematic uncertainties, are shown in Table [I] and Figure ||. 

2.3 H ± Mass Determination in H ± -> tb 

In the tb channel, the full invariant mass can be reconstructed although this channel suffers from the ir- 



reducible ttb background and the signal combinatorial background [186]. The determination of the mass 
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can be done using the likelihood method described above or by fitting the signal and the background. In 
the latter case, one assumes that the background shape and normalization can be determined by fitting 
outside the signal region, thus, the systematic uncertainties include only the scale uncertainty. We as- 
sume a Gaussian shape for the signal and an exponential for the background and fit signal+background 
including the statistical fluctuations and the scale uncertainty. Both methods are in agreement on the 
mass determination. The results are shown in Table |l| and Figure [|. 

Table 1: The overall precisions on the mass determination are better in the tv channel than in the tb channel. This is due to the 
fact that the latter suffers from large tib and signal combinatorial backgrounds (£ = 100 fb _1 ). 



m H ± (GeV) 




T^~V T 


H±- 


» tb 




< m > 


5m 


< m > 


5m 


225.9 


225.9 


2.9 


226.9 


1.8 


271.1 


271.0 


3.9 


270.1 


10.1 


317.8 


319.7 


5.9 


320.2 


11.3 


365.4 


364.9 


8.1 


365.4 


12.1 


413.5 


414.8 


8.0 


417.4 


17.6 


462.1 


460.7 


10.6 


465.9 


24.1 


510.9 


511.4 


15.7 







2.4 Determination of tan (3 

tan/3 can be obtained by measuring the signal rate in the t^v t channel where the backgrounds are 
relatively low. The main systematic error would come from the knowledge of the luminosity, whose 
uncertainty is taken conservatively to be 10%. The error in the rate measurement can be estimated 
as [fT95j] 



A(<r x BR) _ j S + B | (AC 



ax BR 



S 2 



(5) 



where S and B are the numbers of signal and background events respectively. The uncertainty on tan j3 
is computed as 



Atan/3 ~ A(a x BR) 



d(a x BR) 



(6) 



dtan (3 

The production cross-section for gb — ► tH^ and the branching ratio of — ► t^v t above the top-quark 
mass can be written respectively as |Q] 

cr(gb — > tH ) oc m\ cot 2 (3 + m\ tan 2 j3 

and 

m 2 tan 2 (3 



(V) 



BR(H ± — * t v ) ~ 

T 3(m 2 cot 2 (3 + m 2 tan 2 (3) + m 2 . tan 2 (3 ' 

Using the relations ([7]) and (||), the rate in the t^v t channel at large tan (3 is obtained as: 



a x BR oc tan 2 (3. 



From the Equations @ and @, we get 



Atan/3 _ f A(a x BR) 
tan/3 



2 ax BR 



(8) 



(9) 



(10) 



The expected uncertainties on tan (3 determination from the measurement of the rate in the H ± — r / > 
channel are shown in Table and Figure HI 
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Table 2: The overall precisions on tan/3 determination in the H ± — > r v T channel for C = 30, 100 and 300 fb 1 , and for 
m H ± = 250 GeV. 



tan/3 Atan/3/tan/3(%) 





30 ft" 1 


100 fb -1 


300 ft 


20 


15.4 


10.6 


7.4 


25 


12.2 


8.7 


6.5 


30 


10.5 


7.7 


6.1 


35 


9.1 


7.0 


5.7 


40 


8.4 


6.6 


5.6 


45 


7.7 


6.6 


5.5 


50 


7.3 


6.1 


5.4 
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m H+ (GeV) tanfi 



Fig. 4: The expected overall precision on the charged Higgs boson mass and on tan j3 measurements, as a function 
of the charged Higgs boson mass (left plot) and tan (right plot) respectively. For the mass determination, the 
H ± — > r ± ^ r channel gives better precisions than H ± — > tb except at low m H ± . In addition, — ► i^tv allows 
for the determination of tan /3 by measuring the rate in this channel. 



2.5 Conclusions 

In the r^Ur channel, there is no resonance peak, only the transverse mass is reconstructed. A likelihood 
method is used to estimate the expected precisions on the mass measurements. The systematic effects 
include the background shape, the background rate and the energy scale. The overall relative precision 
in this channel ranges from 1.3% at m H ± = 226 GeV to 3.1% at m H ± = 511 GeV for an integrated 
luminosity of 100 ft -1 . At 300 ft -1 , the precision improves to 0.8% at m H ± = 226 GeV and 1.8% at 
m H± = 511 GeV. 

The tb channel offers a resonance peak with a large background from tib and the signal combinato- 
rial. It is possible to use the likelihood method for the mass determination in this channel. Alternatively, 
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a fit of the signal and background can be performed provided the background shape and normalization 
can be determined by fitting outside the signal region. Results from both methods are in agreement. The 
relative precision in this channel ranges from 0.8% at m H ± = 226 GeV to 5.2% atm H ± = 462 GeV for 
100 fb" 1 . For 300 ftr 1 , the precision improves to: 0.5% at 226 GeV and 3.5% at 462 GeV. 

In either channel, the overall uncertainties are dominated by the statistical errors. The tv channel 
offers better precisions on the charged Higgs boson mass determination than the tb channel, except at 
low m H ± where the r^v T channel suffers from a much reduced cut efficiency. 

tan (3 can be measured in the — > t^v t channel (by measuring the rate) where the background 
is relatively low and the discovery reach is extended to high masses compared to — * tb. Assuming a 
10% uncertainty on the luminosity, the relative precision on tan f3 ranges from 15.4 to 7.3% for tan = 
20 to 50, at low luminosity. For an integrated luminosity of 300 fb -1 , the precision improves to: 7.4% at 
tan /3 = 20 to 5.4% at tan (3 = 50. 



3. Boson in the Threshold Region 

In this section, we discuss charged Higgs boson analyzes in the threshold region, i.e., for m H ± ~ m t , 
taking into account the correct description of the charged Higgs boson production and decay mechanism 
in this region. 



3.1 Motivation 

In the MSSM, the LEP 2 limits on the mass of the lightest Higgs boson convert directly (see [§]) 
into a lower bound on the charged Higgs boson mass (at least, at low tan/3, say, around 3): m 2 H± rj 
m^ ± + vr? h Q > (140 GeV) 2 [196]. Whereas the charged Higgs boson mass region just above this value 



is theoretically well understood, the description of the so-called (top) 'threshold region', m H ± ~ m t , 
requires careful considerations when it comes to production and decay in the context of a Monte 
Carlo (MC) simulation, as explained below. The main production mode of scalars with mass strictly 
below the top-quark mass, m H ± < m t , is the decay of the top (anti)quarks themselves, the latter be- 
ing produced via QCD in the annihilation of gluon-gluon and quark-antiquark pairs. So far, standard 
MC programs, such as PYTHIA, HERWIG and ISAJET JT00][T8T1^82| ], have accounted for this pro- 



cess through the usual procedure of factorizing the production mode, gg, qq — > tt, times the decay one, 
i — > bH~, in the so-called Narrow Width Approximation (NWA). However, this description fails to cor- 
rectly account for the production and decay phenomenology of charged Higgs bosons when their mass 
approaches or exceeds that of the top-quark, hence undermining the ability of experimental analyzes in 
pinning down the real nature of these particles (if not detecting them altogether). 



This is particularly a pressing issue at the Tevatron Run 2 [ |188| ], as the collider reach in m H ± 



dips precisely into the threshold region [198]. Here, the use of the 2^3 hard scattering process 



gg, qq — > tbH~ [|199J, in place of the 'factorization' procedure Q200[ ], is mandatory, as one can clearly 
see from Figure |5| where the discrepancies in the shape and normalization come from the Higgs-strahlung 
mechanism and the relative interferences as mentioned earlier. 

Also differential distributions can strongly be affected by an approximated modeling of the produc- 
tion and decay process in the threshold region, as one can appreciate from Figure ||. Here, differences 
are clearly sizeable also for the top quark. However, in this case one should expect the impact to be 
marginal, as this particle is actually unstable and since its three-body decay products are subject to the 
effect of usual detector resolution uncertainties. In contrast, this is no longer true for the bottom quark, 
which fragments directly into hadrons. Besides, the availability of the newly implanted silicon vertex 
detector may render the tagging of 5-quarks a crucial ingredient in detection strategies of charged Higgs 



bosons at Run 2, pretty much along the same lines as established at the LHC [189]. 
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Fig. 5: Cross section for gg, qq — ■> tbH~ and gg, qq — » it — » tfeiJ - in NWA, at the Tevatron with y's = 2 TeV, as 
a function of m H ± for a representative value of tan j3 (the kinematical effects discussed are the same irrespective 
of the latter). Hereafter, charge conjugated rates are always included. Besides, both top and bottom quark masses 



in the Higgs boson Yukawa couplings are non-running and set to 175 and 4.25 GeV, respectively. CTEQ4M 1 197] 
is used for the Parton Distribution Functions (PDFs), with scale m H ± . 




Fig. 6: Transverse momentum distributions of the final state quarks in gg, qq — * tbH~ (solid) and gg, qq — > tt — > 
tbH~ (dashed) in NWA, at Tevatron with y/s — 2 TeV, for m H ± =170 GeV. Again, the actual tan/3 value is 
irrelevant. 



3.2 Analysis 

If one looks at the most promising (and cleanest) charged Higgs boson decay channel, i.e., 



i v T [ [20 1| ], while reconstructing the accompanying top quark hadronically, the prospects of detection 
are rather good. This is made clear in Table ||. Even if one neglects the tagging of the fo-quarks in the 
final state, the final results are very different between the full process and the NWA. In Table |, we have 
reported the signal and dominant (irreducible) background rates (that is, from gg, qq — ► tbW~ + c.c. 
events, yielding the same final state as the signal) after the following sequence of cuts: 

1. Tau-jets are selected if they satisfy the criteria: pTp >15 GeV and \r] T \ < 2.5. 

2. We require j> T > 20 GeV, since the presence of neutrinos from H~ decays and invisible decay 
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products of r's (mainly i s) implies that a significant fraction of the transverse momentum goes 
undetected. 

3. Quark-jets are selected by imposing pip > 20 GeV and \rf\ < 2.5. We require at least one of these 
to be tagged as a 6-jet. 

4. We demand that two un-tagged jets have an invariant mass around m w ±, e.g., \m q g> — m w ±\ < 
10 GeV and that the 6-jet in combination with the other two un-tagged jets produces an invariant 
mass close to m t , e.g., |mfw — m t \ < 15 GeV. 

5. We require that the reconstructed transverse mass, Equation (Q), be above the IF 1 * 1 -boson mass: 
niT > rn w ± « 80 GeV. 

The r's can be tagged as narrow jets in their 'one-prong' hadronic decay modes — see the rela- 
tions @ — which represent 90% of the hadronic decay rate and about 50% of the total. This distinguish- 
ing feature is in contrast to the typical appearance of quark- and gluon-jets, which yield 'multi-prong' 
hadronic topologies in the detectors, typical of QCD backgrounds of the form + jets and Z° + jets. 

Table 3: The signal rates (in fb) for the process qq,gg — » tbH~(~^ t~v t ), at Tevatton with ^/s = 2 TeV, 
for representative values of m H ± and tan j3, after all cuts described in the text. The corresponding rate of the 
background is 0.22 fb independent of m H ± . 



m H ± (GeV) i /tan/3 


3 


6 


40 


150 


6 


3 


52 


160 


2.8 


1.5 


22 


170 


.4 


0.25 


3.5 


175 


.13 


.08 


1.42 


180 


.067 


.061 


1.09 



3.3 Conclusions 

In the end, despite the fact that one should more realistically expect both signal and background rates 
to be further reduced by a factor of 4 or so (r-identification efficiencies are estimated to be of order 



50% [202], similarly for the tagging of any 6-jet [198]), the final message that emerges is that the chances 



of extracting the — > t^v t signal after 15 fb -1 of luminosity at the Tevatron Run 2 are rather good 
for m H ± up to 1 80 GeV or so at large tan j3, while being negligible at low to intermediate tan j3 values. 
Conclusions would obviously be drastically different in the NWA scenario, if one recalls Figure |5[ 

The situation can be improved even further by taking advantage of the r polarization effects as 
explained above. In this respect, the enforcement of the constraint @ reduces the background by a 
factor of 5, while costing to the signal only a modest — in comparison — 50% suppression (for any 
charged Higgs boson mass in the usual interval between 160 and 190 GeV). 

Although we have relied here on a parton-level analysis, it is clear that its main features would 
remain valid even in presence of fragmentation/hadronization effects. In fact, work is currently ongoing 
in order to include the latter, as well as a more realistic detector simulation, to emulate the real potential 
of the Tevatron experiments, by exploiting the mentioned 2^3 description of the production 
dynamics and the spin correlations in r-decays, as they are now both available in version 6.4 []203|] of the 



HERWIG event generator (the latter also through an interface to TAUOLA Q45[|). 

The problematic is very similar at the LHC, if anything more complicated. In fact, at the CERN 
hadron collider, the above 2^3 reaction is dominated by the gg-initiated subprocesses, rather than by 
qq annihilation, as is the case at the Tevatron. This means that a potential problem of double counting 
arises in the simulation of tbH~ + c.c. events at the LHC, if one considers that Higgs-strahlung can also 
be emulated through the 2 — > 2 process bg — * tH~ + c.c. The difference between the two descriptions is 
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Fig. 7: Cross section for gg, qq — > tbH~; gg, qq — ► tt — * i&il - with finite top-quark width; &g — * and the 
combination of the first and the last, at the LHC with ^/s = 14 TeV, as a function of m H ± for a representative 
value of tan (3. 



well understood, and a prescription exists for combining the two, through the subtraction of a common 
logarithmic term: see Refs [204-207]. Figure [7| summarizes all the discussed issues in the context of 



the LHC. The 2^3 process is available in HERWIG and detailed simulations of the r v T channel at 
the CERN hadron collider are now possible also for the threshold region, as already done for other mass 



intervals [185-187, 19C, 193, 206] 



4. Boson in Large Extra Dimensions 

In this section, we discuss the LHC sensitivity to the charged Higgs boson discovery in the channel 
H~ — > t~^v in models with singlet neutrinos in large extra dimensions. The observation of such a signal 
would provide a distinctive evidence for these models since in the standard two Higgs doublet model 
type II, H~ 



t l v is completely suppressed. Details of this analysis can be found in [209] 



4.1 Motivation 

In models where extra dimensions open up at the TeV scale, small neutrino masses can be generated 
without implementing the seesaw mechanism [21C]. These models postulate the existence of 5 addi- 



tional spatial dimensions of size R where gravity and perhaps other fields freely propagate while the 
SM degrees of freedom are confined to (3+l)-dimensional wall (4D) of the higher dimensional space. 
The true scale of gravity, or fundamental Planck scale M#, of the (4 + <5)D space time is related to the 
reduced 4D Planck scale M PU by M pl = R 6 M% +2 , where M P{ = 2.4 x 10 18 GeV is related to the usual 
Planck mass 1.2 x 10 19 GeV = s/SttMpi. Since no experimental deviations from Newtonian gravity are 



observed at distances above 0.2 mm Q21 lQ , the extra dimensions must be at the sub-millimeter level with 
M* as low as few TeV and 5 > 2. 

The right handed neutrino can be interpreted as a singlet with no quantum numbers to constrain it 



to the SM brane and thus, it can propagate into the extra dimensions just like gravity [|86[]. Such singlet 
states in the bulk couple to the SM states on the brane as right handed neutrinos with small couplings 
- the Yukawa couplings of the bulk fields are suppressed by the volume of the extra dimensions. The 
interactions between the bulk neutrino and the wall fields generate Dirac mass m o terms between the 
wall fields and all the Kaluza-Klein modes of the bulk neutrino: 



rriD 



A M« 
y/2Mpi 



(11) 
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where A is a dimensionless constant. The mixing between the lightest neutrino with mass m o and the 
heavier neutrinos introduces a correction N to the Dirac mass such that the physical neutrino mass m v 
is m v = tud/N, where 



n is a vector with 5 integer components counting the number of states and the summation is taken over 
the Kaluza- Klein states up to the fundamental scale M*. The sum over the different Kaluza- Klein states 
can be approximately replaced by a continuous integration. As shown in Table f|, small neutrino masses, 



m u , can be obtained consistent with atmospheric neutrino oscillations [212]. The spectrum of many 



Table 4: The parameters used in the current analysis of the signal with the corresponding polarization asymmetry. 
In general, H~ would decay to t~[ and t^, H~ — > t^P + r£ip, depending on the asymmetry. For the decay 
H~ — > TfiD (as in MSSM), the asymmetry is — 1. The signal to be studied is H ~~ — > r£tp. 







M* (TeV) 


Su,5 


m H ± (GeV) 


tan/3 


A LR 


m u (eV) 


Sig 


■-1 


2 


4,4 


219.9 


30 


~ 1 


0.5 10" 3 


Sig 


.-2 


20 


3,3 


365.4 


45 


~ 1 


0.05 


Sig 


.-3 


1 


5,6 


506.2 


4 


~ 1 


0.05 


Sig 


.-4 


100 


6,6 


250.2 


35 


1 


0.005 


Sig 


.-5 


10 


4,5 


350.0 


20 


1 


0.04 


Sig 


.-6 


50 


5,5 


450.0 


25 


1 


0.04 



extensions of the SM includes a charged Higgs boson state. We consider as a prototype of these models 
the 2-Higgs Doublet Model of type II (2HDM-II). H~ decays to the right handed t~ through the r 
Yukawa coupling: H~ — ► t^v. The H~ decay to left handed r~ is completely suppressed in MSSM. 
However, in the scenario of singlet neutrinos in large extra dimensions, H~ can decay to both right 
handed and left handed r~ depending on the parameters M#, mp, 5, m H ± and tan j3, due to the large 
number of Kaluza-Klein states of the right handed bulk neutrino: H~ — > P + r^if), where ip is a bulk 
neutrino and z/ is dominantly a light neutrino with a small admixture of the Kaluza-Klein modes of the 
order mrji?/|ra|. The measurement of the polarization asymmetry 

_ r(ff-^T-yo-r(g-^T-p) 

can be used to distinguish between the ordinary 2HDM-II and the scenario of singlet neutrinos in large 
extra dimensions. 

The singlet neutrino may propagate into a subset 5 V (5 V < 5) of the S additional spatial dimensions, 
in which case the formalism for the generation of small Dirac neutrino masses is merely a generalization 



of the case b v = 5 [213] 



The charged Higgs boson decay to right handed t, H — > t r v has been extensively studied for the 



LHC Q190| , |193| ]. Here we discuss the possibility to observe H~ — * r^ip at the LHC above the top-quark 
mass. Table |] shows the parameters selected for the current analysis. The cases where the asymmetry is 
+1 are discussed in details. We assume a heavy SUSY spectrum with maximal stop mixing. The present 



analysis is conducted in the framework of PYTHIA6.1 and ATLFAST [10C, 138], and the Higgs boson 



masses and couplings are calculated to 1-loop with FeynHiggsFast [ 214 ]. 
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4.2 Analysis 

We consider the 2 — ► 2 production process where the charged Higgs boson is produced with a top- 
quark, gb — > tH^ as shown in Figure ||. Further, we require the hadronic decay of the top-quark, 
t — ► — > jj& and the charged Higgs boson decay to T-leptons. The major backgrounds are the single 



Table 5: The expected rates (crx BR), for the signal gb — > tH with H~ — > t^D + r£ip and t — > jjfe, and for the 
backgrounds: ITt and ti with IT - — > t~[v and — > jj. We assume an inclusive ti production cross section of 
590 pb. Other cross sections are taken from PYTHIA 6.1 with CTEQ5L parton distribution function. See Table Q 
for the parameters used for Sig.-l, Sig.-2 and Sig.-3. In the last columns, we compare the H~ — > t^D branching 



Process 


a x BR (pb) 


BR 


BR(MSSM) 


Sig.-l 


1.56 


0.73 


0.37 


Sig.-2 


0.15 


1.0 


0.15 


Sig.-3 


0.04 


1.0 


0.01 


tt 

gb ^Wt(p T > 30 GeV) 


84.11 
47.56 







top production gb — > Wt, and ti production with one W + — ► jj and the other W~ —>■ Depending 
on the polarization asymmetry, H~ — > t£ v will contribute as an additional background. In Table we 
list the rates for the signal and for the backgrounds. The polarization of the r ± -lepton is included in this 
analysis through TAUOLA [45]. We consider the hadronic one-prong decays of the r ± -lepton — see the 
relations (0) — which are believed to carry a better imprint of the r ± polarization [191]. 

For the signal in MSSM, right handed t^'s come from the charged Higgs boson decay, H~ — * 
t^v, while in the backgrounds, left handed r^'s come from the decay of the W~(— * t£u). Because 
of the neutrino in the final state, only the transverse mass, Equation @, can be reconstructed. In the 
framework of large extra dimensions, we are interested in H~ — > T£tp where the polarization of the 
r-lepton would be identical to the background case but opposite to the MSSM case. Therefore, the 
polarization of the T-lepton would not help in suppressing the backgrounds. Nevertheless, there are 
still some differences in the kinematics of the signal and of the backgrounds: the different transverse 
mass bounds and the increasingly harder r-jet and missing transverse momenta as the Higgs boson mass 
increases. The reconstructed transverse mass for the signal and the backgrounds are shown in Figure || 



Table 6: The expected signal-to-background ratios and significances calculated for an integrated luminosity of 
100 fb _1 (one experiment). See Table ^ for the parameters used for Sig.-l, Sig.-2 and Sig.-3. In all the cases 
considered, the signal can be observed at the LHC with significances in excess of 5-a at high luminosity. 



Sig.-l Sig.-2 Sig.-3 



Signal events 


41 


215 


16 


ti 


7 


7 


7 


Wt 


3 


3 


3 


Total background 


10 


10 


10 


S/B 


4.1 


21.5 


1.6 


s/Vb 


13.0 


68.0 


5.1 



and the expected signal-to-background ratios and the signal significances in Table ||. The reconstruction 

17 There is no enhancement in the background rate from the contribution W~ — > t£i/j. 
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Fig. 8: The reconstructions of the transverse mass of the signal in MSSM, the signal in models with a singlet neu- 
trino in large extra dimensions and of the backgrounds, for an integrated luminosity of 100 fb _1 . The observation 
of the signal in the transverse mass distribution would not be sufficient to identify the model: the r ± polarization 
effects must be explored further. 



of the transverse mass is not enough to distinguish between the MSSM and the singlet neutrinos in large 
extra dimensions. The differences in these two scenarios are best seen in the distribution of p* j ' E T ~i et , 
the fraction of the energy carried by the charged track which is shown in Figure ||. In the MSSM, this 
distribution peaks near and 1 while in H~ —> r£ip from large extra dimensions and in the backgrounds, 



this distribution peaks in the center. The backgrounds are relatively very small, and as concluded in Q190 



193], the discovery reach is limited by the signal size itself. Therefore the observation of a signal in the 



transverse mass distribution and in the distribution of the fraction of the energy carried by the charged 
track should help determine whether the scenario is MSSM or not. 



4.3 Conclusions 

We studied the observability of the channel H~ — ► t£u in models with a singlet neutrino in large extra 
dimensions at the LHC. Although the observation of a signal in the transverse mass distribution can be 
used to claim discovery of the charged Higgs boson, it is insufficient to pin down the scenario that is 
realized. Additionally, by reconstructing the fraction of the energy carried by the charged track in the 
one-prong r ± decay, it is possible to claim whether the scenario is the ordinary 2HDM or not. The 
further measurement of the polarization asymmetry might provide a distinctive evidence for models with 
singlet neutrinos in large extra dimensions. 
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M, = 2 Iff Ge V, 5 = 3, m D =145 eV, m H _ = 365.4 GeV, fanp" = 45 




p IE 



Fig. 9: The distribution of the ratio of the charged pion track momentum in one prong r decay to the r-jet energy 
for rriA — 350 GeV, tan/3 = 45, A/* = 20 TeV, S = 3 and m„ = 0.05 eV, and for an integrated luminosity 
of 100 fb _1 . In the 2HDM-II, this ratio would peak near and 1 as shown while in other models, the actual 
distribution of this ratio would depend on the polarization asymmetry since both left and right handed r's would 
contribute. In the case shown, the asymmetry is ~ 1 and the ratio peaks near the center of the distribution. 



5. H ± Decays into SUSY Particles 

Thus far, all analyzes have implicitly assumed that the SUSY counterparts of ordinary particles had 
masses much higher than m#±. However, lowering the typical SUSY mass scale may induce new inter- 
actions among H bosons and several of the sparticles, so that the former may be abundantly produced 



in the decay of the latter (e.g., from gluinos and squarks, see Q223Q ) or, alternatively, new Higgs boson 
decay channels into SUSY particles may well open at a profitable rate [215]. We will defer the study of 
the first scenario to section ^J. Here, we will investigate in some depth the second possibility, focusing 
on the 'intermediate' tan [3 region (say, between 3 and 10) left uncovered by the SM decay channels (see 
Figure HQ 



5.1 Motivation 



It was demonstrated in [ [215p that the decays of a charged Higgs boson into a chargino and a neutralino 
could probe regions of the MSSM parameter space where decays to SM particles, such as bi or r~v T 
yield no significant signal — see Figure ||. In particular, intermediate values of tan (3 between ~ 3 and 
~ 10 were in part accessible via — » xfX{2 3} modes resulting in three lepton final states (where 
leptons mean electrons or muons), a hadronically reconstructed top quark (from gg — * btH~ , gb — > tH~ 

18 Depending on the actual rate of the new SUSY decay channels in the complementary tan/3 areas (less than three or larger 
than ten), some rescaling to the discovery reaches via the SM modes of Figure ^ may be needed. This will be addressed in 



Ref. [227] 
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and their charge conjugate production processes) plus substantial missing transverse momentum (from 
neutralino and chargino decays to the stable lightest neutralino, i.e., the lightest supersymmetric particle 
or LSP). 

We refer to the charginos and neutralinos collectively as "inos", which are the mass eigenstate 
mixtures of the electroweak (EW) gauginos and higgsinos. We expand on the parton level analysis 
of [ 215| |: firstly, by studying the signal in a full event generator environment with an improved back- 



ground analysis that includes potential MSSM background processes (the previous study considered 
only SM backgrounds); secondly, by further investigating the possible role of on-shell or off-shell slep- 
tons (the supersymmetric partners of the leptons) in the signals. As noted in the previous study, if there is 
a light slepton, then the leptonic branching ratios (BRs) of the inos can be significantly enhanced (espe- 
cially those of x§> Since both inos and sleptons play key roles in the signal process, a considerable 
number of MSSM input parameters are relevant. We seek to scan this expanse of parameter space, at least 
at the parton level, identifying parameter points or series of points that merit dedicated event generator 
analyzes. 

To correctly delineate the portions of parameter space that can potentially yield viable signals, 
it is necessary first to know where experimental constraints cut off otherwise favorable zones. LEP 2 
experiments have yielded a number of bounds on MSSM particles and parameters that impinge directly 
on our preferred signal regions. Among these, the most crucial are []216|]: the mass limit on the lighter 



chargino, taken as m~± > 103 GeV (which allows for the possibility of a light sneutrino resulting in 
negative interference from the then-significant i-channel in chargino pair production); the mass limits 
on the sleptons, taken as mg 1 > 99.0 GeV, > 91.0 GeV, m fl > 85.0 GeV and mc > 43.7 GeV 
(the last being from studies at the Z° pole); and the regions excluded by searches for signals of the type 
e + e~ — > Z°*h and e + e~ — ► Ah. For the numerical limits just given, it is assumed that m~± — m~o, 

Xi Xi 

m-n — rrir-o > 5 GeV. The constraint from Higgs boson production is somewhat vague owing to an 



estimated uncertainty in the expected mass of the light Higgs boson of 2-3 GeV from un-calculated 



higher order corrections [217] and up to 5 GeV from the error in the measurement of the top quark 



mass, m t = 174.3 ± 3.2 ± 4.0 GeV Q2 1 8Q . A small shift in the light Higgs boson mass translates into a 
substantial shift in the location of the bound seen in the tan f3 versus tua plane. Other LEP 2 limits which 
generally should be less restrictive than those just mentioned are also incorporated into our analysis. 

There are also other processes where charged Higgs bosons (or A, to whose mass that of the is 
closely tied) enter as virtual particles at the one-loop level. These include neutral meson mixing (K°K°, 



D°D°, or B°B°) pE| , g20| |, Z° -> bb (R b ) [|220|,[22l]], and b -► 57 decays p!9| - ^ |, the last of which, 



where restrictions on m H ± are linked to a number of MSSM variables, notably including the masses of 



the lighter chargino and the stops, is generally thought to be the most constraining [221 ]. This constraint 



may well disqualify regions of the parameter space where our signal is strong and otherwise allowed; 



however, the applicability of these low energy constraints is unclear due to a variety of factors (see [ [2 1 5p ) 
and thus they will not be included in this analysis (though we do choose our stop parameters with an eye 
towards attempting to evade potential b — ► 57 bounds). 

5.2 Parameter Space Exploration 

We begin in Figure [Rj (top plot) with a look at the raw cross-section for gg — ► btH~ , gb — > tH~ and 
their charge conjugate production processes in the tan (5 vs itia plane. Other MSSM parameters are fixed 
at the values noted in the figure caption. Here we correctly take into account the subtraction needed to 



avoid double counting between the 2^2 and the 2^3 production processes [204, 205]. The charged 



Higgs boson mass is calculated including radiative corrections as contained in ISAJET Q182[ ] and the 



CTEQ4L [197] structure function set is employed. 

In Figure [R] (bottom plot), we fold in the BRs of a charged Higgs boson into an ino pair (restricted 
here to either xfx® or xfxis) multiplied by the BR of the ino pair into a trio of charged leptons l+fl'^ 1 



100 



(recall that "£" herein denotes either an electron or a muon), where £ and £' may or may not be of the 
same flavor. In calculating the leptonic BRs of the inos all conceivable decay chains are taken into 
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Fig. 10: (top plot) Raw cross-section (in fb) for gg — > btH~ , g& — > and their charge conjugate production 
processes in the tan (3 versus rriA plane. The proper subtraction factor to avoid double counting between the 2^2 
and the 2^3 production processes is included, (bottom plot) Expected number of pp — > XtbH^ , XtH^ — > 
£ + £~£ l:iz + t events per 100 fb -1 at the LHC (with no cuts), assuming the charged Higgs boson decays into X1X2 



or X\ X*3- Here M 2 = 200 GeV and p = —115 GeV. One-loop formulae as found in [182, 228] are used to relate 
m H ± to v/ia- TUt = 175 GeV and m& = 4.25 GeV. Other MSSM input parameters are: rriq = 1 TeV for the first 
two generations, m f = 600 GeV, m f = 500 GeV, mi = 800 GeV, A t — 500 GeV, A b — 0; mz = 150 GeV, 
m i = 771 — 200 GeV and Ai = for all three generations. 



account^. The plot shows that, assuming an integrated luminosity of 100 fb -1 , hundreds to thousands 
of events are expected for ttla < 400 GeV. The possibility of extracting the SUSY signal exists across 
the full range of allowed tan/3 values and for (and m H ±) ~ 400 GeV provided the threshold for the 
— > xf X{2 3} decay is exceeded. Higher values of tan (3 and ~ 240 GeV are optimal choices 
(quite different than when only the raw production rate is considered). The region from the top of the 
plot down to the dashed curve marked as ' 103 GeV is excluded by the LEP 2 limit on the chargino mass, 

"including possible extra minor contributions arising when a neutralino decays into a \i ( or a W ± ) and one lepton and the 
xf then decays into Xi along with the second lepton (or W ± decays leptonic ally). These decay modes yield extra neutrinos 
in addition to the final products of the main decay modes, but should be experimentally indistinguishable. Possible modes with 
x'i — * X2X are also taken into account. 
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and the region below the dotted curve is excluded by Higgs boson production. As shown on the plot, the 
upper bound is fairly sensitive to the m~± limit; and thus in turn very dependent on the values chosen 

for other MSSM parameters, in particular the higgsino mixing parameter fi and M 2 ^|, which are chosen 
to be favorable to our signal in the plot. Also, as mentioned in the previous section, the exact location of 
the latter bound is fairly loose due to uncertainty in the mass of the light Higgs boson. 

The expected number of events is shown again in Figure 111], this time in the M 2 versus fi plane 
with tan/3 fixed at 8 and rriA = 290 GeV {m H ± ~ 300 GeV). The shaded region is excluded by the 
LEP 2 bound on the chargino mass. Again we see that hundreds to thousands of events are possible in 
un-excluded regions of the parameter space; but it is also apparent that small values of are strongly 
preferred. This is a serious restriction which means that the preferred and un-excluded signal region is 
just beyond that region probed by LEP 2 (and thus also the region expected to be searched in the first 
phase of a future e + e~ linear collider). 




Fig. 11: Expected number of pp XtbH ± ,XtH ± -> ± + t events per 100 fb 1 at the LHC (with no 

cuts), seen in the M 2 versus /i plane. Here tan (3 = 8 and ttia = 290 GeV (m H ± ~ 300 GeV). Other inputs as in 
Figure [To] 



Figure 12 shows contours for the raw number of three-lepton events (11 l + + p™ lss + t) 
expected at the LHC for <fs = 14TeV and C = 100 fbr 1 in the M 2 versus mz plane before any 

selection cuts are applied. The soft slepton mass spectrum is fixed by = m~ = m| + 50 GeV 
(soft slepton masses are assumed to be degenerate for all three generations; the Ae's are still kept equal 
to zero) and \i is set to —115 GeV. All other inputs are as in the previous figures. The shaded region 
is excluded by LEP 2 limits on m~± (below) and ra i (left)0. We see that over a thousand events are 
possible. The optimal spot (excluding the sliver region) is at (m^, M 2 ) ~ (110.0 GeV, 195.5 GeV) 
which boasts well over a thousand events. 



20 Mi and M 2 are the U(1) Y and SU(2) L gaugino masses, respectively; Grand Unified Theories (GUTs) predict gaugino 
unification and M t = | t&n 2 8 w M 2 ~ 0.5032 M 2 — |M 2 , as will be assumed in all numerical calculations. 

21 There is a small "sliver" of allowed parameter space that slices through the excluded region to the left. This is where the 
charged sleptons are only slightly more massive than the LSP (taken as within 5 GeV) and thus can evade the LEP 2 searches. 
Note though that then we expect to get soft leptons from the slepton decay and so our signal probably disappears here. 
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5.3 Analysis 

Here, we will show how all SM backgrounds can be completely removed, leaving only MSSM processes 
as irreducible backgrounds in the 3£ + p™ lss + 1 channel. The relevant SM noise is constituted by WZ°, 
Z°Z°, bbZ°, tt, tt£ + £~ (consisting of both ttZ° and tij*, but not their interference), tbW and tbW£ + £' 
(again, via Z° and 7*) production and decay. As for the MSSM backgrounds, one has to deal with pair 
production of sleptons, gauginos or squarks/gluinos and with neutral MSSM Higgs bosons produced in 
association with heavy quark pairs (SM-like tth production is found to merit attention while associated 
production of the heavier neutral Higgs bosons is negligible). 

We simulated the inclusive — > 3/ signal and the MSSM backgrounds at the 5 points in the 
intermediate tan (3 region of the MSSM parameter space listed in Table [7]. These points were chosen in 
the favorable regions obtained in the previous section. In each case we have also chosen m g = 700 GeV, 
ma = 1000 GeV, m- h = 800 GeV, m, = 600 GeV, m f = 500 GeV and A t = 500 GeV. Rather 

1 b R ' t L ' t R <■ 

large gluino and squark masses are chosen to preclude large charged Higgs boson production rate from 
MSSM cascade decays [223], thus leaving the 'direct' production modes of the previous section as the 
only numerical relevant contributors at the LHC [224]. Sleptons are chosen to be light in accordance 
with the discussion in the last section. 



Table 7: Simulated MSSM parameter points. All masses in GeV. The event number is the parton-level result 

H ± -> xfx\ 2 3} -> M x P er lOOfbT 1 . M 1 = \M 2 is assumed. 
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HERWIG 6.3 [181[ is used to generate all hard processes^, fragmentation and showering, with 
adoption of the default settings JT8l|]. The SUSY spectrum was obtained from ISASUSY 7.58 [pi 



"An exception is the ttj* process, which uses a set of separate subroutines [225]. 
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through the ISAWIG interface [ j226[ ]. The detector aspects were simulated using CMSJET 4.801 $3^ , 
which contains a parameterized description of the CMS detector response; however, since none of our 
selection cuts depend strongly on unique performance factors of the CMS detector, we expect the analysis 
here to roughly coincide with what the ATLAS detector environment would yield. The effects of pile-up 
at high luminosity running of the LHC have not been included. 

In order to distinguish between signal and background events, the following selection criteria are 
applied: 

1. We require exactly three isolated leptons (£ = e, fi), with px > 20, 7 and 7GeV, respectively, 

all within \rj\ < 2.4. The isolation criterion demands that there are no charged particles with 
pr > 1-5 GeV in a cone of radius R = 0.3 radians centered around each lepton track and that the 
sum of the transverse energy in the ECAL (Electromagnetic Calorimeter of CMS) crystal towers 
between R = 0.05 and R = 0.3 radians is smaller than 3 GeV. 

2. We impose a Z°-veto, i.e., rejection of all di-lepton pairs with opposite charge and same flavor that 

have an invariant mass in the range mz ± 10 GeV. 

3. We select only events with three or four jets, all with pt > 20GeVand|r/| < 4.5. 

4. Among these jets, we look for the combination of three that are most probably coming from a top 

decay, by minimizing the difference rrijjj — m t , rrijjj being the invariant mass of the three-jet 
system. These three jets must have rrijjj in the range m t ± 35 GeV. 

5. One of the jets must be 6-tagged (the signed significance of the transverse impact parameter must 

be a(ip) = > 2.5). The two other jets must have their invariant mass, Mjj, in the range 

m w ± 20GeV^ 

6. The px of the hardest lepton must be less than 80 GeV^ 

7. The transverse missing momentum is in the following interval: 50 GeV < j> T < 200 GeV. 

8. The p T of the hardest jet is < 180 GeV. 

9. We exploit the 'effective' transverse mass of Ref. [ ]215| ], by imposing 



M eff = J2p^ T (l - cos A(f>) < 140 GeV 
(where pip is the three lepton transverse momentum and A<p the azimuthal angle between pip and 
Pi)- 

10. In addition, the three-lepton invariant mass, m^i, must be < 120 GeV. 

In Table || we summarize the signal and background events that remain after applying these criteria, 
assuming a luminosity of 100 fb _1 and optimizing the selection using the itla = 250 GeV signal. Due to 
the aforementioned small mass differences in the ino sector, leptons coming from the signal events will 
often be soft, which explains much of the difference between the numbers in the "3£ events" column of 
Table || and those from the parton-level analysis in Table [7|. Clearly, it is vital to have a low pt threshold 
for accepting leptons. Conservatively, here we have set this to 7 GeV for both electrons and muons. 

Requiring three tightly isolated leptons in addition to a hadronically reconstructed top quark allows 
us to eliminate most of the SM backgrounds: W^Z , Z°Z°, bbZ° as well as the initially large tt channel. 
The tt£ + £~ background, resulting from ttZ°, ti'j* and tih° production can easily mimic the signature 
of the signal. However, after applying the Z°-veto, the previously dominant tiZ° component of this 
background becomes negligible compared to the intrinsic MSSM backgrounds^ Among the latter, the 
channels that survive the selection are squark/gluino and gaugino pair-production (in which we include 
Xg too) as well as the associated production of a light Higgs boson with a ti pair. In gaugino-gaugino 
production, a gluon jet can fake one of the three jets from the top decay, but in the next stage of the top 
reconstruction these events will always be rejected (see Table |8|). On the other hand, squark/gluino and 

23 This reflects the relatively small mass difference between the parent ino and the daughter ino (the latter is typically the 
LSP). 

24 The single-top counterparts, tbW~ , are sub-leading, as they amount to a ~ 25% correction to the double-top rates. 
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Table 8: Number of events after the cuts mentioned in the text at a luminosity of 100 fb 
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^Includes 6-tagging efficiency for the third jet 



squark-gaugino background events that pass the top reconstruction cut are likely to have a real top in the 
final state, or else a bW + pair. 

The events that pass the top reconstruction cut can be distinguished from the signal by noticing 
that MSSM cascade decays will typically produce harder leptons, more and harder jets, as well as a larger 
amount of missing pr, with respect to the signal. Also, one can further suppress these backgrounds by 
putting an upper limit on the three lepton invariant mass and on the mass made up from the pr of the 
three leptons and the j> r . In particular, these last two variables are very effective if the mass of the Higgs 
boson is just above the chargino+neutralino threshold, i.e., for tha — 250 GeV. For larger charged Higgs 
boson masses, the lepton pr cut-off increases and the discriminative power of this variable is reduced. 



5.4 Conclusions 

The inability to cover the intermediate tan f3 region by exploiting charged Higgs boson decays into SM 
particles prompted us to carefully investigate the potential offered by other channels. In SUSY models 
such as the MSSM, it is natural to explore to this end the interactions between charged Higgs bosons and 
the inos. Assuming that squarks and gluinos are heavy enough so that the production rates of H ± states 
via cascade decays of sparticles are negligible, the obvious place to look is the decay of a charged Higgs 
boson itself into a chargino-neutralino pair. 

The limited XtH~ (and the charge conjugate) production rate precludes exploration for large 
values of m H ± ; however, for m H ± ~ 300 GeV, a signal could well be observed above the background, 
provided that: (1) \i and M 2 are not much above the current LEP restrictions from ino searches; (2) 
sleptons are sufficiently light. We are thus guided to a region of parameter space where both the ino 
and the Higgs boson spectra will be accessible at the LHC, making it a quite reasonable choice for initial 
phenomenological studies. (A more refined analysis encompassing a wider span of the MSSM parameter 
space will be ready in due course [227].) 

Before closing, we would like to remark on possible improvements in our analysis. One such 
improvement would be the inclusion of the leptonic decays of the top, in which case one can look for 
a signature with four leptons and a fo-jet in the final state [228)0. After applying a reasonable set of 

25 In fact, it should be noted that both leptonic and hadronic decays of top quarks (and of inos) have been allowed in the gener- 
ated events of the present analysis, although the former have negligible impact because of the enforced top mass reconstruction 
procedure. 
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cuts on the lepton momenta, we find the SM noise can be completely eliminated. Furthermore, rejecting 
events with more than two jets and limiting the E™ lss eliminate the background from squark/gluino pair 
production. Ino (again, including \g) and slepton pair production as well as associated production of a 
light neutral Higgs boson with a heavy quark pair then constitute the principal backgrounds, though the 
former are largely reduced by requiring one ^-tagged jet. 

Another extension under development is the study of charged Higgs boson decays including the 
higher ino mass eigenstates (xt an d xl) f° r m H ± ~ 300 GeV. For the highest mass inos, the number 
of possible decay chains ending in the LSP can quickly multiply. Such decay modes may have preferred 
regions of parameter space for fi and M 2 quite different from the cases presented here. 

Although limited to some restricted region of parameter space, we have managed with our present 
analysis to cover some portion of the elusive intermediate tan f3 region in the MSSM; and we regard our 
results thus far as encouraging enough to look further into these more exotic charged Higgs boson decay 
modes [227]. 



6. H± from SUSY Cascade Decays at the LHC 

In this section, we analyze the cascade decays of the scalar quarks and gluinos of the MSSM, which are 
abundantly produced at the LHC, into heavier charginos and neutralinos which then decay into the lighter 
ones and charged Higgs particles. We show that these decays can have substantial branching fractions. 
The production rates of these Higgs bosons can be much larger than those from the direct production 
mechanisms, in particular for intermediate values of tan (3. An event generator analysis shows that the 
detection of bosons produced through this mechanism is possible. 



6.1 Motivation 

As previously recalled, for masses m H ± > mt, the two production mechanisms with potentially sizeable 
cross sections at the LHC are the 2 — > 3 and 2 — > 2 processes — shown in Figure [l] — which have to 
be properly combined to avoid double counting [204-207]. However, the cross sections are rather small: 
even for the extreme values, tan (3 = 2 and 40, they hardly reach the level of a picobarn for a charged 
Higgs boson mass m H ± = 200 GeV. For intermediate values of tan (3 and/or larger masses, the cross 
sections are too small for these processes to be useful. For instance, for tan J3 = 10, the cross section is 
below the level of a few femtobarn for m#± ~ 250 GeV. The other mechanisms for if ^ production at 
hadron colliders give even smaller cross sections [206, 224]. 



In a recent paper, Ref. [223], it has been shown that there is a potentially large source of the H ± 
bosons of the MSSM at the LHC: the cascade decays of squarks and gluinos, which are abundantly 
produced in pp collisions, thanks to their strong interactions. These squarks and gluinos can decay into 
the heavy charginos and neutralinos and if enough phase space is available, the latter particles could then 
decay into the lighter charginos/neutralinos, and H bosons, with substantial branching ratios. 

In this section, we summarize the production of H particles through these cascade decays at the 
LHC and describe a Monte Carlo simulation which shows that these final states can be possibly detected 
in some regions of the MSSM parameter space. 



6.2 H ± Bosons from Cascade Decays in the MSSM 

At the LHC the total squark and gluino production cross section is a(q + g) ~ 110 (3) pb for m- g ~ 
rrig ~ 0.5 (1) TeV leading to a large, ~ 3 • 10 7 to 10 6 , number of events with an accumulated luminosity 
of C ~ 300 fb -1 . These squarks and gluinos can decay into the heavy charginos and neutralinos, xt > X3 
and X4 with significant branching fractions, a few ten percent. If enough phase space is available, the 
latter particles could then decay into the lighter charginos/neutralinos, xf > X? an d X2> an ^ bosons, 
with branching ratios of the order of a few ten percent, again. A key point here, is that the coupling of 
the Higgs bosons to chargino and neutralino states is maximal for higgsino-gaugino mixed states Q229[], 
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while the gauge boson couplings to neutralinos are maximal for higgsino-like states. In the gaugino- 
like or higgsino-like regions, this results into the dominance of the decays of the heavier charginos and 
neutralinos into the lighter ones and Higgs bosons compared to the same decays with gauge boson final 
states. 

The total number of charged Higgs particles produced at the end of the chain 

PP -> 99, qq, Q9 -> xf, xi xl + X^ X f, xl, Xi + ^+X (14) 

could be rather large (of the order of a few 10.000 to a few 100.000 events for the high-luminosity option) 
in favorable regions of the parameter space. The interesting and important point to note is that the rate 
of production does not depend very crucially on tan f3 unlike in the other mechanisms mentioned 
above — Note also that bosons could be searched for, if kinematically possible, in the direct decays 
of heavy third generation squarks into their lighter partners or in direct gluino three-body decays. This is 
illustrated below, in two scenarii with the intermediate value tan j3 = 10 and where the universality of the 
gaugino masses has been assumed at the GUT scale, leading to the relation Mi ~ 2M\ ~ M3/3 ~ m g / 3 
at the weak scale. 



c 

q 
O 

Q 
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M 2 = m 5 /3 (GeV) 

Fig. 13: Cross sections times branching ratios for gluinos decaying into squarks and squarks decaying through cascades into 
H bosons as functions of A/2 with fi = 150 GeV and tan (5 — 10. rriq = 900 GeV and nig — 3M2. 




Scenario 1 (Figure [l3|): Here gluinos (with m g = 3M2) are heavier than squarks (nig = 900 
GeV) and therefore g — » qq occurs 100% of the time. The higgsino mass parameter has been chosen 
to be small, [i = 150 GeV, so that all squarks — in particular those of the first two generations — will 
mainly decay into the heavier charginos and neutralinos which are gaugino like with masses rn x + ~ 
m x o ~ 2m x o ~ Mi. For large enough Mi, there is then enough phase space for the decay of the 
heavier gauginos into the lighter higgsino states, with masses m + ~ m o ~ m v o ~ \u\, and 

bosons to occur. For small Mi values, the states \% 4 an d xt are not heavy enough for the decays into 
bosons to occur, in particular for large m H ±. When these decays are allowed, a x BR(— > H^) 
values of the order of 1 pb for m#± ~ 180 GeV and 0.3 pb for m H ± ~ 300 GeV can be reached. 
For increasing values of M2, the gluino mass increases and the cross sections for associated squark and 
gluino and gluino pair production drop and a x BR(^ H^) decreases accordingly; at some stage, only 
the cross section for squark production survives ra, being fixed. The decrease of a x BR(— > H^) with 
increasing M2 is also due to the more suppressed phase space for q — > q'xf, QXa since for large M2, 
m x o,m x ± ~ M2. For even larger M2 values, M2 ~ 650 GeV, the channel X3 ~^ H^X\ opens up, and 
since the phase space is more favorable, because m x o ~ M2/2, a x BR(^ H ) increases again. 

Scenario 2 (Figure |14|): Here the scalar partners of light quarks are heavier than stops which are 
heavier than gluinos, nig = 1.2m g with m g ~ 3M2 = 800 GeV. Gluinos will decay mainly into three- 
body final states mediated by the exchange of top squarks which have a smaller virtuality. Note that the 
cross sections for squark and gluino production are constant and the variation of a x BR(i/ ± ) is only due 
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Fig. 14: Cross sections times branching ratios for squarks decaying into gluinos with the gluinos decaying through cascades 



into H bosons, as as functions of fi for tan (3 — 10. nig 
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to the variation of the branching ratios BR(<7 — » ^QQi X2QQ') an d BR(^3 4, xt ~> xfH^, Xi 2^^)- 
One sees that a x BR(i7 ± ) is relatively large for small values of /j and m H ±, when the gaugino-like 
heavy x states are light enough for the decays xi - > xfH T and xf ~^ Xi 2-^ ± to occur - m the mixed 
region, ~ M2, the mass difference between the heavy and light x states are too small to allow for 
decays in bosons. For large values of /x, a x BR(^f ± ) increases to reach values of the order of ~ 0.1 
pb for m H ± ~ 200 GeV (in particular when the additional channels X3 — ► X\H^ open up) before it 
drops out because of the gradually closing phase space for the decays g — > qqx?, 41 qq'xt- 

Thus, large samples of bosons can be produced in these SUSY cascade decays. In the regions 
of m,}j± and tan (3 values that we are interested in, the dominant decay modes of the charged Higgs 
bosons are @[23(j] H + -> tb (~ 90% sufficiently above the tb threshold) followed by H + r 



(~ 10%). In the simulation which we present below, we will focus on the latter decays which are easier 
to detect in the jet environment of the LHC. 

6.3 An Event Generator Analysis 

Selecting the one-prong hadronic decays of these tau leptons will allow us to exploit the tau polarization 
effects in our analysis. Therefore, the signature of bosons produced in SUSY cascade decays consists 
of one hard r-jet plus additional hard jets (often 6-jets) accompanied by a large amount of missing energy 
due to the presence of the lightest neutralinos and the neutrinos. 

The main SM processes leading to the same signature are top pair production, pp — > ti, with one 
top decaying hadronically and the other one decaying leptonically, and QCD + jets production with 
— > t^v t . As we will argue below, these SM backgrounds can be efficiently suppressed. However, 
a more difficult task will be to distinguish the signal from other SUSY cascade decay processes. As 
was discussed above, all squark and gluino production in colliders will end up in lightest neutralinos and 
fermions through cascade decays via gauginos and sleptons. Events in which these fermions are taus will 
mimic the signal. Hence, a good understanding of the nature of these SUSY backgrounds will be needed 
prior to a search for charged Higgs bosons in cascade decays. 

We performed a Monte Carlo simulation of the production signal and the main backgrounds for the 
following scenario: m H ± = 200 GeV, tan (3 = 10, n = 450 GeV, M 2 = 2M X = 200 GeV, m-g = 800 
GeV, rrig = l.2m g and m ; - = 300 GeV The signal and background events were generated with PYTHIA 



6.152 [10C]. To account for the tau polarization, PYTHIA was interfaced with the TAUOLA [45] pack 



age. The detector aspects were simulated using CMS JET 4.801 []135j], which contains a parameterized 
description of the CMS detector response. The effects of pile-up at high luminosity running of the LHC 
have not been included. The features that will allow to distinguish the signal from the SM and SUSY 
backgrounds are summarized below. 
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Fig. 15: Normalized distribution of E™ lss (left) and normalized Et distribution of the hardest jet in the event 
(right) for the ti background (green-light dashed line), the SUSY cascade background (blue-medium line) and the 
charged Higgs signal (red-dark line). 



In the left-hand panel of Figure 15, the E^ 1BB distribution, normalized to the number of events, 
is shown for the ti background, the SUSY cascade background and the charged Higgs boson signal. 
Demanding a very large E™ ss in the events allows us to eliminate the ti background. Also part of the 
SUSY cascade background can be suppressed relative to the signal due to a slight excess of missing 
energy in cascade decays including charged Higgs bosons. Similarly, making a hard requirement on 
the Et of the hardest jet in the event will help in the background rejection. In the right-hand panel of 
Fig. 15, this Et distribution is shown and it is clear that the SM processes can be efficiently suppressed. 

After eliminating the SM backgrounds, the more difficult task remains to discriminate the 
bosons from the other particles in the SUSY cascade decays, r-leptons in the cascade backgrounds 
originate from charginos and neutralinos in the intermediate state. In our scenario, charginos decay 
predominantly into bosons and neutralinos into Z° bosons. Taus coming from these particles will 
typically have a softer spectrum than the ones coming from a 200 GeV charged Higgs boson. Therefore 
we impose a lower limit on the transverse energy of the r-jet, E T ~ }et > 120 GeV, which is well above 
the mass of the and Z° bosons. As an upper limit on E T ~ iet , we chose the mass of the charged 
Higgs boson. 

The significant presence of W s in the SUSY cascade background also needs to exploit the tau 
polarization effects [191] explained in section ^]: by selecting events in which the fraction of the r-jet 
transverse momentum carried by the charged pion is large, the SUSY background involving W bosons 
can be suppressed relative to the signal. The tau leptons in the background coming from either f and Z° 
decays or MSSM neutral Higgs boson (in particular H and A) decays cannot be suppressed this way. 

These considerations lead us to the following selection criteria in order to distinguish between 
signal and background events: 

i) The transverse missing energy E™ 1 ^ in the event should be larger than 300 GeV. 

ii) The hardest jet in the event should have Et > 400 GeV. 
in) Events with more than five jets are rejected. 

iv) We require exactly one hadronically decaying tau (1 -prong) i.e. we demand a narrow jet within 
|?7| < 2.5 which should contain a hard charged track with pt > 5 GeV in a cone of AR = 0.15 radians 
around the calorimeter jet axis, and it should be isolated i.e. no charged tracks with pt > 2 GeV are 
allowed in a cone of AR = 0.4 radians around the axis. 
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v) The Et of the r-jet, defined as the Et reconstructed in a cone of AR = 0.4 radians around the 
jet axis, should be in the interval: 120 GeV < E T ~ jet < 200 GeV. 

vi) More than 80% of the r-jet transverse energy should be carried by the charged track: see the 
requirement (||) in section ^ 

Events that satisfy conditions (i) and (ii) can be efficiently triggered on using the jet and missing 
energy triggers. If the above selection criteria are applied, we obtain the following numbers for an 
integrated luminosity of 100 fb -1 : 17 events for the boson signal and 21 events for the SUSY 
cascade background, while the ti background is completely suppressed. This results in a significance 
S/\f~B of 3.5, meaning that a 5<r-observation of bosons produced in SUSY cascade decays could be 
made with about 2 years of high luminosity data of the LHC, assuming the above physics scenario and 
provided the SUSY background processes are well understood. 



6.4 Conclusions 

We have shown that charged Higgs boson production from cascade decays of strongly interacting SUSY 
particles can occur with large rates, in favorable domains of the MSSM parameter space. This is true, 
in particular for intermediate values of tan f3 where the standard production processes are disfavored 
because of the smallness of the H~tb Yukawa coupling. We have shown that the SM background to 
the cascade production can be efficiently suppressed. By exploiting the characteristics of the r-jet in the 
final state, the H signal can be made visible above the other SUSY cascade decays. 

Our strategy to distinguish the bosons from other cascade processes depends on prior knowl- 
edge of the properties of the SUSY background. Therefore, before any reasonable search for charged 
Higgs bosons in cascade decays can take place, the general nature of SUSY cascade decays should be suf- 
ficiently understood in order to select the kinematical region in which the signal is enhanced with respect 
to the SUSY background. Due to the complexity of the signature and the dependence on calorimetric 
energy resolutions, these results would need to be confirmed by a more detailed simulation. However, 
our current findings point towards the conclusion that there exists a potential for observing charged Higgs 
bosons produced in SUSY cascade decays at the LHC, provided the nature of the SUSY background is 
well understood. 

A similar analysis, dealing with the production of the neutral Higgs particles of the MSSM through 
the cascade decays of squarks and gluinos is under way [231]. 



7. Summary 

We have investigated the feasibility of detecting various signatures of the charged Higgs boson in order 
to provide a complete coverage for searches in upcoming experiments at the Tevatron and at the 
LHC. 

In the threshold region, i.e., for m H ± ~ m t , the gg — ► tbH^ process has been used to correctly 
account for the production and decay phenomenology in this region of parameter space instead of 
the usual narrow width approximation. It is found that a significant charged Higgs boson signal can 
be extracted in the channel — > t^v t after an integrated luminosity of 15 fb -1 at the Tevatron 
Run 2 while the factorization approach used in the narrow with approximation would not predict such a 
favorable signal. The case of the LHC, which is further complicated by the potential problem of double 
counting when calculating the inclusive cross section, is currently being studied. 

In models with singlet neutrinos in large extra dimensions, the process H~ — ► t£ tp+c.c. — which 
is completely suppressed in the 2HDM — can be enhanced thanks to the large number of Kaluza- Klein 
states of the right handed bulk neutrino. Such a signal can be observed at the LHC with significances 
exceeding 5-a. However, in order to identify the scenario that is realized (2HDM-II or large extra dimen- 
sions) both the transverse charged Higgs mass and the fraction of the energy carried away by charged 
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tracks in the one-prong r decays must be reconstructed. Further evidence for large extra dimensions 
would come the measurement of the polarization asymmetry. 

In the intermediate tan f3 region, the charged Higgs boson decays to SM particles do not yield 
any significant discovery potential. To cover this region of the parameter space, charged Higgs decays 
into chargino-neutralino pairs have been studied. It is demonstrated that, by searching for a three-lepton 
final state with a top quark and a large missing energy resulting from the gg — ► tbH^ where the charged 
Higgs decays to xfx® or xfx*s> a significant charged Higgs signal could be extracted for intermediate 
tan/3 values and m H ± ~ 300 GeV. Further analysis is in progress aiming to extend the coverage to 
higher charged Higgs boson masses. 

A significant source of charged Higgs boson production is the cascade decays of SUSY particles 
(squarks and gluinos) also sensitive to the intermediate tan (5 region where, as previously mentioned, the 
SM charged Higgs boson productions and decays yield no discovery potential. The subsequent decay of 
the charged Higgs boson into the r-lepton has been studied taking advantage of the polarization effects 
in the T-jet final state. It is determined that a significant charged Higgs boson signal can be observed 
through these cascade decays of SUSY particles provided the SUSY background is well understood. 

The charged Higgs boson mass can be determined in — > tb and — > t^v t where the 
precisions range from 0.5% at ~ 200 GeV to 1.8% at ~ 500 GeV for an integrated luminosity of 300 fb -1 . 
By measuring the rate of H ± — > r ± ^ T , tan/? can be determined with precisions ranging from 7.4% at 
tan/3 = 20 to 5.4% at tan/3 = 50 for an integrated luminosity of 300 fb^ 1 and assuming a 10% 
uncertainty on the luminosity. 

The studies discussed here are the necessary continuation of previous work and in the process, 
reveal the enormous potential of the charged Higgs boson: 

• To understand the structure of the Higgs sector through the determination of m H ± and tan/3 and 
to probe the decoupling limit of the MSSM, thus distinguishing between the SM and the MSSM, 
particularly in the — ► r ± ^ r channel. Indeed, the scope of the parameter space covered by this 
channel is comparable to the reach of the A/ H — > tt channel in the neutral Higgs sector. 

• To provide evidence — or not — of large extra dimensions, the existence of which, especially at 
the electroweak scale, would constitute the solution to the outstanding hierarchy problem. 

• To explore the SUSY particle arena and thus providing various signatures whose detection will 
constitute evidence of the existence of these exotic particles. 
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